Measurement of the charge asymmetry 
in top quark pair production in pp coUision data 
at = 7 TeV using the ATLAS detector 



Dissertation 



zur Erlangung des mathematisch-naturwissenschaftlichen Doktorgrades 

"Doctor rerum naturalium" 
der Georg-August-Universitat Gottingen 



vorgelegt von 
Fabian Kohn 
aus Northeim 



Gottingen, 2012 



Referent: Prof. Dr. Arnulf Quadt 
Korreferent: Dr. Carsten Hensel 



Tag der miindlichen Priifung: 07.03.2012 



Measurement of the charge asymmetry 
in top quark pair production in pp coUision data 
at = 7 TeV using the ATLAS detector 



by 



Fabian Kohn 



A measurement of the charge asymmetry in the production of top quark pairs in 
the semileptonic decay channel has been performed. A dataset corresponding to an 
integrated luminosity of 1.04 fb~^, obtained at a centre-of-mass energy of ^/s = 7TeV 
with the ATLAS experiment at the LHC, was used. After performing a selection of 
events with one isolated lepton, at least four jets and missing transverse energy, a 
kinematic fit was performed to reconstruct the tt event topology. 
The charge asymmetry was determined using the differential distribution of the 
reconstructed observable — where yt and denote the top and antitop quark 
rapidities, respectively. An unfolding procedure was applied to correct for detector 
acceptance and resolution effects and to obtain the corresponding distribution at 
parton level. The total charge asymmetry after unfolding was measured to be 



in agreement with the Standard Model prediction. 

In addition, a simultaneous unfolding in \yt\ — and the invariant tt mass, Mf-^, 
was performed. 
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Not only is the universe stranger than we imagine, 
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1. Introduction 



The field of Elementary Particle Physics is concerned with understanding the most fundamental 
building blocks of Nature and their interactions. The Standard Model of Particle Physics [THS] 
is one of the most successful and thorough theories in physics, its predictions being in aston- 
ishing agreement with the observed phenomena to highest precision. Regardless of its success 
in explaining the interactions of all fundamental particles, the Standard Model is not without 
shortcomings. The fine-tuning of radiative corrections to the Higgs boson mass [7l[H], the strong 
evidence for the existence of Dark Matter and Dark Energy, and a missing mechanism to de- 
scribe gravity within the framework of a quantum field theory are only a few of the remaining 
ambiguities. 

According to current knowledge, there exist six quarks in Nature: the up, down, strange, 
charm, bottom and the top quark along with the six known leptons (electron, muon and tau 
together with their corresponding neutrinos). Quarks and leptons are grouped into three subsets 
or generations of a quark and lepton doublet each. 

The top quark is the heaviest known elementary particle and has been the focus of studies 
for several decades, from indirect searches using electroweak precision data to its discovery by 
the CDF [9] and D0 [lOj experiments at the Tevatron in 1995, the observation of the single top 
quark |lllll2j in 2009, to newest precision measurements of its properties at the Tevatron and 
the LHC. Its unique characteristics, namely the large mass of about 173 GeV |13j and its short 
lifetime, provide the opportunity to perform precise measurements of electroweak interactions. 
Due to the affinity of its mass to the electroweak scale and the potential link to the vacuum 
expectation value of the Higgs field, the top quark also allows indirect constraints of the Higgs 
mass in combination with precision measurements of the W boson mass. Furthermore, since 
the top quark is the only quark which has a decay width larger than the hadronisation scale, 
it does not form hadronic bound states. As a result, top quark properties, such as its spin, are 
accessible without being obscured by the process of hadronisation. 

At hadron colliders, such as the Tevatron or the Large Hadron Collider, top quark pairs are 
mainly produced via the strong interaction, either via gluon-gluon fusion or via quark-antiquark 
annihilation: 

q + q t + t, 

g + g t + t. 

In the Born approximation, these production mechanisms are entirely symmetric under the 
exchange of the final state top and antitop quark. Consequently, there is no angular discrimi- 
nation between the top and antitop quark and the resulting predicted differential distributions 
are identical for both particles. 

In the Standard Model, an asymmetry in the production of top quark pairs arises due to 
radiative corrections from virtual and real gluon emission if higher order corrections are taken 
into account. These higher order corrections introduce interferences between amplitudes which 
are odd under the exchange of the final state quark and antiquark. Interference terms of final 
state and initial state gluon bremsstrahlung, and of higher order amplitudes with Born level 
amplitudes contribute to an overall imbalance of the differential distributions of the final state 
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top quark and antitop quark. An additional small contribution originates from the interference 
of different amplitudes in quark-gluon scattering: 

g + q t + t + q'. 

The measurement of the charge asymmetry in top quark pair production provides the oppor- 
tunity to verify perturbative Quantum Chromodynamics and consequently, the Standard Model. 
Moreover, a similar effect is predicted and observed in Quantum Electrodynamics |14ftl6j . where 
radiative corrections lead to an asymmetry in the electroweak production of fermion-antifermion 
pairs. Since this effect has been studied to very high precision, the verification of its counterpart 
in Quantum Chromodynamics would be yet another confirmation of the Standard Model and 
its predictions. 

Furthermore, potential new physics, in particular theories involving the breaking of elec- 
troweak symmetry, could lead to deviations from the Standard Model expectation due to large 
anomalous couplings to the top quark predicted in numerous theoretical models. As a matter 
of fact, current independent measurements performed by CDF [T71[18] and D0 [191120] suggest a 
possible discrepancy between the predicted and observed charge asymmetry in proton-antiproton 
collisions. This effect is observed in particular for high invariant tt masses and high tt rapidity 
differences, which is supported by several models predicting physics beyond the Standard Model. 

Since the predicted charge asymmetry is small at hadron colliders due to the probabilistic 
nature of the initial state parton kinematics, precise knowledge of the detection mechanisms, 
sophisticated analysis methods and detailed understanding of potential systematic effects are 
crucial to accomplish such a measurement. This is in particular true for the Large Hadron 
Collider, where a high centre-of-mass energy and a symmetric hadronic initial state (pp) make 
this measurement even more difficult. The increased fraction of top quark pairs produced via 
(charge symmetric) gluon-gluon fusion lead to a dilution of the measured asymmetry. In addi- 
tion, there is no preferred initial state quark direction in proton-proton collisions and hence no 
resulting forward-backward asymmetry which could be measured directly as it is the case at the 
Tevatron. Consequently, a new analysis concept and new observables have to be considered to 
perform this measurement under the conditions of the LHC. 

The charge asymmetry has been measured in top quark pair production at the Tevatron by 
both the CDF [T7] and D0 [19] collaborations and preliminary results have also been shown 
by the CMS [21] experiment at the LHC. This thesis describes a measurement of the charge 
asymmetry in top quark pair production which has been performed with the ATLAS experiment 
for the first time in 2011 [22] . 

This document is organised as follows: Chapter [2] gives an introduction to the theoretical 
aspects of the Standard Model and theories beyond, paying special attention to the top quark 
and the charge asymmetry in top quark pair production. Chapter [3] covers technical aspects 
of the LHC and the ATLAS detector. Definition and description of the objects taken into 
account for the described analysis and the trigger strategy is given in Chapter H] followed by 
a summary of the event selection performed to increase the fraction of relevant signal events 
with respect to various background processes in Chapter [SJ A summary of the data and Monte 
Carlo samples used, and a description of data driven methods to determine the background 
contributions from Vl^-|-jets and QCD multijets, are given in Chapter A detailed explanation 
of the reconstruction method used to obtain parton level information of the top and antitop quark 
based on measured quantities follows in Chapter [T] An unfolding approach performed to account 
for detector acceptance and resolution is described in Chapter [HI followed by Chapter [9l covering 
relevant systematic uncertainties affecting the analysis. Finally, the results of the analysis are 
presented in Chapter [10] and a summary of this thesis is given in Chapter [TTl 
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2. The Top Quark Charge Asymmetry in 
Context of the Standard Model 



2.1. The Standard Model 

In Nature, all observed matter consists of building blocks which are considered to be elementary, 
the leptons and quarks, shown schematically in Figure 12. 1[ They are classified into three gener- 
ations or families with increasing order of quark masses. Leptons and quarks obey Fermi-Dirac 

^ e (l? @ Higgs 

Figure 2.1.: Summary of elementary particles. Aside from the quarks (upper left box) and 
leptons (lower left box), the gauge bosons (right vertical box) and the hypothetical 
Higgs boson are shown. 

statistics and hence carry a non-integer spin. Leptons are described by their quantum numbers, 
electric charge (Q), the third component of the weak isospin (T3), and the weak hypercharge 
(Yw = 2((5 — Ts)), summarised in Table [2m 





Leptons (spin 






Flavour 


Mass [MeV] 


Q 


T3 


Yw 




< 2 • 10-3 





1 
2 


-1 


e 


0.511 


-1 


1 
2 


-1 




< 0.19 





1 

2 


-1 




105.658 


-1 


1 
2 


-1 


Vr 


< 18.2 





1 
2 


-1 


T 


1776.82 ±0.16 


-1 


1 
2 


-1 



Table 2.1.: Leptons and their properties and quantum numbers, ordered by generation |23j . 

Where no uncertainty on the mass is given, it is negligible at the given precision. 

Similarly, quarks are assigned the quantum numbers charge (Q), the third component of the 
weak isospin (T3), and hypercharge {Y = 2(Q — T3)), describing qualities of related hadronic 
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Quarks (spin s - 


2-' 








Flavour 


Mass [MeV] 


Q 


Ta 


1" 


u 


o.u to ( .u 


2 
3 


1 

2 




4 
3 


d 


1.5 to 3.0 


1 

3 


1 
2 




1 

3 


c 


1250 ± 90 


2 
3 


1 
2 




4 
3 


s 


95± 25 


1 

3 


1 
2 




2 
3 


t 


(173.2 ±0.9) • 10^ [13] 


2 

3 


1 
2 




4 
3 


b 


(4.20 ±0.07) • 10^ 


1 

3 


1 
2 




2 

^3 



Table 2.2.: Quarks and their properties and quantum numbers, ordered by generation |23) . 



bound states. A summary of these properties can be found in Table 12.21 Furthermore, all quarks 
carry a colour charge, denoted by a respective quantum number which can take the values red, 
green and blue. 

The first generation is constituted by the up (n) and down {d) quark doublet (the build- 
ing blocks of protons and neutrons) alongside the electron (e) and the electron-neutrino (fg)- 
The second generation contains the charm (c), strange (s) quarks, the muon (/i) and muon- 
neutrino {v^). Finally, the top {t) and bottom (6) quark compose the third generation, together 
with the tau (r) and the tau-neutrino {vr) in the lepton sector. The top quark assumes a quite 
distinct role among the other quarks due to its large mass of (173.2 ± 0.9) GeV [T^ . 

Within the Standard Model, particle interactions are described by a quantum field theory 
consistent with both quantum mechanics and special relativity, combining the electroweak the- 
ory and Quantum Chromodynamics into a structure denoted by the gauge symmetry group 
SU(3)(^ X SU(2)^ X U(l)y. This structure describes colour charge (C), weak isospin (L) and 
hypercharge (Y) gauge groups. The underlying gauge theory is non-Abelian due to the non- 
commutative nature of the SU(3) and SU(2) field strength tensors. 

Given the matter and gauge fields, and requiring local gauge invariance and renormalisability, 
the Standard Model Lagrangian can be constructed as 

-CSM = -CsU(3) + ^SU{2)xU(l) (2-1) 

^Gauge . ^Matter , ^Gauge . ^Matter , ^Higgs . ^Yukawa /n n\ 

~ ''-SU(3) +^SU(3) + ^SU(2)xU{l) + ''-SU{2)xU{l) + ^SU{2)xU{l) + ''"SU{2)xU{l)- 

The first term, describing strong gauge interactions, is given by 

^Gaugc _ _}_rpy. r'i^'^ n ('9 'X\ 

^SU(3) ~ 

where G^'^ is the gluon field strength tensor and gs is the strong gauge coupling constant. Strong 

interactions are described by perturbative Quantum Chromodynamics (QCD). 

The matter term of the SU(3) Lagrangian contains the gauge covariant derivatives of the quarks: 

/iMatter •- rf.a a\ 
^SV{3) =^^i«-P/3^> (2-4) 

where a,/3 G [1,2,3] are the quark colour indices. A summation over the quark flavour index i 
is implied and it is 

^-^ = 9^5^ + ^550^^. (2.5) 

The first term of the electroweak Lagrangian describes the corresponding gauge interactions 
of the electroweak theory: 

^Gauge — J_TrPF^''PF —B'^'^B (2 6) 

^SU(2)xU(l) ~ Ag'"^ ^ ' 



4 



2.1. The Standard Model 



with the weak isospin and hypercharge gauge field strength tensors W^'^ and B^'^ , respectively. 

The SU(2) and U(l) gauge couplings are represented by the constants g and g' . Through 
mixing of the B and W3 fields, the photon and Z boson are generated: 

^\ = { sin6lM/\ f B\ , . 

zj \-smew cosewjywsj' ^ ' 

where 6\y denotes the weak mixing angle. Similarly, the bosons are generated through 
mixing of the Wi and W2 fields. 

The matter term of the electroweak Lagrangian contains the kinetic energy terms from the 
fermions and their gauge field interactions: 

^^uSxUd) = ^^l^^l + i^R^^R + ^^R^d^R + i^i^^i + i^R^^R- (2-8) 

A summation over the index i is implied. The indices L and R refer to the left and right-handed 
chiral projections 

V'L = (1 - 75)^ and Vfl = (1 +75)-. 
The left-handed quark and lepton fields are represented by the SU(2) doublets 

^i={1,) and li=(''U (2.C 



L 



while the right-handed fields are represented by the singlets u^, and e^. The gauge covariant 
terms, describing the electroweak gauge interactions of the fermions, are given by 

Df^qi = {d^+'-^TW^ + i^B^)ql, 
DA = id>. + jrW^-i^B^)lh 

D^u'n = {d^ + i'^g'B^)u'R, (2.10) 

D.e^R = {d^-ig'B,)e}R. 

Note that there are no mass terms for the fermions in either the SU(3) or the SU(2) x U(l) 
gauge theory since such terms are forbidden by the gauge invariance of the Standard Model. 
A Dirac mass term for a fermion field is not invariant under a chiral transformation and hence 
would violate the requirement of gauge invariance and renormalisability. However, since weak 
interactions are observed to be short ranged, the gauge bosons must obtain non-vanishing masses 
through a different mechanism. Both gauge boson and fermion masses are generated by spon- 
taneous symmetry breaking of the Higgs term in the Standard Model Lagrangian, 

4i?2)xu{i) = {D^<P)^D,4> + A^VV - K<t^H?, (2.11) 
containing the kinetic energy of the Higgs field, which is represented by the complex scalar field 

(2.12) 
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Figure 2. 2. : The Higgs potential^. Elcctrowcak symmetry breaking induces a non-zero vacuum 
expectation value in the minima of the Higgs potential, leading to effective masses 
for the gauge bosons and fermions. 



The first term describes the Higgs field interactions with the gauge fields and the latter two 
denote the Higgs potential, shown in Figure 12.21 

The final term in the Standard Model Lagrangian describes the Yukawa interactions of the 
fermions with the Higgs field: 

^vt^mi) = -^'dlie^y'R - r'Jqiecp4 - Ti^lie^e^ + h.c, (2.13) 

where e = ia2 denotes the two dimensional total antisymmetric tensor which ensures electrical 
neutrality of the individual Yukawa terms and a2 represents the Pauli matrix 

Furthermore, the Yukawa couplings T^, F^ and Fe denote the respective complex 3x3 matrices in 
generation space, describing the interactions between the Higgs doublet and the different fermion 
flavours. As they are not required to be diagonal, a mixing amongst the three generations is 
allowed. 

The electroweak symmetry is spontaneously broken by acquisition of a non-zero vacuum ex- 
pectation value of the neutral Higgs fleld component 

^ - " (2.15) 



2A V2 

which, consequently, generates masses Mz and M^y for the electroweak gauge bosons through 
unitarity gauge, and Mf for the fermions from gauge invariant Yukawa couplings Fj of the Higgs 
boson to fermions: 

Mz = \vg, (2.16) 
Mw = Ivy^JTT^, (2.17) 
Mf = Tf^. (2.18) 



A summary of the gauge bosons and their properties can be found in Table 12.3 



^ Image taken from wikipedia (public domain). 
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Bosons 


(integer spin) 






Particle 


Mass [GeV] 


Charge 


Spin 


(7) 


Electromagnetic force 





n 
u 


1 
1 


(3) 


Strong force 








1 


{W-) 






80.403 ±0.029 


-1 


1 


{W+) 


Weak force 




80.403 ±0.029 


+1 


1 








91.188 ±0.003 





1 




Mass (hypoth.) 


116 


- 127(95% c.L.) mm] 









Table 2. 3. : Gauge bosons and their properties and quantum numbers [13] . 



2.2. The Top Quark 

In the following, the production and decay of top quark pairs and singly produced top quarks 
{single tops) within the Standard Model will be discussed in detail. Furthermore, an overview 
of important properties of the top quark and their measurement will be given. In particular, 
the charge asymmetry in the production of top quarks pairs within the Standard Model and in 
theories beyond will be covered. 

2.2.1. Top Quark Production at Hadron Colliders 

At hadron colliders, tt pairs are mainly produced through strong interactions described by per- 
turbative QCD. Interactions between the quark and gluon constituents of the colliding hadrons 
(either protons or antiprotons) participate in a hard scattering process and produce a top quark 
and an antitop quark in the final state. At Born level approximation, top quark pairs can be 
produced via gluon-gluon fusion (gg) or via the annihilation of quark-antiquark pairs (qq). The 
relevant leading order Feynman diagrams for the contributing processes are shown in Figure 12.31 



5 '^"O'O'O'O'ffY 



g > 0.0 0.0 0. 




Figure 2. 3. : Lowest order diagrams contributing to top quark pair production at hadron col- 
liders. Top quarks are produced via strong interaction, either in quark-antiquark 
annihilation (top) or gluon-gluon fusion (bottom). 

Due to the fact that hadrons are composite particles, consisting of partons with unknown 
fractions x of the initial hadron momenta, the initial state of the parton interaction is not 
precisely known. However, hadron interactions in pp and pp collisions can be described by 
separating the partonic reactions into a short distance and a long distance contribution. 
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The long distance part can be factorised into longitudinal parton momentum distribution 
functions (PDFs) where fip denotes an (arbitrary) factorisation scale describing the 

separation of the long and short distance contributions. An additional renormalisation scale is 
introduced to account for higher order corrections, where ultraviolet divergent terms may emerge 
and a renormalisation approach can be used to absorb such divergences into corresponding 
counter terms. Both scales /i^ and //^ are commonly chosen to correspond to the momentum 
transfer fip = fi'j^ = Q^. Furthermore, for the calculation and simulation of top quark processes, 
is typically chosen such that hf = I^R = Q = fnt corresponds to the top pole mass mt and 
the associated scale variation dependency is studied. 

The PDFs represent the probability distribution of observing a parton of type i at a given 
scale with a longitudinal parton momentum fraction Xj. Since these probabilities cannot be 
universally derived from QCD, they have to be provided from experimental studies of the proton 
structure, mostly from deep inelastic lepton-proton scattering experiments at the HI [26fl29j and 
ZEUS |3UH33| experiments at the HERA electron-proton collider. As an example, the e'^p and 
e~p production cross-sections measured in deep-inelastic scattering experiments at HERA can 
be found in Figure [231 [M] in comparison to the CTEQIO PDF next-to-leading order (NLO) 
prediction [35] . 




Shifted HERA-1 data (circles) 
CTIO NLO tlieory (lines) 



10 



10^ lO"' 

Q- [GeV-; 

(a) e'^p 



10' 



10" 



lO'* 



10 ^^ h 



10" 



T 

D- 10 



1 - 



10' 



10"- 



10' 



x=2- 10". i=l2 
■ ■ - x=5- 10";^ 



x=8- 10", i=10 

x=1.3- 10": i=9 
, B x=2- 10": i=8 

9 9 x=3.2- 10": i=7 



1 » «- 



=5- 10": i=6 
_ x=8- 10": [=5 



M x=0.13, i=4 
__^x=0.18,,=3 



Shifted HERA-1 data (circles) 
CTIO NLO theory (lines) 



10" 



lo-' 



10' 



10-' 



Q- [GeV"] 

(b) e'p 



Figure 2.4-: Comparison of CTEQIO NLO predictions [35] for reduced cross-sections in e+p 
(left) and e~p (right) neutral-current deep inelastic scattering experiments from 
combined HERA-1 data [34], with correlated systematic shifts included. 



The short distance term arises from the hard scattering process of the respective partons, 
denoted by the partonic cross-section for partons i and j, aij. This contribution is characterised 
by high momentum transfer. Hence, it is not dependent on the incoming hadron type or the 
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respective wave functions and can be described by perturbative QCD, as indicated by the leading 
order diagrams in Figure [2T3l 

At a given centre-of-mass energy ^/s and for a top mass parameter mj, the total top quark 
pair production cross-section can be calculated from the short distance and long distance terms 
as 

o"tf (\/s,"it) = dxidxjfi{xi,Q'^)fj{xj,Q^) X aij_,ti{p,ml,Xi,Xj,as{Q'^),Q^) , (2.19) 

where the summation is performed over all permutations of i,j = {q,q,g}. The PDFs of the 
initial state protons are denoted by fi{xi,Q'^) and fj{xj,Q'^), respectively and the parameter p 
is given by 



4mf 



7i' 



(2.20) 



■^XiXjS 

where XiXjS = s denotes the effective centre-of-mass energy in the partonic reaction. 

The probability of a parton i to be carrying a momentum fraction of Xi decreases significantly 
with rising Xj, as can be seen in Figure where two PDFs from the CTEQIO PDF set [3S] 
are shown as an example. The PDFs have been evaluated at scales fi = 5 GeV and p = mt, 




(a) = 5 GeV 




(b) ^ = mt 



Figure 2.5.: CTEQIO parton distribution functions at different momentum transfers for 
gluons and different quark/antiquark flavours. Shown are the PDF sets for p = 
5 GeV (left) and p ^ mt (right). 



respectively, where p = Q. 

The minimal energy carried by the two incoming partons to produce a top quark pair at the 
threshold (i.e. at rest), is given by 

^XiXjS > 2mt (2.21) 



and hence, assuming both partons carrying the same momentum fractions as an approximation, 

(2.22) 



This corresponds to a typical value of x ~ 0.05 at the LHC for a centre-of-mass energy of 
^/s = 7TeV. As shown in Figure [231 the gluon PDFs dominate significantly over any other 
parton in the corresponding range of x. Consequently, the production of top quark pairs at the 
LHC is dominated by gluon-gluon fusion. At the Tevatron (where the typical value of x is of 
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the order of 0.2) the production of top quarks is dominated by quark-antiquark annihilation 
processes, in particular involving up and down valence quarks. Since the centre-of-mass energy 
at the LHC is significantly higher, top quark pairs are typically produced above the threshold, 
but still within the gluon-gluon fusion dominated range of the PDFs. 

The total tt cross-section at the LHC is predicted in an approximate next-to-next-to-leading 
order (NNLO) calculation to be ISStJg pb [36H38] for a centre-of-mass energy of ^/s = 7TeV and 
mt = 172.5 GeV. Preliminary measurements have been performed at both ATLAS and CMS, 
yielding 

attiVs = 7TeV) = 179.0l^;^ (stat.+syst.) ± 6.6 (lumi.) pb [39], (2.23) 
attiVs = 7TeV) = 165.8 ± 2.2 (stat.) ± 10.6 (syst.) ± 7.8 (lumi.) pb 00], (2.24) 

respectively. Both measurements are in agreement with the Standard Model prediction. 

This cross-section is several orders of magnitude lower than, for example, the SM Z and 
W boson production cross-sections or the inclusive QCD multijet production cross-section at 
comparable values of Q^. This can be seen in Figure [21S1 where the total production cross- 
sections for several SM processes are shown as a function of centre-of-mass energy of the colliding 
(anti) protons. Consequently, a sophisticated real-time selection to identify the relevant final 

proton - (anti)proton cross sections 

10' 
10' 
10' 
10' 
10' 
10' 
10' 
10' 

s , 

S ^° 

10 

10"' 
10"' 
10"' 
10"' 
10"' 
10"' 
10"' 

0.1 1 10 

Vs (TeV) 

Figure 2.6.: QCD predictions for hard-scattering cross-sections at the Tevatron and the LHC 
[4T] . The top quark pair production cross-section is denoted by at ■ The disconti- 
nuities in the different curves denote the change from pp collisions at the Tevatron 
to pp collisions at the LHC. For the LHC, different centre-of-mass energies are 
highlighted by three vertical dashed lines, where the leftmost line corresponds to 
a centre-of-mass energy of ^/s = 7 TcV. 

state particles and obtain a good signal to background separation with respect to other SM 
processes and, more importantly, the dominant QCD multijet background, is crucial for all top 
quark related measurements at the LHC. Furthermore, an extensive theoretical understanding 
and modelling of these backgrounds is necessary to facilitate the measurement of top quark 
properties to highest precision and in order to achieve a sensitivity to potential deviations from 
the Standard Model expectations. 
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2.2. The Top Quark 



2.2.2. Production of Single Top Quarks 

In addition to the production of top quark pairs via the strong interaction, single top quarks can 
be produced in electroweak charged current interactions. Three mechanisms for this production 
exist, as shown in Figure 12.71 Single top quarks can emerge in the fusion of W bosons and gluons 
(f-channel process) similar to the production of heavy flavour quarks in deep-inelastic scattering 
via charged current interactions. In addition, they can be produced by annihilation of quark- 
antiquark pairs (s-channel process) and exchange of an off-shell W* , or by Wt production in 
quark-gluon interactions. The corresponding production cross-sections have been approximated 




(a) (b) (c) 

Figure 2. 7. : Electroweak single top production diagrams via ly-gluon fusion (a) , exchange of 
an off-shell W* (b) and Wt production (c) . 



at next-to-next-to-leading-order to be 



t-ch 

t 




7TeV) 


= 64.57l|^}pb 


s-ch 
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(Vi = 


7TeV) 


= 4.63t°;l?pb 


Wt 
t 




7TeV) 


= 15.741}:°^ pb 



for rrit = 172.5 GeV in the t-channel, in the s-channel, and for Wt production, respectively. 
Since in all production channels a top charged current is involved, the respective cross-sections 
behave as 

crt oc \Vtb\'^g'^, 

where \ Vtb\ denotes the relative probability that the top quark decays into a bottom quark via the 
exchange of a W boson. This is described by the Cabbibo-Kobayashi-Maskawa (CKM) matrix 
VcKM [23] which summarises the relative transition probabilities in charged weak interactions. 
The magnitudes of the CKM matrix elements, here denoted by the matrix Mckm? are 



/iKdl \Vus\ \Vub\\ 
McKM = \Vcd\ \Vcs\ Kb\ 

\\Vtd\ \Vts\ \Vtb\J 
/0.97428 ± 0.00015 0.2253 ± 0.0007 
0.2252 ±0.0007 

v 



07345+0-00015 



n nnsfi9+o-ooo26 
0.008b2_o^ooo20 



0.0403' 



"0.0011 
-0.0007 



00347+°-°°°^^ 
0410+°-°°^i 

U.U41U_Q QQQ7 

nqqi52+0-°00030 



(2.25) 



As a consequence of this dependency, measurements of the single top quark production cross- 
section provide an implicit sensitivity to the CKM matrix element Vtb- At the LHC, single tops 
are predominantly produced via the t-channel interaction, followed by Wt production, due to 
the large initial state gluon contribution at the LHC centre-of-mass energy, while the s-channel 
process is suppressed. 
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2. The Top Quark Charge Asymmetry in Context of the Standard Model 



First direct evidence for single top quarks was found in 2006 by the D0 collaboration at the 
Tevatron [U], followed by its observation |45|l46j in 2009 by the D0 and CDF experiments. 
First preliminary measurements of the inclusive single top production cross-section have been 
conducted at the LHC, yielding 

^^^(^/^ = 7TeV) = 90tipb[17! 

and 

cj*-'^^ (Vs = 7TeV) = 83.6 ± 29.8 (stat.+syst.) ± 3.3 (lumi.) pb 08] 

in the t-channel as measured by ATLAS and CMS, respectively. The Wt production cross-section 
has been measured at CMS to be 

* (Vs = 7TeV) = 22+^ (stat.+syst.) pb [39] 

and a limit has been set by a corresponding measurement at ATLAS, corresponding to 

a^* (Vs = 7TeV) < 39pb [50] 

at 95 % C.L. Furthermore, searches for single top quarks produced in s-channel interactions have 
been conducted at ATLAS, limiting the corresponding production cross-section to 

cjf^h ^ 7TeY) < 26.5 pb [SI] 
at 95 % C.L. All measurements are in agreement with Standard Model predictions. 
2.2.3. Top Quark Decay 

Due to the magnitude of the CKM matrix element Vtb being close to unity, the top quark decays 
in almost 100 % of the cases via electroweak charged current interaction into a b quark and W 
boson, which then in turn decays either leptonically into a charged lepton and the corresponding 
(anti) neutrino or hadronically into a quark-antiquark pair. At leading order, the Standard Model 
prediction for the total decay width of the top quark, Fj , is given by 



where Gp denotes the Fermi coupling constant 



2 



Taking into account higher order corrections at next-to-leading order, the total top quark 
decay width becomes 



F. = FV 1 ^ 1 + 2 ^ 



* ^ mj J \ mj 



2as /27r2 5 

1 



Svr V 3 2 



(2.28) 



where terms of order m'^/m'f and (as/7r)Af^/mj have been neglected. At a top mass of 170 GeV 
and as evaluated at the Z scale, this yields an approximate predicted decay width of 

Ft « 1.3 GeV, (2.29) 
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2.2. The Top Quark 



and a corresponding mean lifetime of 

n 0.5 • IQ-^"^ s, (2.30) 

which is significantly lower than the time scale corresponding to the strong hadronisation scale 
-'^QCD ~ 250 MeV. Hence, the top quark decays before being able to form hadronic bound 
states such as the M-quarkonium. Consequently, the top quark spin/polarisation properties are 
preserved in its decay and are transferred to the decay products. 

Since the top quark decays almost exclusively into a h quark and W boson, the resulting final 
state decay channels are well defined and can be separated into three cases, characterised by the 
final state particles: 

• Full hadronic final state (alljets): Both W bosons from the pair further decay into 
quarks, leading to a total amount of six quarks including the h quarks from the initial top 
and antitop decays. 

• Semileptonic final state (lepton+jets): One W boson from the ti pair decays into 
quarks, while the second one decays leptonically, leading to a total of four quarks including 
the h quarks from the initial top and antitop decays, and one charged lepton and its 
corresponding (anti-)neutrino. 

• Dileptonic final state (dilepton): Both W bosons from the ti decay into a charged 
lepton and the corresponding (anti-)neutrino, respectively. In addition, two remaining b 
quarks from the top and antitop decays are produced. 

If the charged lepton is a tau in the semileptonic or dileptonic channel, either a muon or electron 
and the corresponding (anti-)neutrino, or further quarks from the hadronic decay of the tau 
lepton are produced. 

The respective W branching ratios (BR) at leading order (LO) can be found in Table [27^ 
Taking these branching fractions into account, possible ti final states and their approximate 



Decay mode 


LO BR 


Measured BR 


-> qq 


6 
9 


(67.60 ±0.27)% 




1 

9 


(10.75 ±0.13)% 




1 

9 


(10.57 ±0.15)% 




1 

9 


(11.25 ±0.20)% 



Table 2.4-'- Theoretical (LO) and measured W branching ratios |23) . 

relative probabilities are shown schematically in Figure 12.81 
2.2.4. Top Quark Properties 

Several properties of the top quark have been studied in collider experiments such as the Teva- 
tron and the LHC. Amongst them, the top quark mass has been determined with a relative 
uncertainty of only 0.5 % [13] by combining the most recent measurements from D0 and CDF. 
This combination constitutes the most precise (in relative terms) mass measurement of any 
quark so far. Since the top quark does not form hadronic bound states, the top quark mass mt 
is defined as the pole mass in this context and is measured to be 

mt = 173.18 ± 0.56 (stat.) ± 0.75 (syst.). (2.31) 
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2. The Top Quark Charge Asymmetry in Context of the Standard Model 



Top Pair Decay Channels 
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Figure 2.8.: Top quark pair decay channels (left) and branching fractions (right jl. 



In addition, measurements to exclude an exotic top quark carrying a charge of which 
is predicted to be 2e/3 in the Standard Model, have been performed. This model has been 
excluded with approximately 95% C.L. ^21 and 90% C.L. [SH], respectively, in independent 
measurements at CDF and D0. A corresponding measurement at ATLAS excludes a top quark 
charge of 4e/3 at more than five standard deviations |54| . 

Since the top quark does not form bound hadronic states due to its small lifetime, it provides 
the unique opportunity to measure quark properties which are usually concealed by hadroni- 
sation, such as spin correlations between quark-antiquark pairs produced at hadron colliders. 
Since all other quarks depolarise due to QCD interactions before fragmentation, it is not possible 
to gain knowledge about the spin from the final state, while the final state particles in the top 
quark decay preserve significant amount of information about the spins of the initial top and 
antitop quarks to potentially allow their measurement. First studies at the Tevatron using the 
angular distributions of the final state particles indicate a correlation strength C of 

C = 0.57 ± 0.31 (stat.+syst.) (SS] (2.32) 

in the beam basis, compatible with the next-to-leading order Standard Model prediction of 
CsM = 0.781qq4 [56] and excluding a non-correlation hypothesis at 97.7% C.L. Similar studies 
at ATLAS indicate a correlation strength ^heiicity hi the helicity basis of 

Aoiicity = 0.34;0;}f (stat.+syst.) [57], (2.33) 

which is compatible with the corresponding next-to-leading order Standard Model prediction. 

Furthermore, polarisations of W bosons from the top quark decay are predicted to arise due 
to the different possible helicity states of the produced on-shell W bosons. Consequently, the 
W helicity has been measured and the obtained left-handed, longitudinal and right-handed 
polarisations have been found to be consistent with the Standard Model predictions [581159]. 

The charge asymmetry in top quark pair production is key topic of this thesis and will be 
discussed in detail in the following section. 

^Here, e denotes the electron charge. 

^Images taken from http://www-dO.fnal.gov/Run2Physics/top/. 
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2.3. Charge Asymmetry in Top Quark Pair Production 



2.3. Charge Asymmetry in Top Quark Pair Production 



As discussed in Section [2.2.11 top quark pairs at hadron colliders are produced via gluon-gluon 
fusion or quark-antiquark annihilation in Born approximation: 



q + q 
9 + 9 



t + t, 
t + t. 



These leading order processes obviously do not discriminate between the final state top and 
antitop, as can be seen from the respective Born partonic differential cross-sections for the two 
production mechanisms, 



da, 



qq^tt 



dcos 9 



a. 



da, 



gg-^tt 



dcos 9 



7r/3 
'•3sNc 



[l+c2+4m2) 



2s \Ncil-c^) 



3 
16 



1 + + 8mt 



32mf 



(2.34) 
(2.35) 



where 9 denotes the polar angle of the top quark with respect to the incoming parton in the 
centre-of-mass rest frame, Nc = 3, /3 = \/l — and c = {3 cos 9. Consequently, the same 
holds for the full pp — t- tt differential cross-sections. 

If, however, the processes are considered at higher order in the Standard Model, where radia- 
tive corrections from real or virtual gluon emission are introduced, a significant asymmetry can 
be generated in the differential ti cross-section, leading to a charge asymmetry in the production 
of top quark pairs [60]. This effect is caused by the interference of amplitudes which are odd 
under the exchange of the final state top quark and antitop quark. 

The dominant contribution to an overall charge asymmetry stems from interference between 
the leading order amplitude for quark-antiquark annihilation and the corresponding one-loop 
corrections (box diagrams), creating a positive contribution to the total charge asymmetry. The 
two diagrams are shown in Figure! 






Figure 2.9.: Box (left) and Born (right) diagrams contributing to the production of top quark 
pairs through quark-antiquark annihilation. 

In addition, interferences between initial state and final state gluon bremsstrahlung have to be 
taken into account, contributing negatively but typically of lower magnitude than the box-Born 
interference to the overall asymmetry. The respective diagrams are shown in Figure [2.101 






"OooooO 




Figure 2.10.: Final state (left) and initial state (right) bremsstrahlung diagrams contributing 
to the production of top quark pairs through quark-antiquark annihilation. 
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2. The Top Quark Charge Asymmetry in Context of the Standard Model 



The asymmetry arising from the combination of these contributions can be described by 
comparing the colour factor terms arising in the differential cross-sections from the two cut 
diagrams |61j shown in Figure 12.111 after averaging over initial states and summing over final 
states. The respective colour factors Ca and Cb for the two diagrams (a) and (b), respectively, 





"0W0"0"0"> 



(a) 



(b) 



Figure 2.11.: Cut diagrams contributing to the asymmetry in the production of top quark 
pairs through quark-antiquark annihilation, arising from virtual and real gluon 
emission. The respective contributions to the cross-section are odd under the 
exchange of the final state top and antitop quarks. 



can be expressed as [60] 



1 ^ A^A^A^^ A'^A^A^ 
Ca = TTrT^'^^^Tr— — — 



2 2 2 ' 2 2 2 16A^^ 



{fabc + ^abc) 



(2.36) 



and 



- ^Tryyy TYyyy - [-iabc + dabc) , 



(2.37) 

'C ^ ^ ^ ^ ^ ^ J^KJllfj 

where /^^^ = 24 and d^^^ = 40/3. Consequently, the respective contributions duA and das to 
the cross-section are odd under the exchange of the final state top and antitop quarks: 



daA {t, t) = —das {t, t) 



(2.38) 



A small contribution is introduced at the order of through interferences of different terms 
in quark-gluon scattering, 



g + q t + t + q', 



which is shown in Figure I2.12i 




Figure 2.12.: Quark-gluon scattering diagrams contributing to the hadronic production of top 
quark pairs introduced at order of a^. 
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2.3. Charge Asymmetry in Top Quark Pair Production 



The individual charge asymmetric contributions from quark-antiquark annihilation and quark- 
gluon scattering can be found in Figure [2.131 |60| as a function of the partonic centre-of-mass 
energy, quantified by the integrated forward-backward contributions 



da 



d cos 6 



^dcos9 



d(j\ 



1 d cos 6 



-dcos 9, 



(2.39) 



where i denotes the two contributions from quark-antiquark annihilation (qq) and quark-gluon 
scattering (qg). 




Vs (GeV) 

Figure 2.13.: Integrated charge asymmetric parts of the top quark pair production cross- 
section from quark-antiquark annihilation {qq) and quark-gluon scattering [qg) 
initiated processes as a function of the partonic centre-of-mass energy |60j . 



The resulting overall QCD charge asymmetry A^i can be expressed in terms of the asymmetric 
contributions cta and symmetric contributions cjs to the total production cross-section as ratio 



A 



QCD 



0"A 
0-S 



.3.^ + ataf + 



a t a, 



(0) 



(1) 
s 



+ 



.(0) 



(1) 



where 



a, 



^S,gg + ^S,qg + ^S.gg 



(2.40) 



(2.41) 



contains the respective symmetric contributions at a given order i. This asymmetry can also be 
parametrised as 



^QCD 



asA[? + aiA'r' + 



2aW 



(2.42) 



Higher order corrections arising from QCD which are not taken into account in Equation 12.401 
such as A^^\ can be evaluated in next-to-next-to-leading log approximation using soft gluon 
resummation techniques |62fl64j in order to improve the theoretical predictions on the respective 
differential cross-section contributions and the associated uncertainties. 

In addition to the asymmetry arising from the Standard Model QCD terms, further contribu- 
tions originate from mixed QCD-electroweak interference terms to the quark-antiquark annihi- 
lation process [60p65j. The tt colour-singlet configuration of the box diagram in Figure [2T9l can 
interfere with the production of top quark pairs through a photon or a Z boson (and similarly 



for the interference between initial state and final state radiation), as shown in Figure 2.14(a 



Furthermore, interferences of the gluon-7 and gluon-Z box diagrams with the leading order QCD 



amplitude as indicated in Figure 2.14(b) contribute to an additional asymmetry, which together 
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2. The Top Quark Charge Asymmetry in Context of the Standard Model 




Figure 2.1 4-: Mixed QCD-electroweak cut diagrams contributing to the asymmetry in the pro- 
duction of top quark pairs through quark- ant iquark annihilation, arising from 
interferences between the singlet-state box diagram and initial and final state ra- 
diation diagrams with electroweak tt production (left) and from the interference 
of the gluon-7 and gluon-Z box diagrams with the Born diagram (right). 



with Figure 2.14(a) leads to a total increase of asymmetry by a factor of about 1.09 [60] with 
respect to the QCD contributions. 

The overall charge asymmetry, including both the QCD and electroweak contributions, can 
be expressed as 

2~{0) , 3 (1) , 2 I 4 (2) , 

2~(0)^ 2 (0)^ 3 (1)^ 2 ' ^ ' 

a CTg -|- a%a^ + afcg -|- a%aa^ + ... 

(i) (i) 

where a)^ and o"g denote the asymmetric and symmetric contributions from QCD-electroweak 
mixing, respectively. 

Note that a similar effect emerges in Quantum Electrodynamics, where interferences at the 
order of create an asymmetry in the angular distribution of final state particles produced in 
e~^e~ collisions due to virtual radiative corrections and soft and hard photon emission |14H16| . 



2.3.1. Charge Asymmetry Beyond the Standard Model 

Numerous theoretical models predicting the manifestation of physics beyond the Standard Model 
(BSM) exist, several of which can have implications on the charge asymmetry observed in the 
production of top quark pairs at hadron colliders. The most popular models will be explained 
in the following. 



Colour-octet Gauge Bosons 

In chiral colour models |66H76| , SM colour charge is extended to contain a right-handed and left- 
handed contribution SU(3)^ x SU(3)^ to reflect the chirality. This extension implies a breaking 
of the symmetry of the diagonal SU(3)(^ group and the generation of a heavy colour-octet gauge 
boson, the axigluon. The associated coupling to quarks has a pure axial- vector structure and is 
of the same magnitude as the QCD coupling. 

Alternative models with non-chiral structure imply the existence of massive gauge bosons 
with pure vector- like couplings to quarks (colorons) [TTtiSO] or Kaluza-Klein [5T1[52] excited 
states arising from models including extra dimensions [83H91j . 

Further generalisations include the assumption of different coupling strengths for the different 
SU(3) contributions [92VI100] . leading to both vector and axial- vector couplings in the inter- 
actions of the respective colour-octet resonance and quarks. The vector and axial-vector 
coupling strengths, gy and g'^^ , respectively, lead to the following generalised term Cq/ [lOlj in 
the modified SM Lagrangian: 

Cg' = gst^Qi {g'v + QaI^) l^Glqi. (2.44) 
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2.3. Charge Asymmetry in Top Quark Pair Production 



The corresponding leading-order cross-section for top quark pair production in the annihilation 
of quark-antiquark pairs is given by |1021I103| 



dagg^u _ 2^/5/^, 2s {s- m%) 




ml) +mlTl 



o^2 .^A9Wv{C+) + 2g\g\c] (2.45) 



where C± = 1 + it Amf . The vector and axial- vector couplings of the resonances to the light 
quarks and top quarks are given by the constants gy and (7^, and by gy and 5^, respectively. 

Compared to the SM cross-section, an additional asymmetric contribution is introduced by 
terms which are odd in c. Hence, a large positive asymmetry can be generated in models where 
9a9a < or where the term Sgy^^^y^^c is dominant. A negative asymmetry on the other hand 
can be created in flavour universal models where g'^ = 9a- 

Models involving a colour-octet resonance typically require the resonance to be off-shell, either 
heavy |98j . or below the tt production threshold |1UU| . since in the intermediate mass range a 
distinct excess would emerge in the tt mass spectrum, which is not observed. Alternatively, a 
very broad resonance |lUmilU4] would be concealed due to limited statistics in the tail of the 
Mff distribution. 



Extra Weak Gauge Bosons 

Different theoretical models, such as some Grand Unified Theories (GUTs), topcolor or left- 
right |105j models predict the existence of extra weak gauge bosons, such as the W or Z' . 
Furthermore, these states can appear as Kaluza-Klein excitations in extra dimensional models 
[106yi07] . A sizable contribution to the charge asymmetry can only be introduced in W' or Z' 
t-channel interactions |105tll08ltl26j . The s-channel tt production through a Z' is suppressed 
due to the fact that the corresponding amplitudes do not create interferences with the SM 
amplitude [lUlj . 

A potentially significant asymmetry can be created by the introduction of flavour violating 
couplings into the SM Lagrangian by a term Ci^i /z' |101| , such as 

Cw'/Z' = i {g^ + l^Z'^u + i [gf + 5]f' 75) rW^d. (2.46) 

A real Z' contribution is constrained by the absence of like-sign top quark pair produc- 
tion |122j . except for very light Z' resonances. Furthermore, for a Z' leading to a sizable 
positive charge asymmetry, a large corresponding Z' coupling has to be assumed and a corre- 
sponding excess in the tail of the ti mass distribution would be generated. Such excess is not 
observed, however. Furthermore, W and Z' left-handed couplings are disfavoured by precision 
measurements in B hadron systems |127j . 



Coloured Scalars 

In addition to gauge bosons, neutral or charged coloured scalars |128j can be introduced in 
particular in the presence of larger gauge groups such as SU(5) or SO(IO) |129j . These oc- 
cur primarily in GUT models close to the unification scale. However, some of the coloured 
scalar states can manifest at lower energy scales |130lll31j . for example due to gauge coupling 
unification. 
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2. The Top Quark Charge Asymmetry in Context of the Standard Model 



Similar to extra weak gauge bosons, tt production via coloured scalars in the s-channel is not 
affected by any charge asymmetric contributions due to the absence of interferences with the 
SM amplitudes |1U1| . Consequently, only t-channel flavour- violating couplings can introduce 
a significant asymmetry in the production cross-section, such as the exchange of scalar colour 
singlets [TU51fmfTB5] . triplets [TIMfT^fnnil32lll34HTHH] . sextets [r^fn^imfTBU] and octets 
[13^1132) . Furthermore, a light scalar contribution is disfavoured due to the implied existence 
of a highly constrained new top quark decay channel {t — )• S'u) p3 tll40j . 

A generalised contribution from a coloured scalar SU(2) doublet S' to the SM Lagrangian can 
be described by an additional term £5/ |10ipi32j . given by 

Cs' = t''iigs + gp-/5)qy''u, (2.47) 

where gs and gp denote the scalar and pseudoscalar coupling constants, respectively. The 
resulting asymmetry yg/ created in the exchange is hence given by 

ys' = \/gs + 9p- (2-48) 

This contribution is typically negative for a heavy coloured scalar. However, due to the destruc- 
tive interference of the scalar contribution with the SM, a positive overall asymmetry can be 
generated. 

2.3.2. Top Quark Charge asymmetry at Hadron Colliders 

In order to quantify a potential charge asymmetry created in the Standard Model or BSM 
models, the natural choice of observable would be the production angle 6t of the top/antitop 
quarks with respect to the incoming partons from the hard scattering process, as depicted in 
Figure [2.151 The corresponding differential charge asymmetry A(cos0t) at the partonic level is 




Figure 2.15.: Top quark pair production kinematics in quark-antiquark annihilation. The ini- 
tial state partons and their momenta pq and pq and the produced top quarks and 
their momenta pt and pi, respectively, are shown. In addition, the production 
angle of the top quark Of is shown. 

given by 



in the qq rest frame, where 



^(cos Bt 



Nt{coset)-Nt{cos9t) 
Nt{coset) + Nticos Ot)' 



(2.49) 



(2.50) 
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2.3. Charge Asymmetry in Top Quark Pair Production 



and N^{cos9t) = Nt{— cos 9t) due to the symmetry of charge conjugation. Consequently, the 
integrated charge asymmetry A can be quantified such that 



However, in the strong production of top quark pairs at hadron cohiders, the production angle 
as such is not accessible experimentally due to the fact that the initial state of the partonic 
reaction is of probabilistic nature. Since the available information is limited to the hadronic 
initial state {pp or pp) and the PDFs of the protons and/or antiprotons, respectively, different 
methods to measure the charge asymmetry, making use only of the final state information of 
the hadronic collision, must be employed. 

At non-symmetric hadron colliders such as the Tevatron, where protons are brought to collision 
with antiprotons, a charge asymmetry in the production of top quark pairs as introduced in 
Equation 12.491 observed in the tt rest frame, corresponds directly to an equal-sized geometric 
forward-backward asymmetry, Ap-Q, since Nt{y) = Nj:{—y). Since this quantity is experimentally 
accessible in a direct way due to the fact that the initial directions of the proton and antiproton 
are known, the measurement of the underlying charge asymmetry in the laboratory frame is 
possible. 

At pp colliders such as the LHC, no forward-backward asymmetry is visible in the laboratory 
frame due to the intrinsic charge conjugation symmetry of the initial state collisions. However, 
since top quarks are preferentially emitted in the direction of the incoming parton, and quarks 
in the proton on average carry a larger momentum fraction than antiquarks, an excess of top 
quarks in the forward and backward regions is expected in the laboratory frame. Consequently, 
different widths of the corresponding rapidity distributions of top quarks and antitop quarks, 
and hence, the respective decay products, are predicted. The underlying charge asymmetry in 
the ti rest frame can be extracted either from the final state particles directly or by performing 
a kinematic reconstruction of the ti decay signature. 

Since an asymmetry can solely be created from the quark-antiquark annihilation contribution 
to the top quark pair production cross-section, the total charge asymmetry in both pp and pp 
collisions can be significantly diluted due to the (symmetric) gluon-gluon fusion contribution. 
The magnitude of the overall asymmetry depends strongly on the centre-of-mass energy since 
the fraction of soft gluons in the proton/antiproton PDFs and hence the probability of gluon 
interactions in the partonic reaction increases with rising hadron momentum. In addition, the 
top quark pair production cross-section in pp collisions shows a higher contribution from gluon- 
gluon fusion since interactions of sea quarks are dominant. This differs from pp collisions, 
where potential interactions of valence (anti)quarks from the colliding (anti)protons lead to an 
increased contribution from quark-antiquark annihilation in the overall cross-section. 

In order to quantify the charge asymmetry at the LHC, a proper observable, taking into 
account the potential differences in the rapidity distributions of the top and antitop quark in 
the laboratory frame has to be chosen. Different frame-invariant variables based on rapidity or 
pseudorapidity differences of the final state top and antitop quarks are typically used to measure 
the asymmetry. 

Potential observables to quantify the inclusive charge asymmetry in the tt rest frame 



based on rapidities and pseudorapidities of the final state top and antitop quarks include 



A 



Nt{cos 9t>0)- Nt{cos 9t > 0) 
iVf (cos 9t>0)+ iVf(cos 9t>0)' 
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parametrisations such as 



N+ = N{rjt-r]t>0) , N- = N{rjt - rjt < O); 

N+ = N{yt- y-t > 0) , =N{yt- yi < 0) [IZ1[2D1[I1I] ; 

N+ = Ni\7jt\-\7]t\>0) , N- =N{\7]t\-\i]t\<0)m; 

N+ = Ni\yt\-\yt\>0) , N~ = Ni\yt\ - \yt\ < 0) 



Since 



(2.53) 



Ay = 2/t - y-t = 2yf (2.54) 



and correspondingly for the pseudorapidities of the final state top and antitop quarks, it follows 
that 

_ NjAy > 0) - NjAy < 0) 

" NiAy>0)+N{Ay<0)- ^^'^^^ 

At the Tevatron, the predicted charge asymmetry within the Standard Model for the given 
observable Aij-, evaluated for a centre-of-mass energy of 1.96 TeV, is 

Aa = 0.087 ± 0.010 [Il2]. (2.56) 

In the following analysis, the charge asymmetry will be quantified using a parametrisation 
based on the difference of absolute rapidities of the top and antitop quarks, defined as 

, ^ N{\yt\-\yt\>0)-N{\yt\-\yi\<0) 

N{\yt\-\yt\>0) + N{\yt\-\yt\<Oy ^ ' > 

At the LHC, the predicted charge asymmetry within the Standard Model for the given ob- 
servable Ac, evaluated for a centre-of-mass energy of 7 TeV, is 

Ac = 0.0115 ± 0.0006 [Hg. (2.58) 

A summary of the predicted charge asymmetries for Tevatron and LHC measurements for 
various BSM models can be found in Figure 12.161 Potential regions in the phase space of inclusive 
charge asymmetry from new physics, indicated by the variable A^"" at the LHC plotted against 
the associated forward-backward asymmetry Ap*^ at the Tevatron, are highlightecfl. Different 
generalised predictions and parametrisations |116lll27j of Z' and W models, scalar triplet (uj^) 
and sextet {0.^) models, a generalised colour-octet resonance model {G/i) and a colour-singlet 
Higgs-like isodoublet (j) are shown. The same model predictions are shown for a high tt invariant 
mass region. 

Depending on the observed asymmetries at the Tevatron and the LHC, the exclusion of dif- 
ferent theories can be possible. Recent measurements of the charge asymmetry by the CDF and 
D0 collaborations at the Tevatron indicate large partonic asymmetries Aj-t in the tt rest frame 
of 

Ait = 0.201 ± 0.065 (stat.) ± 0.018 (syst.) [IB] (2.59) 

in the combination of the dileptonic and semileptonic decay channel, indicating a 2.9a excess, 
and 

Att = 0.196 ± 0.065 (stat.-Fsyst.) [20], (2.60) 



*Note that ^pB* and A^"" denote only the respective contributions to the overall charge asymmetry and forward- 
backward asymmetry originating from the corresponding BSM model. The predicted Standard Model contribu- 
tion in the respective variables is subtracted (and hence corresponds to ^pB^ = and = 0, respectively). 
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Figure 2. 1 6. : Predicted charge asyrametries at the Tevatron and LHC for various BSM models 
|116[|127j . The inclusive charge asymmetry originating from new physics, A'^'" , 
at the LHC vs. the forward-backward asymmetry at the Tevatron (left) 

is shown. Furthermore, the identical predictions in a high invariant mass region 
where Mti > 450 GeV (right) for the different models in the created phase 
space is shown. The Standard Model prediction corresponds to Ap™ = and 
= 0, respectively. 



respectively, indicating a 1.9a" excess above the Standard Model prediction. Furthermore, larger 
deviations for high tt invariant masses |141j and for high rapidity differences |143| have been 
observed. In particular, an asymmetry of 

Ait = 0.475 ± 0.114 (stat.+syst.) [HI] (2.61) 

for tt invariant masses above 450 GeV has been observed, indicating a 3.4o" deviation from the 
Standard Model prediction, as shown in Figure [2.171 



, tt parton -level 
CDF data 5.3 fb"' 

0.4-1- — tf NLOQCD 



0.2- 



0.0-- 



-0.2- 



450 GeV/c' 



Figure 2. 1 7. : Summary of results obtained in a tt invariant mass dependent measurement of 
the charge asymmetry at CDF [141] . A deviation of 3.4cr is observed in the high 
Mti region. 



Similar measurements have been performed by the ATLAS and CMS collaborations at the 



23 



2. The Top Quark Charge Asymmetry in Context of the Standard Model 



LHC. In previous ATLAS measurements, a charge asymmetry in the variable Aq of 

Ac = -0.024 ± 0.016 (stat.) ± 0.023 (syst.) [22] (2.62) 

has been measured, while a measurement at CMS using the corresponding observable based on 
pseudorapidities instead of rapidities indicates a charge asymmetry of 

Ac = -0.016 ± 0.030 (stat.)+°;°2i (syst.) [H], (2.63) 

in compatibility with the Standard Model prediction. 
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This chapter describes the technical details of the Large Hadron Collider and the ATLAS exper- 
iment, focusing on the detector subsystems and their properties after a short general overview. 
In addition, the ATLAS trigger system is explained in more detail concerning technical and 
functional parameters. 

3.1. LHC and ATLAS Technical Overview 

The Large Hadron Collider (LHC) |144j is situated at CERN, the European Centre for Nuclear 
Research near Geneva, Switzerland, and is the technologically most advanced particle acceler- 
ator so far. It is designed as a proton-proton accelerator with the potential to accelerate and 
collide heavy ions as well. Its construction started in 1999 and physics operation commenced in 
November 2009 with the first proton-proton collisions. The collider has been constructed in the 
former accelerator ring of the Large Electron Positron Collider (LEP), which is about 27 kilo- 
metres in circumference and 100 to 125 metres below ground level, partly making use of the 
already existing infrastructure of the LEP ring as well as two existing caverns for experiments. 
A schematic view of the accelerator ring is shown in Figure 13.11 




Figure 3.1.: Layout of the LHC accelerator complex, including four of the LHC experiment^. 

The accelerator complex incorporates six experiments, two of them being multi-purpose ex- 
periments, ATLAS |145j and CMS |146| . while the other four are designed for more specific 
fields of research. Having two independently designed multi-purpose detectors is vital for cross- 
confirmation of any potential discoveries made and naturally allows to combine the results of both 
experiments. Among the four other detectors, the LHCb |147j experiment focuses on b physics, 
while ALICE [HH] has been designed for heavy ion physics. TOTEM [HH] and LHCf [l50] are 
both designed to conduct studies of forward physics, including soft and hard diffractive processes 
and low-x QCD. 

^ Image taken from wikipedia (public domain). 
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Protons are extracted from hydrogen atoms in a duoplasmatron |151j and then accelerated by 
a hnear accelerator (up to 50 MeV) before being injected into the first circular accelerator, the 
Proton Synchrotron Booster (50 MeV — )• 1.4 GeV). Afterwards, they are injected into the Proton 
Synchrotron (PS, 1.4 GeV — )• 26 GeV) and subsequently into the Super Proton Synchrotron 
(SPS, 26 GeV — t- 450 GeV), before being transferred into the main LHC ring. The two LHC 
proton beams deliver an energy of 3.5 TeV each (7TeV at design specifications) and bunches of 
about 10^^ protons can be brought to collision at a bunch crossing rate (BCR) of 40 MHz within 
one of the various detectors. 

Collision interactions are typically characterised by the instantaneous luminosity C, which 
relates the cross-section o" of a given process to the corresponding event rate A^, given the 
experimental acceptance A and measurement efficiency e: 



C 



N 
aAe' 



(3.1) 



At the design luminosity of £ = 10^^ cm~^s~^ and the given collision parameters, this leads to 
a total of about 23 proton - proton collisions per bunch crossing on average, which implies an 
overall interaction rate in the GHz regime. 

The ATLAS detector is designed to investigate a wide range of physics processes, including 
the measurement of well-known Standard Model processes, the search for the Higgs boson, extra 
dimensions, and particles that could constitute dark matter. 

It is 44 m in length, 25 m in height and 25 m in width, with a total weight of about 7 000 tonnes. 
As a comparison, CMS, the second LHC multi-purpose experiment, weighs about 12 500 tonnes 
while being 21 m long, 15 m wide, and 15 m high. ATLAS features an onion-like structure, which 
will be discussed in detail in Section 13.31 A summary of the ATLAS and the LHC specifications 
compared to other accelerator/detector combinations is shown in Table [3Tl 





Type 




BCR 


ATcH 


Event Size 


Year 


ATLAS, LHC (CERN) 


PP 


7 TeV 


40 MHz 


~ 108 


- ISOOkB 


2009 


CDF/D0, Tevatronf.(FNAL) 


PP 


1.96 TeV 


2.5 MHz 


- 10^ 


- 250 kB 


2001 


ZEUS, HERA (DESY) 




0.318 TeV 


10 MHz 


- 10^ 


- 100 kB 


1992 



Table 3.1.: Comparison of LHC/ ATLAS to other accelerators/detectors. The particle types 
brought to collision, the centre-of-mass energy -/s, the bunch crossing rate (BCR), 
the amount of readout channels (A^ch), the average amount of data per event and 
the year of startup of the accelerator are shown. 



As can be seen from the comparison table, the LHC exceeds the bunch crossing rate of 
the largest particle accelerator up to 2009 (Tevatron in Run II, Fermi National Accelerator 
Laboratory) by a factor of seven. Furthermore, the amount of detector readout channels used 
at ATLAS increased by two orders of magnitude with respect to CDF/D0, together with the 
average amount of data per event increasing by about one order of magnitude. This results in a 
much higher total raw data throughput rate to be processed by the readout and data acquisition 
system. 

As a consequence, a sophisticated trigger system designed for the reduction of data rates is 
needed to facilitate analyses. This can only be achieved by selecting a subset of potentially 

■^Values refer to Tevatron Run II. 
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relevant events out of the large amount of collisions that take place at LHC/ATLAS, rejecting 
a large fraction of raw data. 



3.2. ATLAS Coordinate System 

The ATLAS coordinate system is a right-handed coordinate frame with the x-axis pointing 
towards the centre of the LHC ring and the z-axis being directed along the beam pipe, while 
the y-axis points upwards (slightly tilted with respect to the vertical direction (0.704°) due to 
the general tilt of the LHC tunnel). In this context, the pseudorapidity can be introduced as 

rj = — In tan - (3.2) 

with 6 being the polar angle with respect to the positive y-axis. For massive objects such as 
jets, the rapidity is used, given by 

Li addition, the transverse momentum of a particle in the detector is defined as the 
momentum perpendicular to the z-axis: 



PT = y/pl+pl. (3.4) 
Furthermore, the azimuthal angle (p is defined around the beam axis. 



3.3. The ATLAS Detector Subsystems 

In order to allow for reliable detection of particles and measurement of their properties, the 
ATLAS detector requirements include 

• a good hermiticity with respect to detector acceptance, 

• a high spatial and timing resolution, in particular to minimise occupancy of individual 
detector components, to measure pT with high resolution and to allow for distinction 
between different particles, 

• a low material budget to minimise particle interaction and energy loss with non-active 
detector materials. 

ATLAS has a cylindrical shape with layers of detector components arranged in axial succession. 
Each of these layers is designed to detect different types of particles which are mostly originating 
from the primary interaction point of the proton beams at the centre of ATLAS. As they 
travel throughout the detector, they can be measured by its successive layers. The different 
detector subsystems are shown in Figure [3^2] and constitute, from the innermost to the outermost 
layer, the inner detector, the solenoid magnet, the electromagnetic calorimeter, the hadronic 
calorimeter, the toroidal magnet, and the muon spectrometer. 

The detectors are complementary: charged particles are detected in the innermost layers by 
their hits in the tracking chambers, where the particle trajectory is bent by the magnetic field 
of the superconducting solenoid magnet. Using this tracking information, the momentum of the 
particles can be determined. Around the magnet, the electromagnetic and hadronic calorimeters 
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Muon Detectors Tile Calorimeter Liquid Argon Calorimeter 




Toroid Magnets Solenoid Magnet SOT Tracker Pixel Detector TRT Tracker 



Figure 3.2.: The ATLAS detector subsystemtH. 

are designed to measure the energy of particles. These are brought to a stop in the calorimeter by 
interaction with the detector material [23] (ionisation), thus depositing all of their energy, which 
is measured in the calorimeter cells. Finally, the muon chambers allow for additional momentum 
measurements of muons, which penetrate all other layers of the detector only depositing very 
little energy in the detector material. The measurements in the muon chambers are performed 
using the shape of the tracks, which are bent by the magnetic field of the toroidal magnets. 

3.3.1. Inner Detector 

In the following, the inner detector |145^I152] is described in more detail. It is situated near 
the interaction point to allow for high precision measurement of charged particle trajectories. 
It covers a pseudorapidity range of |r/| < 2.5 and consists of three subsystems, the silicon pixel 
detector, the semiconductor tracker (SCT) and the transition radiation tracker (TRT). The 
innermost layers of the inner detector (three in the cylindrical barrel region, three endcap disks 
on each side of the forward region) constitute the pixel detector, which is designed to measure 
particle vertices and extract track momenta from the reconstructed particle hits in the detector 
layers. Due to its close proximity to the primary interaction point, a very high spatial resolution 
of the pixel detector is required, which is achieved by very small pixel sizes of 50 \mi x 400 |J.m 
and 50 |J.m x 600 |a.m (about 8 • 10'' readout channels in total), with the pixel detector covering 
a total area of 2.3 m^. 

Around the pixel detector, the semiconductor tracker (or silicon strip tracker) measures the 
momentum of charged particles. It consists of four barrel layers and nine endcap wheels on each 
side, covering a total area of 61.1 m^ and making use of about 6.3 • 10^ readout channels. 

The outermost part of the inner detector is constituted by the TRT, which consists of straw 

^Image taken from http://www.atlas.ch/. 
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Figure 3.3.: Cross-sectional view of a quarter-section of the ATLAS inner detector showing 
each of the major detector elements alongside with its active dimensions and 
envelopes [145] . 



tubes with a diameter of 4 mm and a maximum length of 80 cm, filled with an ionisable gas. 
The barrel tubes are divided in two at the centre and read out at each end to reduce occupancy. 
The ionisation charges created by charged particles travelling through the gas filled tubes are 
used to enhance track pattern recognition and to improve momentum resolution of the objects 
identified in the pixel detector and semiconductor tracker with an additional average of 36 hits 
per track. Furthermore, its function is to distinguish electrons and pions making use of the 
different amount of transition radiation emitted by these particles when crossing the boundary 
surface of two media with different dielectric constants (in this case a special radiator foam with 
a large amount of air bubbles to achieve a maximum material transition surface) . The transition 
radiation tracker has a total of 351 000 readout channels. 

The resulting tracking performance of the inner detector subsystems for single particles and 
particles in jets can be found in Table [3^21 

3.3.2. Calorimeters 

Around the solenoid magnet (which is described in more detail in Chapter 13.3. 4p . an electromag- 
netic (EM) liquid argon sampling calorimeter is used to detect and identify electromagnetically 
interacting particles and to measure their energy |1451I153T - I155| . It also allows, in combination 
with the hadronic calorimeters, for reconstruction of hadronic jets and the measurement of the 
missing energy of an event. 

The EM calorimeter covers a pseudorapidity range of \'q\ < 3.2 and comprises several layers 
of accordion-shaped kapton-copper electrodes and lead absorber plates shrouded in stainless 
steel, with the gaps in between filled with liquid argon at a temperature of 87 K. Whenever 
an electromagnetically interacting particle passes through one of the lead absorber plates, it 
creates a particle shower that ionises the liquid argon. Exposed to an electric field, the drifting 
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Track Parameter 


0.25 < \t] 

(Tx(oo) 


1 < 0.50 

px[GeV] 


1.50 < \t] 
crx(oo) 


1 < 1.75 

px[GeV] 


jVIoniGiituni (l/j^x) 


0.34 TeV"^ 


44.0 


0.41 TeV^^ 


80.0 


Aziniuthal angle {(j)) 


70 ixrad 


39.0 


92 ^irad 


49.0 


Polar angle (cot 6) 


0.7 X 10-3 


5.0 


1.2 X 10-3 


10.0 


Transv. impact parameter (do) 


10.0 [im 


14.0 


12.0 nm 


20.0 


Longit. impact parameter (zq x sm0) 


91.0 M-m 


2.3 


71.0 i^m 


3.7 



Table 3.2.: Track parameter resolutions at infinite momentum (Tx(oo) and the transverse mo- 
mentum px for which the intrinsic and multiple-scattering contribution equals the 
intrinsic resolution. The momentum and angles correspond to muons, while the im- 
pact parameters correspond to pions. The values are shown for two rj regions, one 
in the barrel inner detector (where the amount of material is close to its minimum) 
and one in the endcap (where the amount of material is close to its maximum) jl45| . 



Tile Barrel Tile Extended Barrel 




Figure 3. 4.: Overall layout of the ATLAS calorimeters |145j . 



ionisation charges induce a signal in the electrodes due to capacitive coupling. The resulting 
signal is sampled and digitised at 40 MHz, corresponding to the bunch crossing rate. The EM 
calorimeter has a total of 170 000 readout channels, and provides an energy resolution of 



<M _ (lO-l ± 0-4) % 

E [GeV] 



(stochastic) 0.2 ±0.1% (constant) [US], 



as measured in particle test beams. The total module thickness in the barrel region corresponds 
to at least 22 radiation lengths {Xq), increasing from 22 Xq to 30 between \ri\ = and 
I??! = 0.8 and from 2AXq to 33 between \r]\ = 0.8 and \r]\ = 1.3. In the endcaps, the total 
thickness is greater than 2AXq except for |?7| < 1.475, increasing from 24Xo to 38 in the 
outer wheel (1.475 < \r]\ < 2.5) and from 26 Xq to 36 in the inner wheel (2.5 < \r]\ < 3.2). 
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In analogy, the hadronic calorimeter (hCAL) |145lll53lll54j is designed to measure the energy 
of hadronic particles that can penetrate the EM calorimeter. This sampling calorimeter consists 
of iron absorbers for showering which are interleaved with plastic scintillator tiles (thus being 
referred to as tile calorimeter) in the barrel part of the detector (|?7| < 1.7). The scintillator 
tiles emit a shower of photons whenever charged particles pass through them due to excitation 
of the atoms in the scintillating material and subsequent emission of visible or UV photons. 
These light pulses are carried by optical fibres to photomultiplier tubes and converted to electric 
signals. The total number of tile calorimeter readout channels is of the order of 10 000. 

Since the scintillating tiles are very sensitive to radiation damage, liquid argon is used as 
sampling medium together with copper absorbers in the forward endcap regions (1.5 < \ri\ < 3.2) 
in close proximity to the proton beams. This provides improved radiation hardness in the region 
of increased particle flux. The total number of channels for both endcaps is 5 632. For the same 
reason, a high density copper/tungsten absorber liquid argon forward calorimeter (FCAL) covers 
the pseudorapidity region 3.1 < |?7| < 4.9, with an additional 3 524 channels for both forward 
regions together. 

In order to estimate the performance of the hadronic calorimeter, test beam studies were 
conducted, showing an energy resolution of 



HAD 

^ IT 



(56.4 ±0.4)% , , , ^ , , ,— ^ 

— ^ (stochastic) © 5.5 ±0.1% (constant) [US] 



K [GeV] 

for pions, and a radial depth of approximately 7.4 interactions lengths (A) for the tile calorimeter. 
The hadronic endcaps show an energy resolution of 



^HEC ^ (21.4 ±0.1)' 



'0 



(stochastic) |145| 



Ee [GeV] 
for electrons (the constant term being compatible with zero), and 



g^^^ _ (70.6 ±1.5)% 
E^ ~ ^E [GeV] 



(stochastic) © 5.8 ± 0.2 % (constant) [T35] 



for pions. 

The jet energy resolution for the overall calorimeter system is described by the parametrisation 

For central jets in the region 0.2 < \r\\ < 0.4, it is a ~ 60 % VGeV (stochastic), c « 3 % (constant) 
and the noise term h increases from 0.5 GeV to 1.5 GeV from barrel to endcap ranges |145j . 




3.3.3. Muon Chambers 

The muon chambers |145pi54lll56j are designed to detect muons which are able to pass all other 
detector systems depositing only a small amount of energy in the material due to the fact that 
muons in the GeV regime are approximately minimum ionising particles. 

Being deflected by the magnetic field in the detector, it is possible to determine the muon 
momentum and sign of electric charge by measuring its trajectory as it passes through the 
tracking chambers. The muon spectrometer is also designed to trigger on these particles, utilising 
dedicated trigger chambers. The driving performance goal is a standalone transverse momentum 
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resolution of approximately 10 % for 1 TeV tracks, which translates into a sagitta along the z 
axis of about x = 500 \xm, to be measured with a resolution of < 50 [xm. The sagitta x is 
given by 

e , , 

X = ^--Rcos-, (3.5) 

where R is the radius of the track curvature and Q is the angle enclosed by the outermost of 
three equidistant points along the track, as can be seen in Figure [331 

To achieve high spatial tracking resolution, three layers (stations) of drift chambers (precision 
chambers) are employed both in the barrel and in the endcap region. In the barrel these chambers 
are arranged in concentric cylinders, with the radii of the detector layers being at about 5 m, 
7.5 m and 10m, covering a pseudorapidity range of \r]\ < 1.0. Two of these layers are placed 
near the inner and outer field boundary, while the third is situated within the field volume. The 
muon momentum is determined from the track sagitta. In this region, exclusively monitored 
drift tube chambers (MDTs) are used. 

Due to the magnet cryostats in the endcap region, however, placing one station within the field 
is not possible. Hence it is necessary to rely on a point-angle measurement to determine the track 
momentum in this detector region (a point in the inner station and an angle in the combined 
middle-outer stations). The endcap layers are arranged in four concentric discs at 7m, 10 m, 
14m, and 21-23 m from the detector origin, covering a pseudorapidity range of 1.0 < \rj\ < 2.7. 
Here, in addition to MDT chambers, cathode strip chambers (CSCs) are used in the innermost 
layer of the inner station due to radiation hardness requirements. 



magnetic field area 




Figure 3. 5. : Illustration of the track sagitta x and its geometric relation to the bending radius 
R for a track of a particle with a given momentum p travelling through a homo- 
geneous magnetic field, shown for the case of three equidistant track hits along 
the curvature. 

The only exception to the continuous rj coverage of the muon chambers is made at |r/| < 0.05 in 
the R-(j) plane to allow for cable and service outlets for the inner detector, the central solenoid 
and the calorimeters (central gap). Further regions of reduced acceptance can be found at the 
feet of the detector. 

Regardless of the high spatial resolutions, the timing resolution of the precision chambers as 
shown in Table 13.31 is too low due to the drift time to ensure differentiation between muons from 
subsequent bunch crossings for the trigger (where the scale of the required timing resolution is set 
by the bunch crossing interval of 25 ns). Thus it is necessary to employ additional drift chambers 
with high timing resolution (at cost of spatial resolution), the trigger chambers. These provide 
a fast momentum estimate and are primarily used for the trigger (since their spatial resolution 
is too low with respect to the precision chambers), which will be described in Chapter 13.41 

In the barrel region, two out of three resistive plate chambers (RPCs) are placed directly in 
front of and behind the central MDT, while the third is situated directly below or above (ac- 
cording to the mechanical constraints in the respective region) the outermost precision chamber. 
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The RPCs are furthermore used to determine the second coordinate for the MDT chambers (in 
tube wire direction). In the endcaps, three layers of thin gap chambers (TGCs) are located near 
the central endcap MDT layer for triggering. 

A comparison of the different chamber technologies, necessary to allow for both fast triggering 
and high precision measurements, is shown in Table 13.31 





MDT 


CSC 


RFC 


TGC 


z/R resolution [[im] 


35 (z) 


40 (R) 


10 X 10^ (z) 


(2 - 6) X 10^ (i?) 


# readout channels 


354000 


30 700 


373 000 


318 000 


# chambers 


1150 


32 


606 


3 588 


(p resolution [mm] 




5 


10 


3-7 


timing resolution [ns] 




7 


1.5 


4 



Table 3.3.: Comparison of muon chamber technologies. Both the trigger chambers (RPC and 
TGC) and precision chambers (MDT and CSC) and their respective resolution and 
timing performance are shown. 



3.3.4. Magnet System 

Outside the inner detector follows a solenoid magnet [14511157] . which is used to bend the 
trajectories of charged particles on their way through the inner detector, making it possible to 
measure the particle momentum with high resolution. Its axial field strength is about 2 T (peak 
2.6 T) using low-temperature superconducting cables cooled down to 1.8 K with liquid helium 
during operation with a nominal current flow of 8 000 A. 

In addition, a large toroidal magnet |145tll58j consisting of eight superconducting coils in 
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Figure 3.7.: Schematic view of the ATLAS solcnoidal (inner cylinder) and toroidal magnets 
(outer coils) |156j . 

the barrel region and eight more at each of the forward regions extends throughout the muon 
chambers, providing a magnetic field strength of 4T (peak 4.7 T). The whole toroid system 
contains over 70 km of superconducting cable, allowing for a design current of 20 000 A with a 
stored energy of above 1 GJ. To minimise multiple scattering of the muons, the toroid design 
incorporates an air core. 

Similar to the solenoid magnet, its purpose is to bend the trajectories of muons in order to 
measure their transverse momentum (in combination with the tracking information from the 
inner detector). 

3.4. The ATLAS Trigger System 

The design luminosity of 10^^ cm~^s~^, in combination with the bunch crossing rate of 40 MHz 
and the amount of protons contained in each single bunch, leads to a proton-proton collision 
rate in the GHz regime. This corresponds to an extremely high theoretical raw data rate of 
about 1.5PBs~^. Being able to store only a fraction of this amount of data on storage media 
(~ (300 — 500) MBs~^) and only a small fraction of these collisions being useful for analysis, the 
ATLAS trigger system |145yi54tll59| has been designed to reduce the initial data rate by several 
orders of magnitude. Figure \2M shows the total rates of several physics processes in comparison 
to the total interaction rate. As an example, the frequency of the Standard Model ti production 
at y/s = 7TeV is of the order of 1 mHz, constituting only a small fraction of the total amount 
of the raw data rate, making the selection of this and other physics processes a crucial task. 

In order to achieve such a reduction and to select only relevant physics events / processes, 
ATLAS uses a three-level trigger system for real-time event selection (with the last two levels 
being referred to as high-level trigger), while each trigger level refines the decisions of its prede- 
cessor. An overview of the different trigger levels and the global structure is shown in Figure 13.81 
and will be discussed in more detail in the following sections. 

3.4.1. Level 1 Trigger (LVLl) 

The LVLl trigger is completely hardware-based, where highly specialised components are de- 
ployed, including field programmable gate arrays (FPGAs), application specific integrated cir- 
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Figure 3.8.: Block diagram of the trigger/DAQ system. On the left side the typical collision 
and the data equivalent at the different stages of triggering are shown, while 
in the middle section the different components of the trigger system are shown 
schematically |154) . The right side of the graphic gives a short summary of the 
operations and the technologies used at the respective level. 



cuits (ASICs) and reduced instruction set computing (RISC) chips. Since most of this hardware 
is integrated directly into the particular detector components in order to reduce material occur- 
rence from cabling and additional readout electronics, the LVLl trigger system is required to be 
highly parallelised. 

Being based on the muon and calorimeter system only, the LVLl trigger |16U| performs an 
initial selection on the basis of the hits in the muon trigger chambers and calorimeters. In the 
muon chambers, low-pT and high-pT muons are identified by measuring the tracks in the trigger 
chambers (RPCs and TGCs) using coincidence windows for discrimination. Low-px muons have 
a smaller bending radius, thus allowing for detection by two layers in close proximity to each 
other (Moreover, track-hit matching can be difficult if the in-plane hit distance in the coincidence 
layers is too large due to the curvature of the track). In the barrel region, low-pT muons are 
identified by the consecutive layers RPCl and RPC2 (The number after the chamber type 
identifies the layer). In contrast, high-px muons produce an almost straight track and therefore 
the coincidence layers used should be as separated as possible to allow for the measurement of the 
track radius. Consequently, high-px muons are measured by the combination of the outermost 
RPCl and RPCS hits in the barrel. In the endcaps, low-px muons are identified by the TGC2 
- TGCS coincidence window, while TGCl and TGCS are used to identify high-px tracks. This 
is shown in Figure IS.9| where a quadrant layout of the muon trigger chambers is shown. 

Classification of transverse momentum is achieved by using large lookup tables of track hits 
to find an estimate for the track momentum. Different exclusive low and high px thresholds are 
defined (e.g. 10, 11, 20 and 40GeV) and can be modified if necessary. 

The calorimeter trigger selection is based on low resolution information from all ATLAS 
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Figure 3. 9. : Layout of the muon trigger chambers |156j . A quarter cross-section in the bending 
plane with typical low and high muon tracks and the corresponding coinci- 
dence windows in the different layers for barrel (RPCl-3) and endcap (TGCl-3) 
are shown. 



calorimeters and is designed to identify high electrons and photons, hadron jets and the 
total transverse energy alongside with large missing transverse energy fJ^, where this is defined 
by the sum of all vectored energy depositions Ej- in the transversal plane: 

^T = -^^T- (3.6) 

This definition arises from the fact that the initial transverse momentum of the incoming protons 
is approximately zero and due to conservation of energy in the transverse plane, the total vector 
sum of final state transverse energies has to be zero as well. Furthermore, isolation requirements 
based on calorimeter information are available for the calorimeter trigger on this trigger level. 

Since the LVLl trigger operates synchronously with data taking, the latency for a decision is 
about 2.5 [Ls, which is achieved by making use of pipeline memories despite the bunch crossing 
and data taking frequency of 40 MHz. The trigger decision is derived while the information 
from the sensors is kept in the buffer memory for about 100 bunch crossings. At the end of the 
latency time, the readout data is rejected or accepted (after leaving the pipeline memory). 

Only in the latter case the geometric information of the triggered object is forwarded to the 
next trigger level as region of interest (ROI), which includes spatial position (77 and (j)) and 
Et estimates of the embedded objects (e, /u, r, 7 and jet candidates) as well as global energy 
information (Et and ^t) for further analysis. This way, the data rate is reduced to about 
(75-100) kHz at LVLl. 

3.4.2. Level 2 Trigger (LVL2) 

After a LVLl accept, events are read out from the pipeline memories and stored in readout 
buffers (ROBs) until being processed by the LVL2 trigger jl61| . This trigger level uses the 
ROIs from the previous level to further reduce the data rate to about 1 kHz. This is achieved 
by analysing each ROI in the detector system from which it originated, accessing the data 
from other detector subsystems in the ROI in addition, including the inner detector tracking 
information, the full-granularity calorimeter hits and the precision chamber measurements from 
the muon spectrometer. 
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With an increased processing latency of 10 ms and the data rate aheady reduced at LVLl, 
the LVL2 trigger allows for more complex algorithms being applied to the trigger objects and 
the information from the respective ROIs. As a part of the high-level trigger, the LVL2 stage 
is completely software-based and runs on dedicated computing farms. The trigger decision is 
performed by an event-driven sequential selection procedure using the detector data. Despite 
the LVL2 trigger providing more time to conduct the trigger decision than the LVLl stage, the 
selection algorithms still have to be kept simple and efficient. Hence, the sequence of the different 
algorithms / requirements is determined by their complexity, simple ones (with respect to CPU 
time and memory) are executed first, while each algorithm uses the result of its predecessor 
{seeded reconstruction). 

3.4.3. Event Filter (EF) 

In the final execution level of the trigger, which is software based and runs on CPU farms, the 
global event is collected from the RGBs. The EF accesses the complete event information and 
all detector subsystems using full granularity. An event reconstruction at trigger level is possible 
by using similar algorithms as it is the case for offiine reconstruction, accessing calibration and 
alignment information from databases. 

In addition to the reconstruction of the event as a whole, the EF performs extended tasks 
that are not possible at earlier trigger levels, such as vertex reconstruction, final track fitting 
and algorithms requiring larger ROIs than available at LVL2 (e.g. the calculation of global ^t). 

The latency of the EF is of the order of seconds, after which events passing this final trigger 
level are stored permanently for further analysis or, in some cases, may be redirected to special 
storage elements as well (if needed for calibration or alignment exclusively) . 

3.4.4. Trigger Implementation 

The information used to conduct the LVLl decision is given in terms of multiplicities of trigger 
items resembling candidates for physical objects (like electrons, muons, jets,. . . ) detected in 
the calorimeters or muon trigger chambers which have sufficiently high px- These are sent to 
the central trigger processor (CTP) together with threshold information on global energy sums. 
The delivered multiplicities are discriminated against corresponding requirements or conditions, 
leading to logical values for each condition within so-called trigger menus. Examples for LVLl 
trigger items are 

EM25i: electromagnetic, pT > 25 GeV, isolated, 
MU20: muon, pT > 20 GeV. 

The final LVLl event decision is derived from the values of the defined trigger items by applying 
a logical OR. If the event is accepted, the ROIs for all trigger items and the contained trigger 
elements above the defined thresholds are delivered to the high-level trigger, where they are used 
as seeds for further processing. 

The step-by-step execution of trigger algorithms as it is performed on LVL2 and EF level is 
called a trigger chain, consisting of different intermediate trigger signatures (cf. trigger items on 
LVLl), where the successive trigger algorithms keep refining these signatures in the course of a 
trigger chain. Examples for such trigger signatures, as they have been used in this analysis, are 

EF_e20 -medium: electron, E't > 20 GeV, medium object definition, 
EF_mul8: muon, p'l > 18 GeV. 
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The individual decisions can also be logically combined to more complex trigger items. 

As an additional requirement for the high-level trigger, parallel execution of different trigger 
chains without interference is required to be possible in independent slices (e.g. electron slice 
and muon slice) to ensure transparency and scalability of the trigger system. 

The information of individual high-level trigger signatures is organised in trigger streams, 
grouping similar triggers based on the respective purpose, priority, and procedure of processing 
of the respective events. An inclusive model is chosen to allow for trigger items to be contained 
in multiple trigger streams in parallel. Raw data or physics streams contain the data events 
selected for full reconstruction and later analysis and correspond to the respective types of 
trigger items, such as the physics_Egamma and physics_Muon streams which contain events 
triggered by electron/photon or muon triggers, respectively. Further trigger streams exist for 
calibration and monitoring purposes, such as the Express stream, which contains a small subset 
of relevant events for very fast reconstruction in order to allow for real-time monitoring of the 
data taking and trigger system. 

3.5. Underlying Event and Pile-Up 

At hadron colliders such as the LHC, a multitude of partons is involved in the hadronic collisions 
due to the substructure of the colliding particles, while typically only one parton from each of 
the incoming hadrons is involved in the hard scattering process and has sufficient energy to 
create high energy/momentum final state particles. Nevertheless, the remaining partons still 
contribute to the final event signature through low momentum transfer interactions. These 
additional reactions which occur in parallel to the hard scattering event are denoted underlying 
event. 

Furthermore, additional contributions to the final event signature arise from the interactions 
of other protons within the colliding bunches. Despite the fact that the probability of multiple 
hard scatterings occurring within one bunch crossing is relatively low, the likelihood of soft 
interactions between the constituent partons originating from additional proton-proton collisions 
in the same bunch crossing increases significantly with instantaneous luminosity C. In addition 
to this in-time pile-up contribution, final state particles from different bunch crossings can lead 
to additional signatures if the identification of the correct bunch crossing is unsuccessful {out- 
of-time pile-up). 

Both the underlying event and the pile-up contributions play a significant role in the final 
state of the observed collisions and hence have to be included in the simulation of signal and 
background processes to achieve a proper modelling of the data taken. This modelling is in 
particular difficult due to the fact that both the underlying event and the pile-up cannot fully be 
described in perturbative QCD, since the particles created in such processes typically have very 
low momentum/energy. Furthermore, due to time-dependent changes in the LHC environment 
parameters, the simulation has to be adapted continuously to match the respective setup for a 
given subset of data taken. 
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This chapter covers the determination and different approaches of applying trigger and recon- 
struction efficiencies in analyses. In addition, the object definitions used in this analysis are 
described. 

4.1. Data Quality 

In order to ensure that only data taken under well defined and stable conditions is taken into 
account for physics analyses, dedicated online |162j and offline [163] monitoring systems ensure 
data integrity and quality. 

The online data quality monitoring accesses real-time detector status information and makes 
use of events from the EXPRESS trigger stream to provide several low-level quantities and dis- 
tributions. This allows for a quick response to problems with the LHC beam conditions or the 
detector that may arise during operation. 

The data quality offline monitoring uses a first reconstruction performed in order to identify 
and record problems in the detector hardware and the data acquisition and processing. All 
relevant information from the individual detector systems and reconstructed event quantities 
are combined into a small set of key numbers and distributions to allow for both automatic and 
manual monitoring. 

Information from the online and offline data quality monitoring as well as feedback from the 
individual shift crews is combined into a database containing LHC beam conditions, detector 
status and data flow information which can be used to create lists of runs usable for analyses 
(GoodRunsLists) , containing a set of data taking run and luminosity block (LB) information. 

4.2. Trigger and Reconstruction Efficiencies 

Trigger and detector acceptance and response are represented by the respective trigger and 
reconstruction efflciencies, which are typically both estimated in Monte Carlo simulations and 
measured with data driven methods. However, both the reconstruction and trigger simulation 
used for Monte Carlo samples typically do not fully reflect the actual conditions due to limitations 
of the simulation modelling. Hence, discrepancies have to be corrected for in the Monte Carlo. 
This ensures that the actual detector performance is reflected in the corrected Monte Carlo 
distributions of relevant physical quantities. 

Since for Monte Carlo samples the detector simulation/reconstruction does not retain the 
information which links objects from the Monte Carlo generator level to reconstructed objects 
or trigger objects, it can be necessary to perform a matching procedure to identify corresponding 
objects and their afflliation at the different trigger levels and at offline reconstruction level. This 
is achieved by a geometric matching in AR, deflned for example for the matching of objects at 
reconstruction level to the corresponding generated (true) object as 



AR 



V(Ar?)2 + (A</.)2 



V (??true - Tlreco)^ + (0true " ^/^'reco)^ 
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and for the different trigger levels, accordingly. 

A positive match of the given objects is denoted by a AR below a predefined threshold. This 
threshold is typically chosen according to the spatial resolution of the measurement for the 
respective objects. 

Both trigger and reconstruction efficiencies are commonly taken into account in physics anal- 
yses using scale factors, relating the respective efficiencies measured with data driven methods 
to the expected performance from Monte Carlo simulations. Alternatively, in particular trigger 
efficiencies can be applied by implementing a reweighting approach, taking into account detec- 
tor and trigger acceptance by direct application of data driven efficiencies in the form of event 
weights. 

In this context, the reconstruction efficiency ffreco is typically defined as the fraction of objects 
which have been identified by a specific reconstruction algorithm, 

-' ''candidate 

where A/'candidatc denotes the amount of reconstruction candidates (e.g. tracks or energy deposi- 
tions) considered for reconstruction. 

In analogy, the trigger efficiency etrig refers to the fraction of reconstructed objects which have 
been selected by a given trigger item or chain, 

-tng = f^. (4.3) 
"^reco 

These efficiencies are typically parametrised or binned in kinematic quantities of the given ob- 
jects, such as pt and rj. 



4.2.1. Measurement of Trigger Efficiencies 

In order to obtain trigger efficiencies from data, a simple Monte Carlo based counting method 
is not applicable as events rejected by the trigger are typically no longer available offline. Fur- 
thermore, it is not desirable to rely on the trigger simulation. To measure trigger efficiencies, 
there exist several data driven methods, including: 

• orthogonal triggers: the trigger efficiency e\ for a given trigger item A is determined 
with respect to a different, ideally orthogonal (i.e. completely uncorrelated) trigger B, 
whose efficiency sb is known, 

• minimum bias datasets: the trigger efficiency for an arbitrary trigger item is determined 
on a minimum bias dataset, where only a minimal trigger selection has been applied and 
thus a trigger efficiency measurement can be performed with minimal bias, 

• the Tag & Probe method: the trigger efficiency for a given trigger is measured by 
selecting a specific process from a given data sample (if possible) of which the kinematics 
are known, allowing for a decoupling of event selection and trigger efficiency determination. 

In the following, the Tag & Probe method is explained in more detail since it is commonly used 
to determine both trigger and reconstruction efficiencies for physics analyses. 

To measure the trigger efficiency from data, it is possible to select a specific process, e.g. 
Z/^* — )■ /i/i events for the determination of muon trigger efficiencies by application of the Tag 
& Probe data method, which is shown schematically in Figure HTTJ The selection is achieved by 
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Figure 4-1 ■■ The Tag & Probe method. An event selection is performed using the tag object, 
while the actual trigger efficiency is determined only with the probe object to 
ensure independence of the two processes. 



making use of the fact that if one isolated muon is present in a Z/^* — t- /i/i event, there has to 
be a second muon in the observed event that is isolated as well, on which the actual efficiency 
measurement can be performed. The basic concept of this method focuses on decoupling event 
selection and the actual determination of the trigger efficiency. 

Events are selected by first identifying a muon that has triggered the event, the tag muon. 
Moreover, several requirements are applied to the tag muon, e.g. and isolation criteria. To 
ensure the Z/7* — )• sample to be as pure as possible, these requirements should be very tight, 
leading to an enrichment of Z/^* — t- //// events. 

If the tag muon meets all requirements, a second muon (if present) in the event is selected, 
the probe muon, which must comply to requirements that typically correspond to the selection 
performed in the analysis for which the obtained efficiencies are used. 

In addition, the dimuon invariant mass M^^ has to be sufficiently close to the Z pole mass. 
If that is the case, there is a high probability that the probe muon does not originate from a 
background process, but constitutes a muon from the Z/'j* — )• fifi signal. 

The actual trigger efficiency stp for the Tag &: Probe method can be determined by the 
fraction of probe muons that have been triggered: 



N- 



probe & trigger 



^"^P Tijprobe ' ('^•'^) 

Since the invariant mass constraint directly relates the tag and the probe muon, it is obvious 
that the geometric correlation can create a bias in the estimated trigger efficiency as well. This 
can, for example, be the case if one of the two muons produced in the Z decay emerges into an 
acceptance gap of the detector, while the other muon is properly reconstructed. In the context 
of the Tag & Probe method, this muon could pass the tag muon criteria, but since there is 
no second muon in the event, no efficiency can be determined, the event not being taken into 
account at all. Hence, the Tag &: Probe method does not cover events where only one of the 
muons was found, neglecting the corresponding contribution to the overall efficiency. This effect. 
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however, is typically very small. 

4.2.2. Scale Factors and Event Rejection Approach 

If the data driven efficiency edata for a specific trigger item or trigger chain has been measured 
and the corresponding Monte Carlo trigger simulation efficiency emc is known, the respective 
scale factor / is defined by 

/ = (4.5) 

which quantifies the discrepancies between measured and predicted trigger or reconstruction 
efficiency. Typically, both the underlying efficiencies and scale factors are parametrised accord- 
ing to actual dependencies on geometric and/or analysis related quantities such as transverse 
momentum or pseudorapidity. 

The scale factors are applied by rejecting events on MC which do not pass the trigger criteria 
of the trigger simulation and applying the scale factors as an event weight based on all relevant 
objects on MC events passing the trigger simulation. For events requiring exactly one recon- 
structed object, the event weight equals the scale factor corresponding to the respective object 
properties. However, if more than one object is required in the final state, the probability p of 
an event to be triggered is determined by the logical OR of the trigger probabilities of the all 
N objects taken into account, 

N 

P=l-ll{l-e,), (4.6) 

1=1 

where corresponds to the trigger efficiency/probability of the i-th object. Correlations can 
be taken into account by an appropriate parametrisation of the underlying efficiencies and scale 
factors. The corresponding scale factor for events requiring more than one final state object is 
given by 

1 _ n ( 1 - edata\ 
/ = . U^.. !uJ - (4-7) 

i 



i-n(i-^ 



Consequently, this scale factor depends on the individual object properties, even if the respective 
efficiencies e^^'^ and e^*~' do not, since Equation 14.71 does not factorise in these quantities. Due 
to limited statistics, taking into account all possible combinations of parameters for objects 
necessary for an accurate representation of scale factors, is not feasible. In particular, the 
treatment of correlations in the statistical uncertainties of the scale factors becomes increasingly 
difficult both analytically and computationally. 

Furthermore, depending on the magnitude of the underlying efficiencies, a significant amount 
of MC events is rejected when requiring the trigger simulation to have selected events, leading to 
a reduction of available MC statistics and practical loss of CPU time spent on the generation of 
the rejected events. This can be of great importance, in particular if analyses are limited by MC 
statistics. However, since efficiencies of some of the relevant objects, such as electrons, usually 
are relatively high (> 95%) in the used signal regions, the corresponding loss of MC statistics 
can be negligible for those objects, depending on the amount of available generated events. 

4.2.3. Reweighting Approach 

As an alternative to the determination and application of scale factors, a trigger reweighting 
approach, making exclusive use of data driven measurements of trigger efficiencies, can be applied 
to MC events. If the data driven efficiency ffdata for a specific trigger item or trigger chain is 
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known, the corresponding event weights can be determined regardless of the amount of required 
final state objects. The event weight w is calculated directly from the respective parametrisation 
of the object efficiencies: 

TV 

w = l-ll[l-ef^'^), (4.8) 

i=l 

which factorises in the trigger (in)efficiencies. 

Note that only objects following identical selection criteria as used in the efficiency determi- 
nation should be taken into account in the event weight calculation to ensure that the respective 
parametrisation is valid. Furthermore, the choice of a proper parametrisation of the efficiencies 
is crucial due to the fact that neglecting potential dependencies of the efficiencies can lead to 
significant mis-modelling of the trigger response on MC. This effect is typically smaller for the 
scale factor approach, where the weights are typically close to one even if the overall efficiencies 
significantly differ from one (and hence, the impact of the parametrisation is smaller). 

However, the trigger reweighting approach offers several additional advantages over a scale 
factor approach, most prominently due to the fact that no trigger simulation is involved in 
either the determination or application of the trigger efficiencies. Since the simulated trigger 
decision is not taken into account, all MC events are preserved regardless of the magnitude of the 
underlying trigger efficiencies, retaining the full MC statistics. In addition, efficiencies obtained 
from data can be applied to MC independent of potential changes or errors in the implementation 
of the trigger simulation, which would make repeated extraction of scale factors or alternative 
solutions necessary. Since the trigger requirements for a given data period typically do not 
change, the obtained efficiencies have to be determined only once for a set of reconstruction 
algorithm parameters and remain unchanged as well. 

Possible correlations in the statistical uncertainties derived for the object efficiencies can be 
modelled properly in a relatively simple way due to the fact that all objects obtain efficiencies 
from the same parametrisation. Hence, a propagation of the statistical uncertainties to the final 
event weights, event yields or even distributions of object and event quantities is possible. How- 
ever, since the uncertainties on trigger efficiencies are usually relatively small compared to other 
sources of uncertainties in physics analyses, the impact is expected to be small when neglecting 
the respective correlations. Typically, uncertainties on both scale factors and efficiencies are 
regarded as fully correlated in a conservative approximation. 

4.3. Objects 

In the following, a summary of objects used in this analysis and their definition will be given. 
The algorithms used to reconstruct the individual objects and the corresponding performance 
will be discussed. 

4.3.1. Jets 

Jets represent collimated collections of long-lived hadrons, originating mostly from partons after 
the fragmentation and hadronisation process |164| . Jets were reconstructed by associating ob- 
jects from the EM calorimeter and hCAL, and tracks reconstructed in the detector into physical 
jet objects representing the underlying fragmentation and parton shower processes. 

Jet reconstruction algorithms are preferably collinear and infra-red safe, i.e. collinear splitting 
and soft gluon emissions should not change the algorithm response of the final reconstructed 
jet. In this context, the jet reconstruction was performed using the anti-Zcj (R = 0.4) jet 
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clustering algorithm [165^166] . Tracks and energy depositions {topological clusters |145j ) in the 
calorimeter identified at the electromagnetic energy scale |145| (EM scale), calibrated to yield a 
correct energy response for electrons and photons, were associated into combined jets of particles. 
This was achieved by iterative combination of pairs of objects into proto-jets according to the 
respective geometric distance of the individual constituents and the proto-jets until a convergent 
state was reached. 

Jet quality criteria were applied to identify jets which do not correspond to physical in-time 
energy deposits in the calorimeter. Possible sources for such bad jets are for example hardware 
problems, calorimeter showers induced by cosmic rays and beam remnants. 

The jet reconstruction efficiencies were determined from data with a Tag & Probe method, 
using jets from charged tracks in the ID, where the efficiency was defined as the fraction of 
probe track jets matching a corresponding calorimeter jet. The jet energy resolution a{Ei)/ Et 
for jets in the region |Ty| < 2.8 ranged from about 0.2 (for 20GeV jets) to 0.05 (for 900 GeV 
jet^j) |167| . Differences between measured efficiencies and efficiencies obtained in Monte Carlo 
simulations were taken into account by randomly removing jets from events according to a 
Gaussian probability distribution corresponding to the jet reconstruction efficiency uncertainty. 

In addition, a calibration of the jet energy scale |145j (JES) was performed to obtain the 
energy of the final state particle- level jet from the measured energy of the reconstructed jet. 
The differences in response from the EM calorimeter and hCAL and additional effects such as 
energy mis-measurements of particles not coming to a stop in the calorimeter were taken into 
account. 

The following requirements were applied to reconstructed jets (calibrated at the electromag- 
netic + JES (EM+JES) scale unless stated otherwise) after electron/muon overlap removal: 

• Jet transverse momentum: PT(jet) > 25 GeV, 

• Jet pseudorapidity at the EM scale, including JES r] correction: |??EMscaio(jet)| < 2.5. 
b Tagging 

Due to the dominant decay of top quarks into a W boson and a b quark, events containing top 
quarks are characterised by hard b jets, which have the distinction of possessing a long lifetime, 
a large corresponding B hadron mass and a large branching ratio into leptons. Hence, the 
identification of b jets and the discrimination against light quark jets can significantly increase 
the signal to background ratio for top quark decays. 

The identification (or tagging) of b jets is mostly driven by exploiting the increased lifetime 
with respect to light jets and the associated significant fiight path length I, leading to secondary 
vertices and measurable transversal and longitudinal impact parameters do and zq of the fi- 
nal state particles. These denote the spatial distance of closest approach with respect to the 
associated primary vertex in the transversal respectively longitudinal plane. 

For this analysis, the JetFitterCombNN tagging algorithm |168| was used, which determines 
a b tag probability/weight w for a given jet according to a Neural Network combination of the 
weights from the IP3D and JetFitter tagging algorithms [168j . The IP3D algorithm used the 
significance of dQ and zq of each track contained in the respective jet to determine a likelihood 
corresponding to the b jet probability. The JetFitter algorithm implemented a Kalman Filter 
|169j to identify a common line of primary vertex and weak b and c hadron decay vertices within 
the jets. Using the approximated hadron fiight path, a b jet likelihood was obtained based on 
the masses, momenta, fiight length and track multiplicities of the reconstructed vertices. 

^Updated plots at https://twiki.cern.ch/twiki/bin/view/AtlasPublic/JetEtmissApproved2011JetResolution 
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A cut on the b tag weight oi w > 0.35 was apphed to classify b jets, corresponding to an 
approximate b tagging efficiency of 70% in simulated tt decays and a rejection rate of about 
5 for c jets and about 100 for light flavour jets. The performance of the b tagging algorithm 
was determined on specific data samples and compared to the MC prediction. A set of scale 
factors, parametrised in jet pT was determined and applied to MC events for both the b tagging 
efficiency and the light jet mis-tag rate to take into account differences between data and MC. 
The respective event weight was determined as logical OR from the scale factors of all jets taken 
into account. 

4.3.2. Muons 

Muons were identified and selected both at online (trigger) and offiine (reconstruction) level in 
accordance with the respective object quality criteria. 

Muon reconstruction was performed with the Muld |170j algorithm, taking into account the 
tracks of the ID and the MS and creating a combined reconstructed track. Particle hits identified 
in the MS chambers were associated to the corresponding track segments in the individual layers 
of the ID. A combined fit was performed to obtain an optimised reconstructed trajectory from 
the joint information of both systems and to optimise geometrical and transverse momentum 
resolution. The MS momentum resolution in o"(pt)/pt ranges from 0.04 (for 20GeV muons) 
to 0.06 (for 200 GeV muons) in the central detector region and from 0.05 (for 20 GeV muons) 
to 0.20 (for 200 GeV muons) in the forward detector region. The ID momentum resolution in 
a{p^^) ranges from 0.001 (for 20 GeV muons) to 0.0005 (for 200 GeV muons) in the central 
detector region and from 0.004 (for 20 GeV muons) to 0.005 (for 200 GeV muons) in the forward 
detector region |171j . 

In addition, differences in the muon momentum resolution between data and Monte Carlo 
were corrected for by randomly changing the muon pT on Monte Carlo according to a Gaussian 
resolution function to reflect the resolution observed in data. 

Muon trigger and reconstruction efficiencies were measured with data driven methods on 
Z/'j* — )■ /i/i events using the Tag & Probe method and were compared to Monte Carlo simula- 
tions. Both the respective efficiencies and scale factors were parametrised in muon transverse 
momentum, pseudorapidity and azimuthal angle to reflect differences in the detector acceptance 
and the px dependency for both the trigger and reconstruction efficiencies. 

Muons were selected at trigger level by requiring the EF_mul8 trigger chain (corresponding 
to a muon with a > 18 GeV at trigger level), which was seeded at LVLl by the L1_MU10 
trigger item and at LVL2 by the L2-mul8 trigger item, to have fired. This requirement was 
applied to data only, while the trigger efficiency was taken into account for Monte Carlo events 
by determining an event weight based directly on the trigger efficiency associated with the re- 
constructed and selected muons. Hence, no explicit (simulated) trigger decision was required. 
Furthermore, no geometric matching of trigger objects and reconstructed muons has been ap- 
plied due to a software error resulting in incorrect trigger object representation within the used 
data and Monte Carlo samples. The resulting mis-modelling of the trigger response has been 
estimated and was taken into account as an additional systematic uncertainty (c.f. Chapter [9]). 

Reconstructed muons from the Muld algorithm were required to be combined muons, i.e. to 
have a combined ID and MS track passing the respective track criteria. In addition, cuts specific 
to select semileptonic top quark decays were applied in order to reject muons from heavy and 
light flavour decays such as b and c hadrons or in-flight decays of kaons and pions. In particular, 
isolation requirements were applied based on the momentum and energy deposition within a 
Ai? cone of size 0.3 around the muon track, i?T,conc30 (a*) and PT,cone3o(/^)) respectively. The 
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following requirements have been used: 

• Muon transverse momentum: pt{i-i) > 20GeV, 

• Muon pseudorapidity: \r]{fj,)\ < 2.5, 

• Jet overlap removal by requiring AR{fi, closest reconstructed jet) > 0.4, where jets recon- 
structed by the anti-/cj algorithm with R = 0.4 were taken into account, calibrated to the 
EM+JES scale and with pT(jet) > 20GeV, 

• Isolation: -Et.coiicSoIm) < 4GeV and PT,conc3o(^) < 4GeV. 
4.3.3. Electrons 

Similarly to muons, electrons were identified and selected according to the respective electron 
quality criteria. 

Electrons were reconstructed by matching energy cluster hits (or seeds) above a threshold of 
about 3 GeV in the electromagnetic calorimeter to corresponding extrapolated ID tracks, vetoing 
tracks from photon conversion pairs. In addition, a matching of the track momentum and cluster 
energy was performed by application of a cut on E/p. Information from the TRT chambers was 
used to enhance the separation of electron candidates from pions in the reconstruction process. 
The electron energy resolution was about (1.2 it 0.1 (stat.) ±0.3 (syst.)) % in the central detector 
region, about 1.8 % in the endcaps, and of the order of 3% in the forward detector region |172j . 

The electron energy scale and resolution were obtained with data driven methods in kinematic 
regions similar to top quark pair production events, using events from Z/'j* — )• ee decays. On 
data, the energy scale was corrected as a function of electron transverse energy Eic{e) and cluster 
pseudorapidity Irydusterl- On MC, the electron energy was corrected by applying a randomised 
Gaussian resolution function to the electron energy, following a similar parametrisation. Statis- 
tical and systematic uncertainties on both scale and energy corrections were taken into account 
and were assigned to the MC. 

Electron trigger and reconstruction efficiencies were measured with data driven methods on 
Z/^* —7- ee and — )• eug events using the Tag &: Probe method and compared to Monte Carlo 
simulations. The respective electron scale factors were parametrised in electron pseudorapidity 
I fluster I (using cluster quantities) to reflect differences in the detector acceptance for both the 
trigger and reconstruction efficiencies. 

Electrons were selected at trigger level by requiring the EF_e20-medium trigger chain (corre- 
sponding to an electron with a i?T > 20 GeV at trigger level, and medium denoting a specific 
set of object criteria at trigger level, respectively), which was seeded at LVLl by the L1_EM14 
trigger item and at LVL2 by the L2-e20-medium trigger item, to have fired. This requirement 
was applied to both data and MC events. In addition, a geometric matching of reconstructed 
electrons to the corresponding trigger objects in a given event was performed to ensure that 
the reconstructed electron had fired the trigger. Data/MC scale factors were applied as event 
weights to take into account differences in the efficiencies between data and trigger simulation. 

In this analysis, reconstructed electrons were required to pass the tight |173j electron quality 
requirements. In addition, the following requirements were applied to electron candidates: 

• Electron transverse energy: -E'T(e) > 25 GeV, 

• Electron cluster pseudorapidity: |??ciuster(e)| < 2.47, 

• Exclusion of the non-covered transition region 1.37 < |??ciuster(e)| < 1.52, 
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• Jet overlap removal by requiring Ai?(e, closest reconstructed jet) > 0.2, where all jets 
reconstructed by the anti-fet algorithm with i? = 0.4 were taken into account, 

• Isolation: -ET,conc2o(e) < 3.5 GeV after leakage and pile-up correction |173j . 

where -ET,cone20 (e) denotes the energy deposition within a AR cone of size 0.2 around the electron 
cluster. The electron transverse energy £'T(e) was determined from the energy deposited in the 
cluster in the calorimeter -Eciuster(e) and the associated track pseudorapidity r7track(e), 

S^(e) = -^"'"'^^"■•(^^ (4.9) 

cosh (7?track(e)) 



4.3.4. Missing Transverse Energy 

The missing transverse energy is an object-based quantity derived from the topological 
clusters in the calorimeter calibrated at the EM scale, corrected for the energy scale of the 
corresponding object associated to the cluster. Objects taken into account were electrons and 
jets, where jets were separated into high- and low-pT (soft) jets. The electron contribution 
used electrons passing the tight electron quality criteria with > 10 GeV, while jets with 
PT > 20 GeV were corrected to the EM+JES scale and soft jets with 7 GeV < < 20 GeV 
were included at the EM scale. 

In addition, muons were included in the definition using the transverse momentum of the 
corresponding tracks due to the fact that muons typically only deposit small amounts of energy 
in the calorimeter. Muons reconstructed with the Muld reconstruction algorithm with |r/| < 2.5 
were taken into account, distinguishing between isolated (requiring Ai2(//, closest jet) > 0.3, 
where all jets reconstructed by the anti-A;t algorithm with R = 0.4 were taken into account) 
and non-isolated (where the energy deposited in the calorimeter was taken into account in 
the jet term). Furthermore, muons in detector regions with low acceptance (|r/| < 0.1 and 
1.0 < \t]\ < 1.3) were taken into account making use of the corresponding calorimeter response. 

Clusters not associated to any object were included in an additional (CellOut) term, calibrated 
at the EM scale, as well as the energy deposition of isolated muons in the calorimeter. 

The missing transverse energy terms in the x and y directions were calibrated and combined 
into an overall missing energy in the respective dimension: 

^electrons _|_ ^j^ts _|_ ^softjcts _|_ ^muons _|_ ^CellOut j^Q^ 

where i = {x,y}. Consequently, the respective scalar transverse missing energy is given by 

= ^/m+^y■ (4.11) 

The resolution of the and components ranged from about 2 GeV (for a total transverse 
energy of 20 GeV) to 14 GeV (for a total transverse energy of 700 GeV) |174j . 
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The event selection used for the top charge asymmetry measurement was aimed at maximising 
the signal contribution from the tt decay in the semileptonic decay channel. At the same time, 
the background contribution from different sources was minimised. The two distinct observable 
final states, muon+jets and electron+jets, were treated independently. 

As discussed in Chapter [51 top quark pairs are dominantly produced at the LHC by gluon- 
gluon fusion. Both the top and antitop quark then decay into a W boson and a b quark in 
almost 100 % of all cases. The W boson decays into two jets or a charged lepton and a neutrino, 
which can only be measured indirectly as missing transverse energy. Hence, the event selection 
for the semileptonic decay channel was focused on topologies with at least four reconstructed 
jets, exactly one isolated lepton (muon or electron) and missing transverse energy. 

A preselection of the delivered raw data based on a common GoodRunsList (c.f. Chapter 14. ip . 
ensuring stable beam conditions and data quality, has been applied prior to the event selection. 
This selection included global data quality flags, e.g. requiring stable beams at a centre-of-mass 
energy of 7TeV and the LVLl central trigger and luminosity measurement to be functional, 
indicating stable running and data taking conditions of the LHC and ATLAS, respectively. 
Furthermore, the data quality of the individual detector subsystems has been verified. These 
criteria correspond to a data quality selection efficiency of 84.1% in the electron+jets channel 
and 84.3 % in the muon+jets channel. 

In addition to the baseline selection, which will be described in the following, several cor- 
rections have been applied to Monte Carlo samples on an event-by-event basis to account for 
potential mis-matches in the detector simulation with respect to data, as described in Chapter [H 
Muon trigger efficiencies have been taken into account directly by performing a trigger reweight- 
ing using trigger efficiencies obtained with a Tag & Probe method, while discrepancies between 
the simulation and data for electron and muon reconstruction efficiencies and the electron trig- 
ger efficiencies as well as the b tagging efficiencies have been taken into account by applying the 
corresponding scale factors to the Monte Carlo events. Furthermore, each Monte Carlo event has 
been assigned a weight according to the average amount of pp interactions per bunch crossing 
in the respective sample to account for differences in the modelling of pile-up. 

The following event selection has been used to enhance the signal to background ratio in the 
recorded samples. The described selections have been applied to both data and Monte Carlo 
samples unless stated otherwise: 

• The electron or muon trigger was required to have fired. The trigger item used was 
EF_e20 -medium in the electron+jets channel for both Monte Carlo and data. In the 
muon+jets channel, the EF_mul8 trigger item was required for data, while a trigger 
reweighting using EF_mul8 trigger efficiencies obtained with data driven methods was 
applied to Monte Carlo events, made necessary by a mis-modelling of the trigger simula- 
tion in the used Monte Carlo samples. 

• A primary vertex with at least five tracks associated to it was required to improve rejection 
of non-collision background from the underlying event, pile- up and cosmic radiation. 
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• Exactly one isolated lepton (one electron and no muon or vice versa) passing the respective 
object quality criteria with > 25GeV (electron +jets) or > 20GeV (muon+jets), 
respectively, was required. 

• The selected lepton was required to match the object of the fired trigger (electron+jets 
only, since due to a software problem, the muon trigger matching requirement was dropped 
from the selection). 

• Any event where a reconstructed electron and muon share a common track was rejected. 

• Any event containing a bad jet (c.f. Chapter 14.3. ip with pT > 20 GeV at the EM+ JES 
scale was rejected. 

• A missing transverse energy > 35 GeV (electron+jets) or > 20 GeV (muon+jets) 
was required. 

• A cut of 

mT(VF)0> 25 GeV (electron+jets) or a triangular cut of + mT{W) > 60 GeV 
(muon+jets) was applied in order to suppress the QCD multijet background contribution, 
since these events typically have low rriTiW) and low ^t- 

• At least four jets with pT > 25 GeV passing the jet quality criteria were required. 

• Any event where a jet was found in an area of LAr calorimeter defects was rejected and 
electrons which were affected were removed from the respective events, correcting the 
measured fJ^ accordingly. For Monte Carlo, a randomised subset of events was dropped 
according to the relative fraction of the data sample affected by these defects in order to 
correct for the created mis-match between data and simulation. 

• At least one jet which has been b tagged using the JetFitterCombinedNN algorithm with a 
weight w > 0.35 (corresponding to an overall b tagging efficiency of about 70 % in simulated 
tt events) was required to further improve the signal to background ratio. 



'^In this context, mT{W) denotes the W boson transverse mass, defined as 




(5.1) 



where pf^" and (j)'^'^ describe the lepton and neutrino transverse momentum and azimuthal angle, respectively. 
The neutrino information is represented by the measured missing transverse energy, ^t- 
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This chapter describes the data sample used in this analysis and the Monte Carlo samples 
generated to simulate the signal contribution and most of the backgrounds. Furthermore, the 
data driven estimations of the dominant background contributions, VF+jets and QCD multijet 
production, are described. 

6.1. Data Sample 

A set of ATLAS data taken in the course of the year 2011, corresponding to an integrated 
luminosity of 



has been analysed, after the preselection of the delivered raw data using the corresponding 
GoodRunsList. The data has been recorded between March, 22^'^, 2011 and June, 28*^, 2011. A 
peak instantaneous luminosity of about C = 1.3- 10^^ cm~^s~^ was reached and a bunch spacing 
of 50 ns was used. 

6.2. Signal and Background Monte Carlo Samples 

Several Monte Carlo samples have been generated to facilitate this analysis, including nominal 
samples for the signal contribution and various background processes as well as several additional 
samples used in the evaluation of systematic uncertainties. All samples correspond to the mclOh 
production commonly performed for all ATLAS analyses using a generalised set of parameters 
to match the data taking conditions during the time period that was considered. 

In particular, the contribution from in-time and out-of-time pile-up was added to all generated 
Monte Carlo events in the parton showering simulation process after the generation of the initial 
hard scattering. A fixed configuration of average proton-proton interactions per bunch crossing 
was used, while the actual data taking conditions with respect to pile-up changed over the course 
of time. In order to correct for the mis-match between individual data taking periods and the 
simulated pile-up contribution in the MC, an event-based reweighting was performed, taking 
into account the expected and observed distribution of the average number of interactions per 
bunch crossing. This pile-up reweighting was performed for all MC samples. 

The ti signal process has been simulated using the mc@nlo generator |175j (v3.41) which in- 
corporates the CTEQ6.6 |176| parton distribution function set and makes use of a next-to-leading 
order calculation approach for QCD processes. Both the parton showering and fragmentation 
processes, and the underlying event have been modelled using the HERWIG v6.510 |177] and 
JIMMY dZH] generators utihsing the CTEQ6.6 and AUETl [179] tunes to match the ATLAS 
data, respectively. The inclusive ti cross section has been estimated to approximately next-to- 
next-to-leading order using the Hathor tool |180j to be 165^Jg pb [36H38] for nit = 172.5 GeV 
and the MC has been scaled accordingly. For this analysis, only semileptonic decays of the 
top quark pairs were considered. The respective cross-section, taking into account the proper 
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branching fractions, was 89.3 pb, including a fe-factor of 1.117 to rescale the next-to-leading 
order perturbative QCD cross section in MC@NLO to the approximate next-to-next-to-leading 
order cross-section. The signal sample contained 15 000 000 simulated events, corresponding to 
an integrated luminosity of about 150 fb~^. 

The electroweak single top production was simulated using the mc@nlo and jimmy generators 
and the respective cross-sections have been calculated at approximately next-to-next-to-leading 
order to be 64.57110^ pb in the t-channel, 4.63+o;i7 pb in the s-channel, and 15.74^]^Qg pb for 
Wt production, as introduced in Chapter 12.2.21 

The background contribution from the production of heavy gauge bosons with additional 
jets was modelled using the leading order alpgen generator |181j . interfaced to herwig and 
JIMMY for the purposes of parton shower and hadronisation simulation. The CTEQ6.1 |182j 
parton distribution functions and the AUETl tune were employed for proper ATLAS data 
matching for both the matrix element evaluations and the parton showering. The production of 
additional partons was taken into account by generating different subsamples with different final 
state parton multiplicities, where additional partons can be either light {u,d,s) partons (simu- 
lated in VF+jets and Z/7+jets light flavour samples) or heavy quarks (simulated in l^+c+jets, 
W+cc+iets, W+bb+jets, and Z+66+jets samples, respectively). Since the inclusive VF+jets 
and Z/7+jets samples included contributions from both light partons and heavy quarks in the 
matrix element and parton shower simulation, the created overlap in phase space between the 
inclusive samples and the heavy quark contribution was taken into account by removing double 
counted events from the respective samples. The production cross-sections of the used alp- 
gen samples were normalised to the corresponding approximate next-to-next-to-leading order 
cross-sections using A;-factors of 1.20 (Vl^+jets) and 1.25 (Z+jets), respectively. Furthermore, 
the relative fractions of the individual W^+jets heavy quark contributions to the overall VF+jets 
sample have been determined in data driven methods [18311184] . and were accounted for by 
applying corresponding scale factors to the individual samples. They have been found to be 
1.63 ± 0.76 for W+cc+jets and l^+66-|-jets, and 1.11 ± 0.35 for the M^+c+jets contribution, 
respectively. The VF+jets light quark contributions were scaled accordingly to conserve the 
overall predicted cross-section. 

Contributions from diboson (VFVl^, WZ and ZZ) production and decays was simulated using 
herwig at leading order, and the corresponding production cross-sections were normalised to 
the next-to-next-to leading order predictions, using A;-factors of 1.48 (VFW^), 1.60 (WZ) and 
1.30 {ZZ), respectively. Each sample was inclusive and has been filtered to include only events 
containing at least one lepton (electron or muon) with > lOGeV and |?7| < 2.8 at parton 
level. The /c-factors have been determined such that the unfiltered herwig cross sections agree 
with the next-to-next-to-leading order calculations. 

For the evaluation of the systematic uncertainties of the various generators and the simu- 
lation of hadronisation and fragmentation, alternative samples have been used for the signal 
contribution and the Z+jets background contribution. The tt production and decay has been 
simulated with the powheg generator |185j . and the corresponding parton showering and frag- 
mentation processes have been modelled using both herwig and jimmy (as used for the nominal 
MC@NLO sample), and pythia in order to evaluate systematic effects from the parton show- 
ering. Furthermore, additional MC@NLO signal samples using different top mass hypotheses of 
170 GeV and 180 GeV have been used in order to quantify systematics arising from the uncer- 
tainties on the top mass prediction. Finally, several samples using different strengths of initial 
state radiation (ISR) and final state radiation (FSR) have been generated at leading order using 
the acermc generator |186j . corresponding to different contributions of ISR and FSR based on 
observations in data |187| . An alternative inclusive Z+jets background modelling was performed 
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using the SHERPA |188j generator and the CTEQ6.6 parton distribution functions. Details on 
the evaluation of systematic uncertainties can be found in Chapter [9l 

All background samples and the samples used for the evaluation of systematic uncertainties 
corresponded to an integrated luminosity of about 10 — 30fb~^ before analysis specific selection. 

6.3. Data Driven Estimation of the QCD Multijet Contribution 

The identification of top quark pairs decaying semileptonically relies on the identification of one 
lepton in the final state, carrying a large transverse momentum. Hence, mis-identified leptons 
(fake leptons), which can originate from various sources, pose a non-negligible background to 
the identification of tt signal events. Potential sources for mis-identified leptons include 

• semileptonic decays of b quarks into b jets containing leptons with mis-identified isolation 
properties, 

• long lived weakly decaying particles such as vr^ or K mesons, 

• vr^ mesons which are mis-identified as electrons, 

• direct photon conversion and reconstruction of electrons produced in the process. 

These processes are most dominant in regions where the contribution from real leptons are small, 
most prominently for the background contribution from QCD multijets. 

Despite the fact that object and event selection are designed to ensure the suppression of 
these backgrounds by requiring stringent criteria, the QCD multijet production cross-section 
is orders of magnitudes higher than the top quark pair production cross-section. Since the 
simulation of these backgrounds in Monte Carlo simulations is highly difficult and several of the 
described contributions are detector dependent, data driven methods are necessary to obtain 
reliable estimates for the fake lepton background contribution from QCD multijet events. 

In order to estimate the contribution from QCD multijet fake muons and electrons, a data 
driven method, the Matrix Method, was applied to data using different control regions for 
electrons and muons dominated by QCD multijet processes. The Matrix Method allows to sta- 
tistically separate two contributions of a data sample based on the impact of a defined selection 
on the respective subsamples. 

The Matrix Method defines two subsets, A'^ioosc and A'tight) of the data sample before and after 
application of a particular requirement applied in the event selection, typically with a large 
discrimination power between signal (in this case real leptons) and background (in this case fake 
leptons) contributions. The number of events in the loose sample, A^ioosei is given by the sum of 
the signal and the background contributions N^^^ and TV^'^kc -^^ given sample: 



After the application of the defined requirement, which is passed by signal events with a proba- 
bility (or efficiency) of e**'^ and by background events with a probability of e^'^'^°, the number of 
events A'^tight in the tight sample follows by taking into account the respective probabilities for 
the imposed requirement for the signal and background contributions: 



Moose = N''^ + AT' 



fake 




fake 



(6.3) 



This situation is illustrated in Figure 16. 1[ 



53 



6. Samples and Process Modelling 




Figure 6.1.: An illustration of the Matrix Method and the effect of the applied selection on 
the underlying sample subsets. 



The linear system of two equations with two unknown variables can be rewritten as a matrix 
equation: 

' Moose \ f I I \ f iV"^ 



Solving the matrix equations yields the signal and background contributions in the data sample 
prior to the isolation cut: 



g-iakc gSig ' ^ ' 

jysig _ ^^^^'^^loose — -Might 

g-fake gsig ' \ ' J 



In order to determine the background contribution in the signal region (which corresponds to 
the tight selection by construction), the fraction of A^f^^° in the tight sample can be calculated 
using Equation 16.51 as 



AT-fake fake /vrfake 

^^tight — e 

g.fake 



fake 



(e'^^Moosc - Might) • (6.7) 



If the selection probabilities e^^^ and e^^^^ for signal and background, respectively, are sufficiently 
different, the overall contribution of the QCD multijet background can be used to determine 
event based weights for the used data sample in order to obtain the distributions of the QCD 
multijet background contribution in arbitrary variables. This is done by assigning a weight 
to each data event based on the chosen requirement and the corresponding signal and fake 
probabilities of the objects taken into account for a given event. If it passes the loose selection 
only, the event weight is given by setting Moose = 1 > Might = in Equation 16.71 yielding 

g-siggfake 

"^^oose = ^sig _ ^fake " (6-^) 

Similarly, if both the loose and tight requirements are fulfilled, the event weight is given by 
setting Moose = 1, Might = 1 in Equation 16.71 



(^gsig _ l) ^fake 



^tight = ,si. .fake • (6-9) 



This approach allows for a purely data driven estimation of both the normalisation and the 
shape of the QCD multijet background in semileptonic decays of top quark pairs. The individual 
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parameters that have been used in the estimation for both the muon and electron channels in this 
analysis will be covered in the following sections. In addition, detailed studies of the performance 
and stability of the methods used will be demonstrated for the muon+jets channel alongside 
with an approach to obtain well-defined statistical and systematic uncertainties on the estimate. 



6.3.1. Muon+jets Channel 



In order to estimate the contribution from QCD multijet fake muons, the matrix method was ap- 
plied to data in the QCD multijet-enriched low-mT(VF) control region. Furthermore, an inverted 
triangular cut was imposed to achieve orthogonality to the signal region in the determination of 
the fake probabilities e^^^'^ used in the Matrix Method: 

mT{W) < 20GeV and + mT(W) < 60GeV. 

The impact of the described requirements on the QCD multijet estimate and the simulated 
tt signal is shown in Figure | 
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Figure 6.2.: hiipact of the requirement of tot(W^) < 20 GeV and the triangular cut, 

+ m'T{W) < 60 GcV, on both the QCD multijet estimate and the tt signal con- 
tribution, taken from Monte Carlo simulations. The imposed cuts are illustrated 
by the black lines. 

The signal probabilities were determined in the signal region using a Z/'j* — t- fin Tag & 
Probe method in order to select prompt muons from Z decays. 

The loose selection was identical to the full selection applied in the signal region (for details on 
the individual requirements, refer to Chapter H] and Chapter [S]), except for the muon isolation. 
The tight selection in addition requires isolation criteria based on both momentum and energy 
depositions around the muon tracks, 

PT,conc30 < 4.0 GeV and ET,cone30 < 4.0 GeV, 
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corresponding to the full analysis event selection. 

The fake probabilities have been obtained separately both with and without explicitly re- 
quiring at least one b tagged jet (> 6 tags and > 1 b tags, respectively). In addition, the 
fake probabilities have been determined requiring at least two b tagged jets, for completeness. 
For the latter cases, the JetFitCombNN b tagging algorithm with a working point oi w = 0.35 
(corresponding to an overall b tagging efficiency of 70 %) has been used in accordance with 
the signal region event selection. Furthermore, the signal muon contribution from Pt^+jets and 
Z+jets in the control region was obtained from Monte Carlo and subtracted to obtain a purer 
QCD multijet estimation. This contamination was of the order of 1.7% (> 6 tags) and 1.8% 
(> 1 5 tags). 

In order to verify the stability of the method, the signal and fake (both without and with the 
requirement of at least one b tagged jet) probabilities are shown as a function of the relative run 
number (full dataset with GoodRunsList applied) in Figure [6^ As can be seen, neither signal 
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Figure 6.3.: Integrated signal and fake (without and with the requirement of at least one b 
tagged jet) probabilities, shown as a function of relative run number (full dataset 
with GoodRunsList applied). The observed probabilities were stable over the 
whole data taking period under consideration. 



nor fake probabilities showed any significant trend with respect to run number and therefore, to 
instantaneous luminosity and different pile-up conditions. 

In order to take into account dependencies on object kinematics, the signal and fake proba- 
bilities have been determined as a function of muon pseudorapidity r] to reflect the dependency 
on the muon detector acceptance. Furthermore, they have been parametrised as a function of 
the leading jet transverse momentum ptUi) in order to take into account the effects of hard 
jets and hence increased hadronic activity on the muon isolation. The respective projections for 
muon r] and PT(ji) can be found in Figure [631 

The event yields and fractions in the signal region, obtained on a dataset corresponding to 
1.04 fb~^ for different amounts of reconstructed jets required in the event selection (jet bins) can 
be found in Table [6Tl 

The following sources of systematic uncertainty have been taken into account in the estimate 
and were combined to a single uncertainty on the event weight, which can hence be propagated 
into a bin-by-bin normalisation and shape uncertainty on an arbitrary variable: 

• Statistical uncertainty on the signal and fake probabilities: Takes into account 
the uncertainties on both the signal and fake probabilities. The resulting uncertainty on 
the obtained event weights was evaluated using Gaussian error propagation. Note that 
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Figure 6.4-: Signal and fake probabilities as a function of muon rj and leading jet ptOi)- 



the assumption of a symmetric probability distribution function was valid for the signal 
probabilities as well despite the closeness to unity, owing to the high statistics available 
from Z/^* — )■ /iyU and the resulting small impact on the overall statistical uncertainty. 

• Systematic uncertainty from second control region: A second control region (high 
do significance) was used for the fake probabilities to determine a bin-by-bin systematic 
discrepancy which was quoted as additional uncertainty on the event weight. 

• Systematic uncertainty due to choice of control region cut: The low transverse 
W mass control region cut was varied by 5 GeV up and down to estimate the impact on 
the obtained probabilities and event weights. 

• Uncertainties on W/Z-\-jets Monte Carlo normalisation: For the 1 jet inclusive 
bin, the M^/Z-|-jets Monte Carlo normalisation uncertainties (25%) were used to quantify 
the effect on the subtraction of real leptons in the control region. 

The resulting overall normalisation uncertainties are shown in Table 16.11 for the signal region 
(> 4 jets selection). 

6.3.2. Electron-|-jets Channel 

The QCD multijet contribution in the electron-|-jets channel has been estimated in analogy to 
the muon-l-jets channel, following slightly different criteria corresponding to the loose and tight 
selections [183j . The signal probabilities used in the Matrix Method were determined in the 
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6. Samples and Process Modelling 



Jet bin 


> 6 tags 


> 1 6 tags 


> 2 b tags 


Yield 


Fraction 


Yield 


Fraction 


Yield 


Fraction 


1 jet 


36820.5 


5.6% 


3534.3 


15.2% 


0.0 




2 jet 


14915.9 


8.9% 


2633.3 


17.0 % 


106.6 


9.5% 


3 jet 


4141.8 


9.6% 


1003.8 


11.1% 


80.5 


4.0% 


4 jet 


1109.3 


8.2% 


348.7 


6.4% 


34.8 


1.6% 


4 jet (incl.) 




7.7% 


c: 1 Q q+36.2% 
oiy.o__3g 2% 


5.6% 


79 7+83.3% 
'^■-83.3% 


1.8% 



Table 6.1.: Event yields and fractions in the signal region for a dataset corresponding to 
1.04 fb~^. Uncertainties correspond to the overall normalisation uncertainty due to 
statistical uncertainties from the obtained signal and fake probabilities, a second 
control region, the systematic shift of the mT:{W) cut and the MC uncertainty on 
the W /Z+]eis normalisation used in the subtraction of real leptons in the con- 
trol region. The tagged probabilities correspond to the JetFitCombNN b tagging 
algorithm with a working point oi w = 0.35. 

signal region with a Z/^* — )• ee Tag & Probe method in order to select prompt electrons from 
the Z decay. 

The fake probabilities have been determined in the QCD multijet-enriched low-^x control 
region: 

5GeV < < 20GeV. 

Similar to the evaluation in the muon+jets channel, the loose selection was identical to the 
full selection applied in the signal region, except for the electron isolation criteria, which was 
modified in the loose selection: 

^T,cono20 < 6.0 GeV, 

while the tight selection requires 

^T,conc20 < 3.5 GeV. 

In addition, slightly less stringent track quality criteria and the corresponding definition 
were applied to the reconstructed electron candidates in the loose sample with respect to the 
tight sample. 

Furthermore, a subtraction of the real lepton contribution in the control region has been per- 
formed using the corresponding Monte Carlo predictions, in analogy to the muon+jets channel 
treatment. 

6.4. Data Driven Estimation of the VF+jets Contribution 

Due to the fact that the parton density of u quarks in the protons brought to collision in the 
LHC is on average higher than the parton density of d quarks (which can be observed already 
at lower momentum transfers, as depicted in Figure [23]) . a higher rate of than of W~ is 
expected at the LHC. Since the production rates and the asymmetry in the production of 
and W~ events has been determined to a higher theoretical precision |189yi90] than the overall 
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6.4. Data Driven Estimation of the W+jets Contribution 



Vl^+jets rate at the LHC, the observed /W~ asymmetry in data can be used to obtain an 
estimate for the rate of the VF+jets background contribution in the signal region [183] . 

Assuming that all processes other than the M^+jets production are symmetric in the final 
state lepton charge, the total number of VF+jets events, Nw, can be extracted from the amount 
of observed data events passing the selection criteria described in Chapter [S] (except for the 
requirement of at least one b tagged jet) with a positively (negatively) charged lepton, given by 
D'^ and D~ , respectively: 



The fraction 




(6.10) 



has been evaluated on Monte Carlo, based on the same event selection, and has been determined 
to be 1.56 lb 0.06 in the electron+jets channel and 1.65 ±0.08 in the muon+jets channel, respec- 
tively. The dominant contributions to the overall uncertainty were due to PDF and jet energy 
scale uncertainties, and by the uncertainties on the heavy quark contribution fractions (i.e. the 
relative contributions from Wbb+jets, M^cc+jets and VFc+jets). 

Furthermore, the obtained overall Vl^+jets rate N\y has been extrapolated to the full event 
selection by determining the relative fraction of events passing the requirement of at least one b 
tagged jet after requiring exactly two reconstructed jets, /2,>ibtag) on data, and by determining 
the ratio k2^>A of the same fraction for the sample where at least four reconstructed jets were 
required with respect to /2,>ibtagi using the T^+jets Monte Carlo prediction |184j : 

Nw,>lb tag = Nw ■ /2,>lf) tag • k2^>4- (6.12) 

The fraction f2.>ib tag has been measured to be 0.063 it 0.005 in the electron+jets channel and 
0.068 lb 0.005 in the muon+jets channel, including statistical and systematic uncertainties. The 
extrapolation factor k2^>4 has been determined to be 2.52 ± 0.36 in the electron+jets channel 
and 2.35 ± 0.34 in the muon+jets channel, respectively. The uncertainties include contributions 
from the statistical limitation of the used Monte Carlo samples and systematic uncertainties. 
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7. Kinematic Event Reconstruction 



As described in Chapter [5l the signature of a tt event in the semileptonic decay channel at leading 
order is the observation of four reconstructed jets, one isolated lepton and missing transverse 
energy. A reconstruction of the full tt final state was performed following a likelihood approach. 
A probability for the observation of a set of measured quantities under the assumption of a 
specific model and a set of model parameters was assigned. In this particular case, the model 
described the tt decay and the input quantities were the measured energies of the four jets, the 
measured energy of the lepton, and the missing transverse energy. The fit parameters of the 
likelihood were the parton energies, the lepton transverse momentum and the three neutrino 
momentum components. The likelihood was used to assign the measured jets to the decay 
products of the tt system. For this study, all permutations with four out of the five leading 
jets (if exist) were taken into account for the event reconstruction to increase the probability of 
identifying the proper combination in the presence of additional jets (e.g. from ISR or pile-up). 
Moreover, the (non-Gaussian) partonic energy resolution (the resolution of the particle jets with 
respect to the partons) of the final state objects were taken into account through the use of 
object specific transfer functions in order to obtain the final likelihood: 

L = B{Ep^i,Ep^2\mw,Tw) ■ B{Ei,E^\mw,T'w) ■ 

B{Ep^i, Ep^2, Ep^3\mt,Tt) ■ B{Ei, E^, EpA\mt,'i^t) ■ 

W{ErnP.,u) ■ WiEf^'lPy,,) ■ W(^lepl^lep) • 
4 

II W(£^jct,i|^p,i) • P{b tag I quark), (7.1) 

i=l 

where: 

• the Bs represent the Breit-Wigner parametrisation of the parton (from which the associated 
jets originated) energies Ep^i and lepton energies E'lep with respect to the fitted ones, 

• the Ws are the transfer functions associating the measured jets/leptons to the partonic 
objects, where the mapping functions of the objects are parametrised with a double Gaus- 
sian, 

• the mw and Tw denote the W boson mass and its decay width. The parameters are fixed 
to miy = 80.4 GeV and T\y = 2.1 GeV, respectively, 

• the nit and Tt denote the top quark pole mass and its decay width. The parameters were 
fixed to nit = 172.5 GeV and Tt = 1.5 GeV, respectively, 

• the X are the partonic object quantities and X their corresponding measured values, 

• P{b tag I quark) is a 6 tag probability or rejection efficiency, depending on the quark flavor. 

The probability P{b tag | quark) was used to take into account the tagging efficiency and rejec- 
tion rate of the used b tagging algorithms at a specific working point. 
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7. Kinematic Event Reconstruction 



The most probable event topology hypothesis was chosen by iterating over all possible per- 
mutations of reconstructed jets, the lepton and the missing energy and by maximising the 
logarithmic likelihood, logL. The permutation with the highest event probability was used for 
all further studies. The reconstruction efficiency, obtained in simulations, is shown for both the 
muon+jets and electron -|-jets channel in Figure [T?T1 
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(a) muon+jets channel (b) electron+jets channel 

Figure 7.1.: Reconstruction efficiencies for the muon+jets (left) and electron+jets (right) 
channel. The indicated efficiencies denote the probability of reconstructing the 
correct (or true) combination of objects (only matched events taken into account). 
The bars marked pure statistical indicate the efficiency which is expected by choos- 
ing a combination at random. 

As can be seen, the overall efficiencies for the reconstruction of the correct event topology 
{All Correct) in both channels were 62 % (74 %) with a fixed mass parameter, without (with) b 
tagging information taken into account. In order to associate the reconstructed objects with the 
corresponding truth quarks and leptons, a simple AR matching was applied, using cone sizes of 
0.3 for jets and O.f for leptons. An event was considered matched if all truth partons originating 
from the hard scattering process could successfully be identified with reconstructed jets and the 
truth lepton was matched to a reconstructed one. For the shown performance evaluation, only 
events where the four reconstructed jets and the lepton were successfully matched to correspond- 
ing truth level objects (contributing positively to the reconstruction efficiency) were taken into 
account. The matching efficiency on simulated tt events was about 30 % in both channels. 

Examples for the transfer functions used in the likelihood can be found in Figure 17.2! and 
Figure 1731 where the fit functions in different energy regions for b jets in the pseudorapidity 
range < 0.8 and for electrons in the pseudorapidity range 0.8 < \r]\ < 1.37 are shown. 

A double Gaussian function was used in the fit of the transfer functions: 

-. _(AB-pi)2 _(Ai3-p4)2 

Ty(i?truc,i?rcco) = 7^^— r(e '"''^ + P3e '"'^ ), (7.2) 

where the parameters pi, p2, Ps, pn and p^ are functions of the true energy of the respective 
particle and /S.E = Etme — -E^reco- 
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Figure 7.2.: The transfer functions mapping the measured 5 jets to the corresponding partonic 
objects, in the range \r]\ < 0.8, for different energies. 
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Figure 7. 3. : The transfer functions mapping the measured electrons to the corresponding par- 
tonic objects, in the range 0.8 < |r;| < 1.37, for different energies. 
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8. Unfolding 



Any measured observable is influenced by imperfections of the used measurement apparatus and 
procedure itself, such as limited resolution of the detector response, the detector acceptance and 
possible object and event selections which are applied to the data. Due to these distortions, 
any measurement of such observable does not fully represent the original (or true) quantity. 
Mathematically, the actual measurement can be considered to be a convolution of the true 
quantity with a function representing the overall detector and selection acceptance and the 
detector response. 

Let the true quantity be represented by a vector t (with entries tj and j = 1,2, ...,n() de- 
scribing the bin contents of a histogram, and the measured or reconstructed distribution by a 
corresponding vector k (with entries ki and i = 1,2, ...,nfc), respectively. The underlying detec- 
tor resolution effects can be described by a transition or response matrix R^°^, which contains 
the individual transition probabilities and hence the migrations between the observed elements 
of the distribution and the corresponding true values. 

Furthermore, the detector acceptance and applied selection can be quantified by an additional 
weight factor for each element of R^^, taking into account the probability of events from a 
particular entry of t being observed at all in the measurement process. This additional correction, 
which can be described by a second matrix R^'^^, together with the response matrix describing 
the resolution effects, yields the overall response matrix R: 

R = (8.1) 

denoting the transition probabilities between the observed distribution and the true distribution: 

k = Rt, (8.2) 

where 

Rij = P(observed in bin i \ expected in bin j) (8.3) 
= P{h\tj). (8.4) 

The concept of unfolding is illustrated in Figure [STU where an example distribution for an 
arbitrary variable x is shown at the different stages in the measurement process. 

In order to find an estimator for the true distribution given the measured distribution, an 
unfolding method |191j can be applied to correct for the respective acceptance and resolution 
effects. In this process, the response matrix is derived from an arbitrary reference sample, typi- 
cally using Monte Carlo simulations. This procedure is called the training step of the unfolding. 
The obtained response matrix has to be inverted in order to allow the unfolding of any mea- 
sured distribution to its corresponding true distribution. Since the response matrix represents 
the full resolution and acceptance information of the underlying measurement, the unfolding 
procedure can be performed model-independently, assuming that the detector simulation used 
in the training sample is sufficiently accurate. 

However, in most situations where unfolding is applied, an exact and unique inverse response 
matrix R~^ does not necessarily have to exist, such that 

RR~^ = I, (8.5) 
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8. Unfolding 



unfolding 




acceptance resolution 

Figure 8.1.: Schematic overview of the measurement and unfolding process. The true distribu- 
tion of an arbitrary variable x (left) is affected by acceptance effects (centre) and 
resolution effects (right) in the measurement process. The unfolding procedure 
attempts to reverse these effects to obtain the most probable true distribution 
corresponding to the given measured distribution. 



where / is the unity matrix. Hence, approximations are needed to perform the above matrix 
inversion to acceptable accuracy. 

Limited statistics in the reference sample and the resulting statistical fluctuations can lead to 
additional and inadvertent bin migration effects in the response matrix, which do not represent 
the underlying resolution and acceptance effects. Consequently, these contributions have to 
be suppressed in the matrix inversion process, achieved by applying a regularisation procedure 
in order to limit the propagation of statistical fluctuations into the unfolded distribution or 
quantity. This regularisation typically involves a cut-off or weight parameter representing the 
sensitivity of the unfolding approach with respect to short-ranged bin-by-bin changes. Hence, 
the regularisation can be regarded as a constraint on the smoothness of the response matrix and 
hence the unfolded distribution. 

The obtained approximate inverse matrix is applied to the distribution measured from 
data, m, and the respective unfolded distribution, -u, is obtained as estimator for the true 
distribution based on the given measurement: 

u = R-^m. (8.6) 

Several procedures have been developed to perform the inversion of the response matrix and 
the necessary regularisation. These unfolding methods will be briefly explained and evaluated 
with respect to their value and applicability for the measurement of the charge asymmetry using 
the observable Ac in the following. 

• Bin-by-bin correction - Neglecting bin migrations, a simple reweighting can be per- 
formed by defining a correction factor for each bin of the distribution measured on data 
based on the true and measured reference distributions |192j : 




Consequently, the unfolding is performed by application of the respective correction factors 
to the corresponding bins of the distribution rh obtained from the data measurement to 
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obtain the unfolded distribution u: 

Ui = Ci-mi. (8.8) 

Despite the simphcity of this method, it is prone to biases due to the reference Monte Carlo 
used in the training step, in particular the shape of the respective distributions. Since a 
priori no bin migrations are taken into account, it relies strongly on the correct description 
of the underlying physics and hence does not provide a model-independent approach. 
Furthermore, ambiguities in the determination of the statistical uncertainties may arise 
for cases where ki > tj, where the obtained relative uncertainty can be underestimated, 
being smaller than the expected uncertainty for an ideal detector (i.e. where nii would 
itself be an estimator of ti). 

• Singular value decomposition - This method employs a singular value decomposition 
(SVD) approach |193j in order to perform the inversion of the response matrix in the 
unfolding procedure. Singular value decomposition of a given real mxn matrix R involves 
a factorisation such that 

R = USV^, (8.9) 

where U and V denote m x m and n x n orthogonal matrices, respectively, such that 
UU'^ = U^U = Im and VV^ = V^V = In (with Im and I„ being the corresponding 
mxm and nxn unity matrices). Furthermore, S denotes an m x n diagonal matrix with 
non-negative diagonal elements, i.e. Sij = for i ^ j and Sa = Sj > 0. The matrix entries 
Si are called singular values of the matrix R. 

This form can be used in order to decompose the given measured distribution k and the 
unknown true distribution t into a series of orthogonal and normalised functions of the 
respective m and n classes by performing an appropriate rotation of the respective vectors. 
Consequently, the given initial system of linear equations as shown in Equation 18.21 is 
reduced to a diagonal system of equations, such that 

U^k = SV^t ^ d = Sz, (8.10) 

where d and z denote the rotated measured and true vectors, respectively. This procedure 
is in particular effective if the respective singular values Si are small and the statistical 
uncertainties on the entries of the measured distribution are large, since in such a case 
any exact inversion solution is dominated by statistical fluctuations and thus physically 
meaningless. 

By transformation of the given system of linear equations into the form of a weighted least 
squares minimisation, taking into account measurement uncertainties, a regularisation of 
the unfolding procedure can be achieved by the addition of a corresponding regularisation 
or stabilisation term to the expression to be minimised |194H196] . This introduces prior 
knowledge of the given problem and involves the requirement of minimal curvature of the 
obtained unfolded solution (i.e. the smoothness of the resulting distribution), eliminating 
statistical bin-to-bin fluctuations similar to the suppression of high-frequency harmonics in 
Fourier analysis. A regularisation parameter r defines the relative weight of the additional 
regularisation term with respect to the terms originating from the given system of linear 
equation in the minimisation: 

(^Rt- k^ ^ (^Rt- Pj+T {Ctf Ct, (8.11) 
where C denotes a matrix representing the prior condition on the solution. 
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8. Unfolding 



Bayesian iterative unfolding - A procedure based on an iterative approach to perform 
the inversion of the response matrix foUowing Bayes' theorem is appHed |197| . This ap- 
proach allows incorporating new knowledge to update a prior probability of observation of 
a given event |198H200] in an iterative procedure. 

In order to obtain the inverted response matrix, the posterior probability of obtaining the 
true distribution t given the measured distribution k is calculated accordingly, assuming 
prior knowledge Poitj) for the individual components of t based on the true distribution 
obtained in the Monte Carlo training step of the unfolding: 

Note that in this context, the probability P{ki\tj) is identical to the transition probability 
Rij contained in the response matrix. 

The obtained posterior probability distribution function is used as a prior in the next 
iteration, consecutively updating the existing knowledge about the respective probabilities 
with increasing number of iterations: 

Pi(tj) oc ^P(ij|A;,)-A:, oc^P(A;i|tj)-Po(ti)-A:i (8.13) 

i i 

P2{tj) a J^^(^i|ij)-A(ti)-A:, (8.14) 

i i 

Psitj) « ^P(tj|fci)-^"*«l^^(fci|ij)-^'2(ti)-A:^ (8.15) 



Regularisation of the Bayes iterative unfolding procedure can be achieved naturally by 
restricting the number of iterations Nn such that the underlying true distribution is re- 
covered within the statistical uncertainties, and bin-to-bin fluctuations which are of purely 
statistical nature are suppressed. For a large number of iterations, a convergent state is 
reached, yielding the true, but strongly fluctuating inverse of the response matrix (thus 
minimising any remaining systematic bias of the unfolded distribution at the cost of larger 
statistical uncertainties). The number of iterations necessary to reach convergence depends 
on different conditions, including the choice of binning, the strength of bin migrations in 
the response matrix (i.e. the magnitude of its off-diagonal elements), and the choice of 
prior. The optimal choice of A^it is case dependent and must be determined following a 
well-deflned procedure, balancing remaining bias and statistical uncertainty of the obtained 
result. 

Due to the fact that the bin-by-bin unfolding cannot account for bin migrations and does 
not allow for a regularised unfolding procedure, and hence is expected to be heavily model- 
dependent, this procedure was not eligible for the usage in this analysis. It only allows for a 
comparison of the measured asymmetry with the Standard Model prediction (unless the un- 
folding is performed for different model hypotheses), which can be achieved using the measured 
asymmetry directly as well. 

The SVD unfolding approach, despite providing a well deflned methodology, cannot perform 
an unfolding in more than one parameter in its current technical implementations, which would 
limit the analysis to an inclusive unfolding in a single parameter. 
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Given the requirements of the measurement, a Bayesian iterative unfolding was performed 
in this analysis in order to recover the inclusive \yt\ — \yt\ distribution and the resulting charge 
asymmetry observable Aq at the parton level. In particular, this approach allowed for a simulta- 
neous unfolding in multiple observables due to the fact that Bayesian unfolding is independent of 
the ordering of the classes/bins. Since many of the BSM models summarised in Chapter [2] pre- 
dict different magnitudes of the charge asymmetry for low and high M^j regions, a simultaneous 
unfolding of \yt\ — \yt\ and the invariant tt mass, Mjf, has been performed, taking into account 
bin migrations in both dimensions. Semileptonic tt events generated with the mc@nlo gener- 
ator have been used as reference sample to obtain the response matrix based on the detector 
simulation. Furthermore, a simple extraction of the covariance matrix of the unfolded distribu- 
tions |201| was possible. 

The SVD unfolding was performed for the inclusive measurement as a cross-check to verify 
the stability and consistency of the Bayes iterative approach (c.f. Appendix [C|) . 

The techniques utilised were available in the RooUnfold package |2U1| . which provided simple 
interfaces and efficient implementations for all three mentioned unfolding methods. 
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9. Systematic Uncertainties 



In addition to the statistical uncertainty originating from limited data statistics, there were a 
multitude of systematic effects that can have an impact on the performed measurement. These 
effects were studied individually and a corresponding systematic uncertainty on the measurement 
result was assigned for each contribution. 

In order to evaluate the individual effects, the analysis was performed for each systematic 
under consideration with a modified response matrix and/or background contribution depending 
on the modelled effect. The changes typically corresponded to an uncertainty of one or more 
parameters (e.g. a shift of the muon trigger efficiencies according to the respective uncertainties) 
or an alternate model (e.g. a different MC generator used for the simulation of ti events). The 
uncertainty was extracted in each case based on the shift in the measurement central value when 
changing the parameters accordingly, and was symmetrised. 

In order to suppress the statistical component of the obtained uncertainty inherent in the re- 
evaluation of the central value by changing different parameters of the performed measurement, 
the requirement of at least one b tagged jet in the event selection was replaced by a reweight- 
ing method. The \yt\ — \yt\ distribution for the Monte Carlo based background contributions 
(VF+jets, Z+jets, single top and diboson background) was obtained by direct application of the b 
tag weights to the events passing the nominal event selection without the requirement of at least 
one b tagged jet. This approach is very similar to the trigger reweighting approach described in 
Chapter 14.2.31 The same procedure was applied for the simultaneous unfolding in |yf| — and 
Mjf. Furthermore, the normalisation of the resulting distributions was modified to match the 
distribution obtained using the nominal event selection to avoid the introduction of a potential 
bias in the background subtraction. Control plots showing the agreement of the two approaches 
can be found in Appendix [El Since the weighted distributions were compatible within statistical 
uncertainties with the nominal distributions, no additional systematic uncertainty was assigned. 
The bin- by-bin statistical uncertainty in the background distribution is reduced by up to 25 %. 

The following systematics were considered and evaluated for this analysis. All contributions 
were assigned to the Monte Carlo prediction (signal and/or background, were applicable). 

• QCD multijet uncertainty - The QCD multijet background was estimated with data 
driven methods in both the electron and muon channel. Since both the shape and the 
overall normalisation of the estimation can only be verified to a limited extend, a very 
conservative systematic uncertainty of 100 % was quoted despite the availability of more ad- 
vanced estimates (c.f. Chapter 16.3. ip . following the recommendations of the performance 
groups. Results following a less conservative approach are discussed in Appendix [Dj The 
QCD multijet background normalisation was shifted up and down by 100% to quantify 
this uncertainty. Since the QCD multijet background contribution is intrinsically charge 
symmetric, but enters the background subtraction, the normalisation and not shape was 
the dominant source of systematic uncertainty. 

• Jet energy scale - The jet energy scale (JES) uncertainty was derived using information 
from test beam data, LHC collision data and simulation and was taken into account by 
scaling up and down the energy of all considered jets by 1 cr of the associated transverse 
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9. Systematic Uncertainties 



momentum uncertainty based on the jet pT and r/ |202j . The full event selection and 
kinematic reconstruction has been re-run with the scaled jets. In addition, the missing 
transverse energy has been re-evaluated, taking into account the scaled contributions of 
the jets in and py. For jets within the acceptance range, the JES uncertainty varied 
from about 2.5 % for high jets in the central detector region to about 14 % for low 
Pt jets in the forward region. 

• Pile-up (JES) - Depending on the instantaneous luminosity and the amount of vertices, 
an event weight was assigned to match the pile-up contribution in the simulation to data. 
An additional systematic uncertainty of 5% (7%) for low jets or 2% (3%) for high 
Pt jets in the \r]\ < 2.1 (2.1 < \7]\ < 4.5) region was added to the JES uncertainty in 
quadrature |203j . 

• b jet energy scale - In order to account for the difference of the energy scale for b 
jets with respect to light quark jets, all b jets (i.e. jets which have a matched truth b 
quark in simulations) were scaled by an additional fraction ranging from 2.5 % in the low 
Pt jet region to 0.76% in the high jet region and added to the JES uncertainty in 
quadrature |202j . 

• Jet reconstruction efficiency - The jet reconstruction efficiency (JRE) was evaluated 
by randomly dropping jets from events with a probability of about 2% |167j . The resulting 
difference with respect to the nominal case was symmetrised and quoted as JRE systematic 
uncertainty. 

• Jet energy resolution - A smearing of the jet transverse momentum corresponding to 
a resolution of about 10 % was applied to Monte Carlo events as systematic to reflect the 
difference between the jet energy resolution (JER) observed on data and Monte Carlo |167j . 
The resulting discrepancies were symmetrised and quoted as JER systematic uncertainty. 

• Muon efficiencies - In order to account for the trigger and reconstruction efficiencies 
for muons [20411205] . global and object based scale factors and efficiencies were taken into 
account and a systematic uncertainty was assigned on an event-by-event basis (for global 
scale factors) or on an object basis. These were combined into an overall muon efficiency 
uncertainty of the order of 1%. In addition, a one-sided uncertainty of 1.5% (events 
with or 1 6 tagged jet) or 2.2 % (events with more than 1 b tagged jet) was assigned to 
the muon trigger efficiency to account for a mis-modelled and as a consequence discarded 
trigger object matching algorithm in the Monte Carlo samples. 

• Muon scales and resolution - Since the Monte Carlo muon momentum scales and 
resolution differed from the ones observed in data, the muon momentum was smeared and 
a scaling of up to 1.5 % was applied on object level to account for this discrepancy |171j . A 
systematic uncertainty at the sub-percent level was assigned by scaling up and down both 
the momentum scaling and smearing by 1 o" according to the respective uncertainty. In 
addition, the missing energy was re-evaluated with the modified four-vector information. 
The full event selection and kinematic reconstruction was performed for the different scales, 
resulting in a symmetrised systematic uncertainty based on the comparison of the different 
results of the measurement. 

• Electron efficiencies - In order to account for the trigger and reconstruction efficiencies 
for electrons, global and object based scale factors were taken into account |172j . A 
systematic uncertainty was assigned on an event basis (for global scale factors) or on an 
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object basis, which were combined into an overall electron efficiency uncertainty of the 
order of 1 %. 

• Electron scales and resolution - In order to take into account discrepancies between the 
electron energy resolution on Monte Carlo and data, a Gaussian smearing procedure was 
applied to the electron energy for Monte Carlo events to reflect the resolution in data jl72j . 
In addition, the electron energy in data was corrected to account for a scaling mis-match 
between data and Monte Carlo. Both systematic uncertainties were of the order of 1 % to 
2 % and were assigned to the Monte Carlo prediction for consistency. 

• b tag scale factors - Due to discrepancies in the b tagging efficiencies and fake rates 
between data and simulation, all Monte Carlo jets were assigned a specific weight to 
account for this effect [16811206] . The obtained b tag weights for each jet were combined 
into an event weight by multiplication (corresponding to a logical AND of all jets taken 
into account). The provided scale factors contained uncertainties which result in small 
shape variations. In order to determine the deviation in the shapes from the nominal 
case due to the b tag scaling and to quantify the corresponding systematic uncertainty 
on the measurement, the resulting samples were shifted up and down by the provided 
uncertainties. These were of the order of 8%, depending on the jet pT and rj. 

• PDF uncertainty - For the mc@nlo signal Monte Carlo, CTEQ6.6 PDFs were utilised 
to model the incoming partons to the hard scattering process, as described in Chapter 16.21 
The impact of the choice of PDFs was evaluated by varying the eigenvalues of the CTEQ 
parametrisation |2U7| or by comparison with the respective MRST2001 parametrisation 
[208] . using the Lhapdf tool |209j . Event weights were determined and the variations in 
the resulting pseudosamples were taken into account as PDF systematic uncertainty. 

• LAr defects - Parts of the LAr calorimeter readout electronics were inoperative during 
a significant time period of data taking due to a technical problem. Having occurred 
after the production of the used Monte Carlo samples, it was necessary to correct for 
the resulting mismatch between data and simulated events at the analysis level. Monte 
Carlo events were dropped with a probability corresponding to the relative amount of data 
affected (84.0 %) if an electron or a jet entered the region of degraded acceptance (taking 
into account the isolation requirements). A systematic uncertainty was assigned based 
on different transverse momentum requirements for the jets taken into account for this 
procedure after symmetrisation. 

• ISR and FSR - In order to take into account initial and final state radiation, which 
can introduce additional jets in the observed events, different Monte Carlo samples with 
varying ISR and FSR contributions |187j (ISR and FSR contributions scaled up and down 
independently and in combination) were evaluated by replacing the nominal signal sample 
in the measurement. The systematic uncertainty was quoted as the maximum relative 
deviation from the nominal leading order sample observed in these variations and applied 
to the MC@NLO prediction. The parameters were varied in a range comparable to those 
used in the Perugia Soft /Hard tune variations |210j . 

• ti modelling - The impact of using different MC generators for the signal process mod- 
elling was studied. In addition to the mc@nlo generator, the powheg generator was 
used for comparison and a symmetrised systematic uncertainty was assigned based on the 
variations in the measurement results for the alternate modelling. 
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9. Systematic Uncertainties 



• Parton shower / fragmentation - In addition to the matrix element level MC gen- 
erator, the effect of different showering models was taken into account by comparing the 
results for the powheg generator with showering performed by PYTHIA and by herwig, 
and a symmetrised systematic uncertainty was assigned based on the variations in the 
measurement results for the alternate shower modelling. 

• Top mass - Since the top mass parameter was considered fixed, the uncertainty on the 
measurement of the mass was taken into account separately. Two different Monte Carlo 
samples generated with different mass parameters (scaled up and down to 180 GeV and 
170 GeV, respectively) were used and the observed deviations were linearly interpolated 
according to the actual uncertainty of 0.5 % of the top mass measurement [13] and a 
symmetrised systematic uncertainty was quoted. 

• VF+jets background uncertainty - The W+jets background normalisation was es- 
timated with a data driven method and a systematic uncertainty based on the limited 
statistics available and several systematic contributions to the method have been evalu- 
ated (see Chapter 16.41 for details). An overall T^+jets normalisation uncertainty of 22.4% 
and 22.7% was obtained in the muon+jets and electron+jets channel, respectively. In 
addition, the Vl^+jets background shape uncertainty has been evaluated by modifying sev- 
eral generator parameters such as the renormalisation scale or the functional form of the 
factorisation scale compared to the nominal alpgen parameters, based on the leading jet 
Pt- a symmetrised systematic uncertainty on the charge asymmetry measurement has 
been assigned based on the deviations for two different sets of parameters with respect to 
the nominal results. 

• .Z-|-jets background uncertainty - In order to quantify the uncertainty on the Z+jets 
contribution normalisation, a Berends-Giele scaling uncertainty j211j was calculated, cor- 
responding to an overall normalisation uncertainty of 34%. In addition, the .Z^+jets back- 
ground was determined independently from both the alpgen and sherpa generator to 
quantify the shape uncertainty, which was quoted based on the symmetrised discrepancy of 
the results obtained with sherpa with respect to the nominal case, in which alpgen sam- 
ples were used. 

• Charge mis-identification - Since the detector momentum resolution is finite, and the 
lepton charge was identified by taking into account the bending radius of the particle 
track, a certain probability for mis-identifying the lepton charge exists, especially for high 
transverse momentum leptons due to their almost straight trajectories. This probability 
was evaluated on Monte Carlo and on data to be of the order of 0.2 % to 0.5 % in the central 
detector region, and up to 2.5 % in the forward region, respectively, in the electron+jets 
channel [183j . In the muon+jets channel, the probability was found to be below 0.003% 
in all cases. A corresponding symmetrised systematic uncertainty on the measurement of 
the charge asymmetry was determined. 

• b tag charge - A dependency of the b tag efficiencies on the b quark charge can lead to 
a bias in the measurement due to the requirement of at least one b tagged jet. Hence, 
a simple Monte Carlo study on parton level was performed by simulating a difference in 
the b tag efficiency of 5 % between b and b quarks. The resulting impact on the charge 
asymmetry on parton level was studied and the difference with respect to the nominal 
case (assuming identical tagging efficiencies for b and b quarks) was quoted as systematic 
uncertainty. 
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• MC generator statistics - Since the signal Monte Carlo sample entered directly into the 
unfolding response matrix and statistical fluctuations in the bins of this matrix can have 
an impact on the unfolding process, an additional ensemble test was performed by fluc- 
tuating the obtained nominal response matrix on a bin-by-bin basis following a Gaussian 
probability distribution (since the statistics from the mc@nlo sample were very high a 
Gaussian model could safely be assumed). Uncertainties of the order of 0.3% to 3% have 
been obtained, depending on the statistics in each bin of the response matrix. 

• Unfolding convergence - Based on the convergence evaluation which was used to de- 
termine the optimal amount of regularisation in the unfolding process (for details refer 
to Chapter 110. 3p . a remaining absolute uncertainty of l%o, corresponding directly to the 
choice of convergence criterion, was assigned, representing the potential remaining change 
with respect to further increase in regularisation. 

• Unfolding bias - Closure tests have been performed using ensembles of pseudodata to 
quantify any remaining bias from the unfolding at the chosen regularisation strengths 
by obtaining pull distributions for the measured asymmetry, normalised with respect to 
the respective unfolding statistical uncertainty. The corresponding distributions can be 
found in Figure lA. 81 in Appendix [Al The remaining bias was extracted from the residuals 
in the pull distributions in the closure tests and was taken into account as an additional 
systematic uncertainty. The residual bias after unfolding was extracted from the respective 
pull distributions corresponding to the regularisation used in the individual cases and taken 
into account as an additional relative uncertainty on the unfolded results, which was of 
the order of 1% to 11%. 

• Other backgrounds - For the small backgrounds from single top and diboson production, 
only normalisation uncertainties were considered. For the single top contribution, the 
uncertainties of the approximate next-to-next-to-leading order prediction (for details, refer 
to Chapter 12.2. 2p were taken into account (corresponding to an uncertainty of about 11 %), 
while for the diboson production, an overall uncertainty of 5 % was assumed. 

• Luminosity - The relative uncertainty on the measurement of the integrated luminosity 
of the used data sample was estimated to be 3.7% |212j and was taken into account for 
the measurement. 

• Pile-up - In order to take into account the difference in pile-up conditions between Monte 
Carlo and data, an event based reweighting was performed to reflect the distribution of 
average observed number of interactions per bunch crossing on data. The impact of the 
pile-up conditions on the measurement before unfolding is shown in detail in Chapter 110. 3[ 
Since no significant pile-up dependency was observed, no additional systematic uncertainty 
was assigned. 
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10. Results 



10.1. Event Yields & Control Plots 

The final expected and observed number of events in both the muon+jets and the electron+jets 
channel after performing the event selection can be found in Table IIO.H both without and with 
the requirement of at least one b tagged jet. 



Channel 


muon- 


hjets 


electron+jets 


Selection 


> 6 tags 


> 1 6 tags 


> b tags 


> 1 b tags 


tt 

Single top 


7187 ± 588 
458 ± 40 


6342 ± 519 
366 ± 32 


4832 ± 395 
320 ± 28 


4256 ± 348 
256 ± 22 


W^+jets (data) 

Z+jets 

Diboson 


8617 ± 1173 
940 ± 333 
134 ± 7 


1390 ± 311 
134 ± 47 
22 ± 2 


5407 ± 768 
756 ± 267 
80 ± 5 


883 ± 200 
105 ± 37 
13 ± 1 


QCD multijct (data) 


1515 ± 757 


519 ± 519 


944 ± 472 


247 ± 247 


Total background 


11663 ± 1436 


2431 ± 608 


7508 ± 941 


1504 ± 321 


Total expected 


18850 ± 1552 


8774 ± 799 


12340 ± 1020 


5760 ± 473 


Observed 


19639 


9124 


12096 


5829 



Table 10.1.: Observed number of data events in comparison to the expected number of Monte 
Carlo signal events and different background contributions for the event selection, 
both with and without the requirement of at least one h tagged jet. The QCD 
multijct and W^+jets contributions were estimated using data driven methods 
(c.f. Chapter in]). Uncertainties are statistical and include the respective sys- 
tematic uncertainties on the normalisation for QCD multijct and Ty+jcts and 
the cross-section uncertainties on all other contributions. For the QCD multijct 
background, a conservative 50 % (100 %) overall normalisation uncertainty for the 
selection without (with) requiring at least one b tagged jet was assumed. 

Control plots for the full event selection as described in Chapter [5l showing a comparison 
of the observation and expectation for several object and event quantities, can be found in 
Figure [10. II and Figure [1021 for both the muon-l-jets and the electron-|-jets channel, respectively. 
The uncertainties on the expectation include statistical and the leading systematic uncertainties, 
that is QCD multijet background and VF-|-jets normalisation, luminosity, jet energy scale, h tag 
scale factors and ti cross-section uncertainty. 

Due to the tighter definition of selection criteria in the electron -|-jets channel, the number of 
events in the electron-l-jets channel was significantly lower than it was the case for the muon-l-jets 
channel. This was necessary in order to reduce the contribution from the increased number of 
electron fakes originating from the QCD multijet background with respect to the muon channel, 
where the expected fake rate was significantly lower. The overall agreement between Monte 
Carlo prediction and data was very good in both channels. Additional control plots for the 
same quantities without the explicit requirement of at least one h tagged jet can be found in 
Figure [X?3l and Figure of Appendix [XI for completeness. 
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Figure 10.1.: Control plots for the muon+jets channel. From the top left to the bottom right, 
the transverse momentum px, the pscudorapidity 77 and the azimuthal angle 4> of 
the selected muon arc shown. Additional plots show the transverse momentum 
of the leading jet (ptOi)), the W transverse mass mT(W^) and the transverse 
missing energy _^t- Uncertainties are statistical and for T4^+jets also include 
systematic uncertainties on normalisation. For the QCD multijet background, a 
conservative 100 % systematic uncertainty was assumed. In addition, the uncer- 
tainties on luminosity, jet energy scale, h tag scale factors and tt cross-section 
are shown. 
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10.1. Event Yields k Control Plots 





(W) [GeV] [GeV] 



Figure 10.2.: Control plots for the electron+jets channel. From the top left to the bottom 
right, the transverse momentum px, the pscudorapidity ry and the azimuthal 
angle (j) of the selected electron are shown. Additional plots show the transverse 
momentum of the leading jet (pt(ji)), the W transverse mass m'T[W) and the 
transverse missing energy _^t- Uncertainties are statistical and for ly+jets also 
include systematic uncertainties on normalisation. For the QCD multijet back- 
ground, a conservative 100 % systematic uncertainty was assumed. In addition, 
the uncertainties on luminosity, jet energy scale, b tag scale factors and tt cross- 
section are shown. 
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Figure 10.3.: Control plots for the tt event reconstruction, on the left for the muon+jets, 
on the right for the electron+jets channel. The top row shows the logarithmic 
likelihood logL of the kinematic fit, followed by the invariant mass Mfi and the 
transverse momentum of the ti system. Uncertainties arc statistical and for 
VF+jets also include systematic uncertainties on normalisation. For the QCD 
multijet background, a conservative 100 % systematic uncertainty was assumed. 
In addition, the uncertainties on luminosity, jet energy scale, b tag scale factors 
and ti cross-section arc shown. 
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10.2. Measurement of the Charge Asymmetry 



In addition to the control plots showing basic object and event kinematics, several more 
complex quantities based on the kinematic event reconstruction, as described in Chapter [71 
can be found in Figure [IU.3[ The respective distributions for the logarithmic likelihood of the 
kinematic fit, and the invariant mass and transverse momentum of the reconstructed system 
in both the muon+jets and the electron+jets channel are shown. The agreement between data 
and prediction is very good, indicating a proper modelling of the ti signal and background 
kinematics in the various Monte Carlo samples and data driven background estimates. 

10.2. Measurement of the Charge Asymmetry 

Events passing the described event selection were taken into account to determine the differential 
and integrated charge asymmetry based on the distribution of — A subtraction of the 
predicted background contributions was performed on the distributions measured on data to 
obtain an estimate for the tt signal contribution only. 

The corresponding distributions are shown in Figure 110.41 for both the inclusive measure- 
ment of the \yt\ — \yt\ distribution and the corresponding measurement for the two cases of 
Mjf < 450 GeV and M^j > 450 GeV. For the M^f dependent measurement, an additional re- 
quirement on the event reconstruction logarithmic likelihood of logL > —52 (see Figure [IU.3|) 
has been applied in order to improve the resolution in the tt invariant mass. The respective 
relative resolution in M^j for both channels before and after the application of the additional 
requirement on logL can be found in Figure lA. II and Figure [A. 21 in Appendix lAl The relative 
resolution improves from 28.8% to 18.4% in the muon+jets channel and from 28.8% to 17.7% 
in the electron+jets channel. 

The resulting charge asymmetries for the observable A^^^° obtained after background subtrac- 
tion can be found in Table 110.21 for both the muon+jets channel and the electron+jets channel. 
The observed results are shown alongside the predicted tt charge asymmetries obtained with 
MC@NLO for comparison. The uncertainties on the prediction correspond to the limited Monte 
Carlo statistics in the used sample after the applied selection. As can be seen in the table. 



Measurement 


Observed 


Predicted (mc@nlo) 


Muon+jets Channel 


Inclusive 

Mu < 450 GeV 

Mu > 450 GeV 


-0.010 ± 0.015 (stat.) ± 0.008 (syst.) 
-0.012 ± 0.024 (stat.) ± 0.007 (syst.) 
-0.011 ± 0.022 (stat.) ± 0.007 (syst.) 


0.0021 ± 0.0011 (stat.) 
-0.0043 ± 0.0018 (stat.) 
0.0064 ± 0.0015 (stat.) 


Electron+jets Channel 


Inclusive 

Mff < 450GeV 

Mii-> 450GeV 


-0.034 ± 0.019 (stat.) ± 0.010 (syst.) 
-0.060 ± 0.033 (stat.) ± 0.019 (syst.) 
-0.032 ± 0.028 (stat.) ± 0.015 (syst.) 


-0.0003 ± 0.0014 (stat.) 

0.0007 ± 0.0022 (stat.) 
-0.0010 ±0.0017 (stat.) 



Table 1 0. 2. : Measured values of the charge asymmetry observable A^^'^° for the muon+jets and 
clcctron+jets channel after subtraction of the various background contributions. 
The results for the inclusive measurement and the respective measurements for 
Mtt < 450 GeV and Mu > 450 GeV are shown together with the mc@nlo pre- 
dictions. 



the obtained results were compatible with the mc@nlo Standard Model prediction in all cases 
within the statistical and systematic uncertainties. However, in the electron+jets channel, a 
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Distributions for \yt\ — on the left for the muon+jets, on the right for the 
electron+jets channel The top row shows the distribution for the inclusive mea- 
surement, while the lower rows show the corresponding distributions requiring 
Mil < 450 GeV and Mjf > 450 GeV, respectively. Uncertainties are statistical 
and for VF+jets also include systematic uncertainties on normalisation. For the 
QCD multijet background, a conservative 100 % systematic uncertainty was as- 
sumed. In addition, the uncertainties on luminosity, jet energy scale, b tag scale 
factors and tt cross-section arc shown. 
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10.3. Unfolding 



tendency towards more negative asymmetries was observed, in particular for the measurement 
of for < 450 GeV. A breakdown of the individual systematic uncertainties taken into 
account can be found in Table [BTT] and Table [B?2] of Appendix [B] for completeness. 

In addition, the integrated charge asymmetry after background subtraction as a function 
of the invariant tt mass, M^j is shown in Figure 110.51 together with the Standard Model predic- 
tion obtained from simulated tt events. The observed asymmetry values show no statistically 
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Figure 10.5.: Integrated charge asymmetry after background subtraction as a function 

of Mil- The observed asymmetries (red, dashed) and the asymmetries predicted 
by MC@NLO (blue) are shown. Uncertainties are statistical only. 

significant deviation from the Standard Model prediction. Since, in addition, no systematic 
uncertainties were included, the observed discrepancies are negligible. 



10.3. Unfolding 

In order to perform the Bayes iterative unfolding and to correct the measured differential and 
integrated charge asymmetry for detector resolution and acceptance effects, the corresponding 
response matrix has been obtained using the full set of ti signal Monte Carlo events generated 
with MC@NLO (15 000 000 simulated events). A unified binning for the true and reconstructed 
distributions (and consequently for the response matrices) in the variable |yt| — has been used, 
employing six bins with variable width to ensure sufficient statistics in the tails of the respective 
distributions. Bin edges at {—3.0, —1.2, —0.6, 0.0, 0.6, 1.2, 3.0} for the inclusive unfolding and 
{-3.00,-0.96,-0.48,0.0,0.48,0.96,3.00} for the simuhaneous unfolding in \yt\ - \yi\ and M^i 
were chosen. 

Furthermore, for the simultaneous unfolding in \yt \ — |yf| and M^j, an additional requirement 
on the event reconstruction logarithmic likelihood of logL > —52 (see Figure [10. 3p has been 
applied in order to improve the resolution of the ti invariant mass in the same way as for the 
measurement of the charge asymmetry before unfolding. 

The response matrix representation used in the inclusive measurement of the distribution 
of \yt\ — \yt\ and of is shown in Figure [10.6I for both the muon+jets and electron+jets 
channel. Similarly, the corresponding matrices for the simultaneous unfolding in \yt \ — \yt\ and 
Mjf can be found in Figure 110.71 The transition probability information of the corresponding 
two-dimensional distributions (with 2x6 bins) was encoded in the matrix by linearisation of 
the respective histograms, concatenating the bins from the two M^i regions into a single one 
dimensional histogram (with 1 x 12 bins). Hence, the associated response matrix contained 2x2 
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Figure 1 0. 6. : Unfolding response matrices for the inclusive unfolding of the charge asymmetry 
distribution |yt | — Ij/tl- The bin migration probability corresponds to the box sizes 
displayed in the response matrix, and is shown independently for the muon+jets 
channel (left) and the electron+jets channel (right). 



quadrants, corresponding to the transition probabilities for the true and reconstructed |yt| — 
distributions in the low and high M^j regions of phase space, respectively. 
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Figure 10.7. 



Unfolding response matrices for the simultaneous unfolding of the distribution 
in \yt \ — lutl and M^j showing the discrete bin correspondence of the truth distri- 
bution with respect to the distribution obtained after event selection and recon- 
struction. The bin migration probability corresponds to the box sizes displayed 
in the response matrix, and is shown independently for the muon-f-jets chan- 
nel (left) and the electron-f-jets channel (right). The four quadrants represent 
the respective transition probabilities for the true and reconstructed \yt \ — \yt\ 
distributions in the low and high Mf-i regions of phase space. 



A closure test has been performed to verify that the Bayesian iterative unfolding approach 
can be used to recover an arbitrary asymmetry present in the true distribution. The full 
MC@NLO tt signal sample was used as a basis to verify the correct unfolding response for different 
artificially injected truth level asymmetries. This was achieved by performing an event-by-event 
reweighting of the ti signal sample to asymmetries of -10 %, -5 %, %, 5 % and 10 % in the vari- 
able j4^"^ by systematically increasing the weights of events with — < and decreasing 
the weights of events with — > by the same amount. 

Pseudoexperiments, using statistically independent sets of pseudodata corresponding to the 
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10.3. Unfolding 



statistics expected in the used data sample (after background subtraction), were conducted 
based on a Poissonian fluctuation of the respective distribution in \yt \ — \ yt\, taking into account 
the additional event weights determined in the asymmetry reweighting. Ensemble tests were 
performed to confirm the linearity of the unfolding response in the true value of the chosen 
charge asymmetry observable Aj^"'^. Furthermore, the dependency of the obtained results on 
the regularisation of the unfolding procedure was studied. 

The average unfolded value of Aq^^ obtained from the sets of pseudodata as a function of the 
injected true value of ^^"^ can be found for both the inclusive unfolding and the simultaneous 
unfolding in \yt\ — \yt\ and M^j in Figure llU.8l for different regularisation strengths, using A^it = 5, 
10, 20, 40 and 80 iterations in the Bayesian unfolding. 

A straight line fit using a parametrisation of the form 

A^"^ = a • ^g^^ + 6 (10.1) 

has been performed, where a and b denote the slope and offset parameters, respectively. The 
procedure has been repeated for each of the individual choices of Nn to obtain the calibration 
curves in order to verify the linearity and determine the respective slopes and offsets of the 
fit. These provided a measure for the average remaining bias in the unfolded differential and 
integrated asymmetry after unfolding for different regularisation strengths. As can be seen in the 
fits, a slope close to one could be achieved in all cases, indicating a proper average correspondence 
of the unfolded asymmetry to the injected true value, independent of the strength of the injected 
asymmetry. 

In addition, the expected statistical uncertainty on after the unfolding procedure is 

shown in Figure lA. 51 in Appendix lAl as a function of A^it. As expected, the uncertainty increases 
with the number of iterations used in the regularisation due to the increasing sensitivity to 
bin-to-bin statistical fluctuations in the inversion process of the response matrix. The expected 
uncertainties are independent of the injected value of ^g™. 

In order to achieve a convergent state in the iterative unfolding process, the dependency of 
the unfolded result as a function of the number of iterations has been studied with respect to a 
defined convergence criterion. For individual ensembles of pseudodata, the unfolding procedure 
was considered to be converged if the absolute change in the unfolded asymmetry in terms of 
A^^, AA^^, for a given amount of iterations A^it = i with respect to the previous amount of 
iterations Njt = i — 1 was lower than l%o, i.e. if 

AA^'^ = |A^"^(iVit = i)- ^(iVit = i - 1)1 < 0.001. (10.2) 

Figure [A. 61 in Appendix IX\ shows the percentage of ensembles which have reached the defined 
convergence criterion as a function of regularisation strength, parametrised by the amount of 
iterations A'^it. At A^it = 40 (inclusive unfolding), and A^it = 80 (simultaneous unfolding in \yt\ — 
and -/V%) iterations, all ensembles of pseudodata were ensured to have reached convergence 
for all injected value of Aq^^. Hence, these regularisation strengths were chosen for the inclusive 
unfolding and the simultaneous unfolding in \yt\ — \yt\ and M^f, respectively. An additional 
systematic uncertainty of l%o, corresponding directly to the choice of convergence criterion was 
assigned to the unfolded value of A'^^ (c.f. Chapter [9]). 

In addition, a further cross-check was performed to verify that the unfolding procedure reached 
a convergent state by determining the average standard deviation of the variable A^™^ for the 
used set of ensembles as a function of A^it, which shows a monotonous falling behaviour in all 
cases as expected. The corresponding additional control plots can be found in Figure IA.7I in 
Appendix [XI for completeness. 
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Figure 10.8.: The obtained overall inclusive asymmetry after unfolding as a function of in- 
jected true asymmetry for different regularisation parameters A^it for 
both the muon+jets channel (left) and electron+jets channel (right). The 
top row shows the respective distribution for the inclusive measurement, while 
the lower rows show the corresponding distributions for Mfi < 450 GeV and 
> 450 GeV, respectively. For the simultaneous unfolding in \yt \ — \yt\ and 
Mtf, a cut on the event reconstruction likelihood log L was applied to improve 
the Mtt resolution of the selected events. 
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Despite the large amount of iterations necessary to reach convergence, in particular for the 
simultaneous unfolding in \yt\ — and M^j and the associated expected statistical uncertain- 
ties, this approach allowed performance of the unfolding following a well-defined and model- 
independent procedure. Furthermore, the remaining expected bias in the unfolding procedure 
was minimised, as can be seen in the calibration curves in Figure [IU.81 for the used choices 
of Njt = 40 and Njt = 80, respectively. The corresponding slope and offset parameters ex- 
tracted from the straight line fit in the calibration can be found in Table [TO^ for completeness. 
Additional closure tests have been performed using ensembles of pseudodata to quantify any 



Unfolding 


electron+jets channel 


muon-l-jets channel 


a (slope) 


b (offset) 


a (slope) 


b (offset) 


Inclusive 


1.009 ± 0.004 


-0.0032 ± 0.0002 


1.002 ± 0.003 


-0.0037 ± 0.0002 


< 450 GeV 


1.041 ± 0.001 


0.0030 ± 0.0005 


1.027 ± 0.008 


-0.0064 ± 0.0004 


> 450GeV 


0.984 ± 0.004 


-0.0063 ± 0.0003 


0.992 ± 0.004 


0.0007 ± 0.0003 



Table 1 0. 3. : Slopes and offsets from the linear fit in the unfolding calibration. The parameters 
were obtained for a linear fit of the average unfolded value of as a function 
of the true value, obtained from sets of pseudoexperiments for 40 (inclu- 

sive unfolding) and 80 (simultaneous unfolding in \yt \ — \yt\ and M^j) iterations, 
respectively. For the simultaneous unfolding in \yt\ — \yi\ and M^^, a cut on the 
event reconstruction likelihood log L was applied to improve the Mtt resolution 
of the selected events. The shown fit parameter uncertainties are statistical only. 

remaining bias from the unfolding at the chosen regularisation strengths and a corresponding 
systematic uncertainty was assigned to the unfolded result accordingly (c.f. Chapter [9]), while 
no correction of the unfolded results for the obtained calibration was performed. 

The unfolded integrated asymmetries in A'^^ are shown in Table 110.41 for both the inclusive 
unfolding and for the simultaneous unfolding in \yt \ — \yt\ and M^j, using 40 and 80 iterations in 
the unfolding process, respectively. In addition. Figure [10.91 shows the obtained distributions 



Unfolding 


Observed 


Predicted (mc@nlo) 


Muon+jets Channel 


Inclusive 

Mu < 450 GeV 

Mtt > 450 GeV 


-0.002 ± 0.036 (stat.) ± 0.023 (syst.) 
-0.002 ± 0.084 (stat.) ± 0.049 (syst.) 
-0.003 ± 0.045 (stat.) ± 0.034 (syst.) 


0.0056 ± 0.0003 (stat.) 
0.0024 ± 0.0004 (stat.) 
0.0086 ± 0.0004 (stat.) 


Electron+jets Channel 


Inclusive 

Mtt < 450 GeV 

Mtt > 450 GeV 


-0.047 ± 0.045 (stat.) ± 0.028 (syst.) 
-0.196 ± 0.119 (stat.) ± 0.091 (syst.) 
-0.016 ± 0.055 (stat.) ± 0.035 (syst.) 


0.0056 ± 0.0003 (stat.) 
0.0024 ± 0.0004 (stat.) 
0.0086 ± 0.0004 (stat.) 



Table 10. 4-: Unfolded values of the charge asymmetry observable for the muon+jets 

and electron+jets channel. The results for the inclusive measurement and the 
respective results for the simultaneous unfolding in \yt \ — \yi\ and M^f for M^t < 
450 GeV and M^t > 450 GeV are shown. For the simultaneous unfolding in — 
\yt\ and M^t, a cut on the event reconstruction likelihood logL was applied to 
improve the M^j resolution of the selected events. Furthermore, the respective 
MC@NLO predictions are shown. 
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Figure 10.9.: Unfolded distribution of \yt \ — \yt\: normalised to unity for both the muon+jcts 
channel (left) and clcctron+jets channel (right). The top row shows the distribu- 
tions for the inclusive measurement, while the lower rows show the corresponding 
distributions for M^i < 450 GeV and Mfi > 450 GeV, respectively. For the si- 
multaneous unfolding in \yt \ — \yt\ and Mft, a cut on the event reconstruction 
likelihood log L was applied to improve the M^f resolution of the selected events. 
The uncertainties include both statistical and systematic shape components. 
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after unfolding for both the inclusive unfolding and the simultaneous unfolding in \yt\ — |yf| 
and Mjf. The unfolded distributions have been normalised to unity and the shape uncertain- 
ties obtained from all systematic effects described in Chapter and the bin-by-bin statistical 
uncertainties have been included. 

Furthermore, the covariance matrices corresponding to the unfolded \yt\ — \yt\ distributions 
are shown in Figure 110.101 and 110.111 for the inclusive unfolding and the simultaneous unfolding 
in \yt\ — \yt\ and M^f, respectively. 
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Figure 10.10.: Covariance matrices corresponding to the unfolded — |j/t| distribution for the 
muon+jets channel (left) and the electron+jets channel (right) for the inclusive 
unfolding. The numbers inside the boxes represent the values and the sign of 
the correlation among the different bins and have been scaled by a factor of 
1000 for readability. 
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Figure 10.11.: Covariance matrices corresponding to the unfolded — \yt\ distribution for 
the muon-|-jets channel (left) and the electron-t-jets channel (right) for the 
simultaneous unfolding in \yt \ — \yt\ and Mti- The numbers inside the boxes 
represent the values and the sign of the correlation among the different bins 
and have been scaled by a factor of 1000 for readability. 



A summarised list of all systematics and their contribution to the overall systematic uncertain- 
ties can be found in Table [T03] for the inclusive unfolding and in Table [T0T6l for the simultaneous 
unfolding in \yt\ — \yt\ and M^j. 

All systematic uncertainties based on a replacement of the unfolding matrix which were of 
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Absolute systematic uncertainty 




Muon Channel 


Electron Channel 


QCD multijct 


0.001 


0.011 


Jet energy scale 


0.006 


0.012 


b tag jet energy scale 


0.001 


0.001 


Pile-up jet energy scale 


0.001 


A A AO 

0.002 


Jet reco efficiency 


0.003 


0.001 


Jet energy resolution 


A A A O 

O.OOo 


A AA1 

0.001 


Muon efficiencies 


0.001 


(n.a.) 


Muon scales / resolution 


< 0.001 


< 0.001 


Electron efficiencies 


(n.a.) 


0.001 


Electron scales / resolution 


A AA1 

0.001 


A A AO 

0.002 


b tag scale factors 


0.002 


0.004 


PDF uncertainty 


< 0.001 


< 0.001 


LAr hole uncertainty 


A AA A 

0.004 


A A A 1 

0.001 


TOT* 1 T^ri"r» 

IbR and FbR 


U.UiU 


u.uuy 


tt modelling 


(0.005) 0.011 


0.011 


Parton shower / fragmentation 


(0.004) 0.010 


0.010 


Top mass 


0.007 


(0.002) 0.007 


W+jets normalisation 


0.005 


0.007 


ly+jets shape 


< 0.001 


0.003 


Z+jets normalisation 


< 0.001 


< 0.001 


Z+jets shape 


0.001 


0.005 


Single top 


< 0.001 


< 0.001 


Diboson 


< 0.001 


< 0.001 


Charge mis-identification 


< 0.001 


< 0.001 


b tag charge 


0.001 


0.001 


MC statistics 


0.005 


0.006 


Unfolding convergence 


0.001 


0.001 


Unfolding bias 


< 0.001 


0.004 


Luminosity 


0.001 


0.001 


Combined 


0.023 


0.028 



Table 1 0. 5. : List of all systematic uncertainties taken into account for the unfolding procedure 
in the measurement of the top charge asymmetry. The numbers in brackets 
denote the uncertainties before using the larger uncertainty of both channels as 
conservative estimate. 
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Absolute systematic uncertainty 




Muon Channel 


Electron Channel 




Mtt < 450 GeV 


Mtt > 450 GcV 


Mtt < 450 GeV 


Mtt > 450 GcV 


QCD multijet 


0.018 


0.004 


0.054 


0.015 


Jet energy scale 


0.013 


0.010 


0.034 


0.006 


b tag jet energy scale 


0.001 


0.001 


0.003 


0.001 


Pile-up jet energy scale 


0.002 


0.002 


0.004 


0.002 


Jet reco efficiency 


0.002 


0.002 


0.011 


0.001 


Jet energy resolution 


0.016 


0.001 


0.036 


0.002 


Muon efiiciencies 


< 0.001 


0.001 


(n.a.) 


(n.a.) 


Muon scales / resolution 


0.002 


0.001 


0.001 


0.001 


Electron efficiencies 


(n.a.) 


(n.a.) 


0.002 


< 0.001 


Electron scales / resolution 


< 0.001 


0.001 


0.007 


0.001 


b tag scale factors 


0.001 


0.003 


0.012 


0.002 


PDF uncertainty 


0.001 


0.001 


0.001 


0.001 


LAr hole uncertainty 


0.008 


0.001 


0.008 


0.003 


ISR and FSR 


0.027 


0.024 


(0.014) 0.019 


0.019 


tt modelling 


0.018 


(0.001) 0.018 


0.032 


0.015 


Parton shower / fragm. 


(0.002) 0.006 


0.006 


(0.005) 0.006 


0.006 


Top mass 


0.019 


(0.002) 0.019 


(0.009) 0.011 


0.011 


VF+jets normalisation 


< 0.001 


0.006 


0.014 


0.004 


VF+jets shape 


< 0.001 


< 0.001 


0.002 


0.001 


Z+jets normalisation 


< 0.001 


< 0.001 


0.003 


0.001 


Z+jets shape 


0.004 


0.001 


0.017 


0.009 


Single top 


< 0.001 


< 0.001 


0.002 


< 0.001 


Diboson 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


Charge mis-identification 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


b tag charge 


0.001 


0.001 


0.001 


0.001 


MC statistics 


0.012 


0.007 


0.017 


0.009 


Unfolding convergence 


0.001 


0.001 


0.001 


0.001 


Unfolding bias 


< 0.001 


< 0.001 


0.006 


0.002 


Luminosity 


< 0.001 


0.001 


0.003 


< 0.001 


Combined 


0.049 


0.034 


0.091 


0.035 



Table 10.6.: List of all systematic uncertainties taken into account for the unfolding procedure 
in the measurement of the top charge asymmetry. The numbers in brackets 
denote the uncertainties before using the larger uncertainty of both channels as 
conservative estimate. 
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the same order of magnitude or lower than the respective MC statistics uncertainty could not 
be resolved to full extent due to the inherent fluctuations from limited statistics in the response 
matrix. In those cases, the larger of the systematic uncertainties in both channels was used 
for the final systematic uncertainty on the unfolded charge asymmetry. The resulting total 
systematic uncertainties were 0.023 in the muon+jets channel and 0.028 in the electron+jets 
channel for the inclusive unfolding. For the simultaneous unfolding in \yt\ — and M^j the 
resulting total systematic uncertainties were 0.049 (M^f < 450 GeV) and 0.034 (Mjf > 450 GeV) 
in the muon+jets channel, and 0.091 (Mjf < 450 GeV) and 0.035 {Mu > 450 GeV) in the 
electron+jets channel, respectively. 

The overall uncertainty in the individual channels was dominated by the statistical uncer- 
tainty. The systematic uncertainties were dominated by the contributions from ISR/FSR, top 
mass and jet energy resolution in the muon+jets channel, and by the uncertainties originating 
from QCD multijet background, tt modelling and parton shower / fragmentation in the elec- 
tron+jets channel. As described in Chapter O the QCD multijet background contribution has 
been estimated very conservatively, assuming a 100 % normalisation uncertainty. Most of the 
other large contributions can be traced to the available Monte Carlo statistics in the used sam- 
ples. Since only 3 000 000 {ti modelling, parton shower / fragmentation), or 1 000 000 (top mass, 
ISR/FSR) Monte Carlo events were available for the respective samples used in the evaluation 
of the systematic uncertainties (as opposed to 15 000 000 for the nominal signal sample), the 
statistical component in the evaluation of the response matrix uncertainty and the unfolding 
procedure was larger by factors of two to four with respect to the nominal case. 

In addition to the evaluation of the systematic uncertainties in the unfolding process, the 
effect of pile-up on the measured quantity before the unfolding was studied in order to ensure 
the stability of the method for different pile-up conditions. Figure [A. 91 in Appendix [Al shows 
the measured integrated asymmetry before background subtraction, AjJf*^, as a function of the 
number of primary vertices and of the bunch timing, i.e. at which relative position in the 
respective bunch the pp collision corresponding to the respective event occurred. As there was 
no statistically significant dependence of the measured asymmetry on the bunch timing or on 
the number of primary vertices in either channel, no additional systematic uncertainty due to 
pile- up was assigned. 

Both the differential and the integrated asymmetries after unfolding were in agreement with 
the MC@NLO Standard Model prediction within the estimated uncertainties. The provided 
differential distributions alongside with the respective covariance matrices can directly be put 
into context with analogue measurements at other experiments and theoretical predictions. 

10.4. Combination 

A best linear unbiased estimator (Blue) method [21311214] has been used to combine the results 
from the muon+jets and electron+jets channel after unfolding, taking into account systematic 
uncertainties and the associated correlation^. All systematic uncertainties have been consid- 
ered to be fully correlated between the muon+jets and the electron+jets channel except for 
the contributions from the Monte Carlo statistics of the response matrix, LAr defects (since 
the treatment of electrons does not affect the muon+jets channel), QCD multijet and VF+jets 
normalisation (since both have been determined with data driven methods based on orthogonal 

^Note that since the exact correlations are unknown for most of the systematic uncertainties contributing to the 
overall result, correlation coefficients between the muon+jets and electron+jets channel of either zero or one 
have been assumed. 
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datasets), charge mis- identification, and the systematic uncertainty from unfolding convergence 
and remaining bias. No correlation was assumed for the statistical uncertainties. 
For the inclusive unfolding, a combined value of 



was obtained, where the relative weight of the muon+jets channel result was 63.9%. 

The combination of the results for the simultaneous unfolding in \yt\ — lyj | and Mjf yielded 



where the relative weights of the muon+jets channel were 73.7 % and 59.5 %, respectively. Note 
that the combined systematic uncertainty was slightly lower than the systematic uncertainties 
in both the muon+jets and the electron+jets channel. This effect is inherent to the Blue 
uncertainty propagation and is due to the assumed correlations of the individual contributions, 
which can be regarded as additional prior information in the propagation of the uncertainties. 

The combined results were compatible with the mc@nlo Standard Model expectation (c.f. 
Table [TO^ within the estimated uncertainties, not indicating any significant deviation. 



A' 



unf 



0.018 ± 0.028 (stat.) ± 0.023 (syst.) 



A^'^XMti < 450 GeV) 
A^^'^iMa > 450 GeV) 



0.053 ± 0.070 (stat.) ± 0.054 (syst.) 
0.008 ± 0.035 (stat.) ± 0.032 (syst.) 
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11. Summary &; Conclusion 



A measurement of the charge asymmetry in the production of top quark pairs at the ATLAS 
experiment was performed, using a dataset corresponding to an integrated luminosity of 1.04 fb~^ 
taken over the course of 2011 at a centre-of-mass energy of ^/s = 7TeV. 

An object and event selection was employed in the lepton+jets decay channel in order to 
identify events with a signature corresponding to a semileptonic decay of a tt pair, given by one 
isolated lepton (muon or electron) with large transverse momentum, at least four reconstructed 
jets and large missing transverse momentum. The selection was furthermore chosen such that 
various background contributions, including the production of single top quarks, heavy gauge 
bosons in association with jets, the contribution from diboson production and fake leptons pre- 
dominantly produced in QCD multijet events, were reduced. A set of background contributions 
was estimated using Monte Carlo simulations, while the W^+jets and QCD multijet normalisa- 
tion was determined using data driven methods. 

A kinematic fit was performed in order to reconstruct tt events based on the measured objects 
after application of the selection, yielding the most probable object kinematics at the parton 
level under the assumption of a semileptonic top quark decay event topology. 

The reconstructed objects were used to obtain the differential distribution of the difference of 
absolute rapidities of the reconstructed top and antitop, |?/(| — A subtraction of the various 
background contributions was performed in order to obtain the integrated charge asymmetry of 
the ti signal contribution at reconstruction level, Aq^°. Furthermore, A^^^° was determined as 
a function of the invariant tt mass, M^j. 

An unfolding procedure was applied to the reconstructed \yt\ — \yt\ distribution in order to 
correct for resolution and acceptance effects and to obtain the corresponding distribution at par- 
ton level. A Bayesian iterative unfolding procedure was used and cross-checks and calibrations 
were performed to verify the linearity, convergence and stability of the approach. The unfolding 
was performed both for the inclusive \yt\ — \yt\ distribution and in addition simultaneously in 
\yt\ ~ Ivtl ^'^d M(f. A combination of obtained results in the muon-l-jets and electron+jets chan- 
nel was performed using the Blue method, yielding a combined integrated charge asymmetry 
after unfolding of 

yl^'if = -0.018 ± 0.028 (stat.) ± 0.023 (syst.) , (11.1) 

while the MC@NLO prediction was 0.0056 it 0.0003 (stat.). For the simultaneous unfolding in 

\yt\ — lutl and Mjf, combined values of 

A'^^'^Mtt < 450 GeV) = -0.053 ± 0.070 (stat.) ± 0.054 (syst.) , 
A'^^Mu > 450 GeV) = -0.008 ± 0.035 (stat.) ± 0.032 (syst.) 

were obtained. The mc@nlo predictions were 0.0024±0.0004 (stat.) and 0.0086±0.0004 (stat.) , 
respectively. A summary of the results obtained in the simultaneous unfolding in \yt\ — |yf| and 
Mjf can be found in Figure 111.11 

The obtained results were in agreement with the Standard Model prediction. However, both 
for the measured asymmetry after reconstruction and the corresponding unfolded values a ten- 
dency to more negative integrated asymmetries was observed in the electron channel. This effect 
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Figure 11.1.: Unfolded asymmetries in two regions of Mu eompared to the prediction from 
MC@NLO |183] . The error bands on the mc@nlo prediction include uncertainties 
from parton distribution functions and renormalisation and factorisation scales. 



is most dominant for the simultaneous unfolding in the region where M^f < 450 GeV. Neverthe- 
less, the discrepancy is covered by the estimated overall uncertainties including systematics. 

The estimated uncertainties were dominated by the statistical uncertainty, in particular for 
the simultaneous unfolding in \yt\ — \yt\ and Mjf. The largest contributions to the systematic 
uncertainties originated from ISR/FSR, top mass and jet energy resolution in the muon+jets 
channel, and from the uncertainties from QCD multijet background, tt modelling and parton 
shower/fragmentation in the electron+jets channel. Most of these contributions, except for the 
QCD multijet background uncertainty, can be attributed to large parts to the statistical com- 
ponent inherent in the evaluation of the systematic uncertainties, which involves a replacement 
of the unfolding response matrix with a corresponding matrix obtained from a different sample. 
Since only a fraction of events of about 10 % with respect to the nominal ti sample was available 
in those cases, statistical fluctuations in the response matrix dominated over the actual system- 
atic shift to be evaluated. This also explains the discrepancies observed in the comparison of 
the individual contributions between the muon+jets and electron+jets channel. Furthermore, 
the QCD multijet normalisation uncertainty was conservatively estimated to be 100 % following 
the recommendations of the performance groups, which led to a large contribution from this 
uncertainty, in particular in the electron+jets channel. 

A summary of the obtained results alongside with recent results obtained by D0, CDF 
and CMS in comparison to several theoretical models beyond the Standard Model is shown 
in Figure [11. 21 Potential regions in the phase space of inclusive charge asymmetry from new 
physics, indicated by the variable at the LHC plotted against the associated forward- 

backward asymmetry at the Tevatron, are highlighted (c.f. Chapter 12.3. 2p . The mea- 

surement performed in this analysis together with the results from other experiments already 
puts pressure on several of the proposed models. This is most prominent for Z' models, which 
are disfavoured at 2(T to 3(T by this measurement and a similar measurement by CMS, while 
being favoured in the high invariant mass region by CDF. Furthermore, disagreements between 
favoured regions from the measurements at the Tevatron and at the LHC are observed, in par- 
ticular for the high invariant ti mass region. These measurements of the charge asymmetry at 
hadron colliders provide the first step towards a better understanding of Quantum Chromody- 
namics and possible extensions of the Standard Model. 

Despite the limitation of the measurement by the statistical uncertainty, in particular for the 
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Figure 11.2.: Measurement results and predicted charge asymmetries at the Tevatron and 
LHC for various BSM models [11611127] . The inclusive charge asymmetry from 
new physics A^q"" at the LHC vs. corresponding the forward-backward asymme- 
try at the Tevatron (left) and the identical predictions in a high invariant 
mass region where Mu > 450 GeV (right) for the different models in the created 
phase space is shown. The labelled lines indicate the central values of the results 
measured at different experiments, while the dashed lines indicate the 1 a uncer- 
tainty regions (gray area denotes the ATLAS measurement uncertainty). The 
Standard Model prediction corresponds to Ap'^ = and A^"" = 0, respectively. 

simultaneous unfolding in \yt\ — |yf| and M^f, an increased dataset from the ATLAS experiment 
will quickly reduce this uncertainty, significantly increasing the sensitivity. This will allow to 
make much more precise statements especially in regions where an increased asymmetry is 
expected, making it possible to fully exclude several of the available BSM models with sufficient 
significance. Furthermore, non-excluded models could be further constrained in their respective 
parameters. 

Further improvements can be achieved for the systematic uncertainties, where methods will 
evolve and in particular existing estimates can be replaced by more precise studies and methods. 
As an example, the QCD multijet systematic uncertainty has been evaluated in Appendix [D] us- 
ing more advanced methods, as described in Chapter 16.3.11 However, one of the most prominent 
sources of uncertainty is the limited statistics of the used Monte Carlo samples, in particular 
for the evaluation of the systematic uncertainties. These could be reduced significantly if larger 
samples would be available. 

Optimisations in the unfolding procedure could lead to further reductions of the statistical 
uncertainties. Further studies of the unfolding procedure and the associated statistical and 
systematic uncertainties are discussed in Appendix [Dl 

Finally, an additional cross-check has been performed for the inclusive charge asymmetry, 
comparing the obtained results to an SVD unfolding approach, using otherwise identical param- 
eters. Comparable values have been obtained. 
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A. Additional Control Plots 
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Figure A.I.: Relative resolution of Ad^i for the kinematic event reconstruction without any 
requirement on the reconstruction logarithmic likelihood log L for the muon+jets 
channel (left) and the electron+jcts channel (right). 



d 400 - ATLAS Work in progress 
i 350^ 1'-'" = 1-04 



n + > 4 jets {> 1 b tag)r 



300 r 
250 r 
200 r 
150r 
100r 
50 r 



Mean = -0.033 — 
c = 0.184 



-^.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



° 220 
o 

« 200 
I 180 

UJ 

160 
140 
120 
100 
80 
60 
40 
20 



_ ATLAS Work in progress 

Ldt = 1.04fb' 



- AIL 



e + > 4 jets (> 1 b tag) : 



Mean = -0.033 
c = 0.177 



i.5 -2 -1.5 -1 -0.5 0.5 1 



1.5 2 2.5 



Figure A. 2.: Relative resolution of Alfi for the kinematic event reconstruction after application 
of a requirement on the reconstruction logarithmic likelihood of log L > —52. The 
relative resolution improves from 28.8 % to 18.4 % in the muon+jets channel (left) 
and from 28.8% to 17.7% in the electron+jets channel (right). 
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A. Additional Control Plots 
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Figure A. 3.: Control plots for the muon+jets channel prior to the requirement of at least 
one h tagged jet. From the top left to the bottom right, the transverse momen- 
tum pt, the pseudorapidity rj and the azimuthal angle 4> of the selected muon 
are shown. Additional plots show the transverse momentum of the leading jet 
(pT(ji)), the W transverse mass rmiW) and the transverse missing energy _^t- 
Uncertainties are statistical and for M^+jets also include systematic uncertainties 
on normalisation. For the QCD multijet background, a conservative 100 % sys- 
tematic uncertainty was assumed. In addition, the uncertainties on luminosity, 
jet energy scale and ti cross-section are shown. 
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Figure A. 4-: Control plots for the electron+jets channel prior to the requirement of at least 
one b tagged jet. From the top left to the bottom right, the transverse momen- 
tum pt, the pseudorapidity rj and the azimuthal angle of the selected electron 
are shown. Additional plots show the transverse momentum of the leading jet 
(ptOi)), the W transverse mass mT(VF) and the transverse missing energy _^t- 
Uncertainties are statistical and for VF+jets also include systematic uncertainties 
on normalisation. For the QCD multijet background, a conservative 100% sys- 
tematic uncertainty was assumed. In addition, the uncertainties on luminosity, 
jet energy scale and tt cross-section are shown. 
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A. Additional Control Plots 
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Figure A. 5.: Expected statistical uncertainties on the unfolded value of as a function 
of the regularisation parameter A^it for different injected true asymmetries in 
the observable A^"^ for both the muon+jets channel (left) and electron+jets 
channel (right). The top row shows the respective distribution for the inclu- 
sive measurement, while the lower rows show the corresponding distributions for 
M^-t < 450 GeV and M^t > 450 GeV, respectively. For the simultaneous unfold- 
ing in \yt \ — \yt\ and M^^, a cut on the event reconstruction likelihood logL was 
applied to improve the Mff resolution of the selected events. The arrows indicate 
the chosen values for A^it. 



102 



^ l^ftL^Tr'ueAjJ-ollOs' 
■g -.-True A^ = -0.051 

Di - ^True A,, - -0.026 
£ - True A^ = -0.000 
g 120- True A„ = 0.025 

^ - ♦True A„ = 0.050 
*True A. = 0.101 



100 



-«-True A, 

i 



ATLAS Work in progress 

H + > 4 jets {> 1 b tag) 



60- 



40- 



20- 



qI ' i T 

T) 10 20 30 40 50 60 70 80 90 

N„ 




20 30 40 50 



60 70 80 90 
N„ 



^ -.-True = -0.051 

p - .«.True A^ = -0.026 
% - True Ap = -0.000 
g 120- True A. = 0.025 
^ - ♦True A. = 0.050 
' ♦True Aj, = 0.101 

100 

80- I • 

60- • 

40- 
20- I 



ATLAS Work in progress 

H > 4 jets (> 1 b tag) : 



I 



T 

T) 10 20 30 40 50 60 70 80 90 

N„ 



1 

1 120 
u 
Z 

100 
80 
60 
40 
20 



_ ♦True A^ = -0.103 

♦ True A^ = -0.051 

♦ True A^ = -0.026 
True Ap = -0.000 
True Aj. = 0.025 

♦ True A^ = 0.050 

♦ True Ap = 0.101 



ATLAS Work in progress 

e + > 4 jets (> 1 b tag) 



10 20 30 40 50 



60 70 80 90 
N„ 



^ ♦True A^ = -0.051 

g> - ♦True A^ = -0.026 
% - True Ap = -0.000 
i 120- True A. = 0.025 
^ - ♦True A. = 0.050 
' ♦True A. = 0.101 



ATLAS Work in progress 

M > 4 jets (> 1 b tag) : 



100- 
80- 
60- 
401, 
20-1 



140 
120 
100 



_ ■•■True A,, = -0.103 
-c-True A,, - -0.051 
■•■True A,, - -0.026 
True Ap = -0.000 
True Aj. = 0.025 
-•-True A,, = 0.050 
■•■True Ap - 0.101 



10 20 30 40 50 60 70 80 90 



ATLAS Work in progress 

e + > 4 jets (> 1 b tag) 




Figure A. 6.: Distributions showing the convergence behaviour of ensembles of pscudodata as 
a function of regularisation for different injected true asymmetries A^"® for both 
the muon+jets channel (left) and electron+jets channel (right), parametrised by 
the number of iterations A^it in the unfolding procedure. The relative amount 
of ensembles which have reached a convergent state at given TVh (100% being 
indicated by the red dashed line) are shown. The top row shows the respective 
distributions for the inclusive measurement, while the lower rows show the corre- 
sponding distributions for M^^ < 450 GeV and Mj^ > 450 GeV, respectively. For 
the simultaneous unfolding in \yt \ — \yi\ and Mjf, a cut on the event reconstruc- 
tion likelihood logL was applied to improve the M^i resolution of the selected 
events. The arrows indicate the chosen values for A^it. 
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Figure A. 7. 



Distributions showing the standard deviation of the variable AA™^ for ensembles 
of pseudodata as a function of regularisation for both the muon+jets channel 
(left) and electron+jets channel (right), parametrised by the number of itera- 
tions Nit in the unfolding procedure. The applied convergence criterion of 0.001 
is indicated by the red dashed line. The top row shows the respective dependen- 
cies for the inclusive measurement, while the lower rows show the corresponding 
dependencies for Mjj < 450 GeV and Mjj > 450 GeV, respectively. For the si- 
multaneous unfolding in \yt \ — \yt\ and Af^j, a cut on the event reconstruction 
likelihood log L was applied to improve the Mtt resolution of the selected events. 
The arrows indicate the chosen values for iVit . 
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Figure A. 8.: Pull distributions showing the difference between true and unfolded asymmetries, 
normalised to the unfolding respective uncertainty, obtained from pseudoexper- 
iments for both the muon+jets channel (left) and electron+jcts channel (right), 
parametrised by the number of iterations iVit in the unfolding procedure for a 
true asymmetry of A^^'' = 0.0. The top row shows the respective dependen- 
cies for the inclusive measurement, while the lower rows show the corresponding 
dependencies for M^j < 450 GeV and M^t > 450 GeV, respectively. For the si- 
multaneous unfolding in \yt\ — \yt\ and Mtt, a cut on the event reconstruction 
likelihood log L was applied to improve the M^j resolution of the selected events. 
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A. Additional Control Plots 
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Figure A. 9.: Measured charge asymmetry A'^f^'^ before background subtraction and unfolding 
for both the muon+jcts channel (left) and elcctron+jets channel (right) as a 
function of the number of primary vertices in the respective events (top row) and 
the bunch position (bottom row). Uncertainties are statistical only. 
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B. Systematic Uncertainties 





Absolute systematic uncertainty 










^ 001 
\ u.uux 


OOB 




U.UUo 


U.UU4 


b tag jet energy scale 


< U.UUi 


U.UUi 


Pilc-up jet energy scale 


<r 001 


<^ 001 


jcb lecu L,iiieieiic_y 


0.001 


0.002 


Jet energy resolution 


0.004 


0.001 


Muon efficiencies 


0.001 


(n.a.) 


Muon scales / resolution 


< 0.001 


< 0.001 


Electron efficiencies 


(n.a.) 


0.001 


Electron scales / resolution 


< 0.001 


0.001 


b tag scale factors 


0.002 


0.003 


LAr hole uncertainty 


0.002 


0.001 


W+jets normalisation 


0.003 


0.004 


W+jets shape 


< 0.001 


0.002 


Z+jcts normalisation 


< 0.001 


< 0.001 


Z+jcts shape 


< 0.001 


< 0.001 


Single top 


< 0.001 


< 0.001 


Diboson 


< 0.001 


< 0.001 


Charge mis-identification 


< 0.001 


< 0.001 


6 tag charge 


0.001 


0.001 


Luminosity 


0.001 


0.001 


Combined 


0.008 


0.010 



Table B.I.: List of all systematic uncertainties taken into account for the measurement of the 
top charge asymmetry before unfolding. 
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B. Systematic Uncertainties 





Absolute systematic uncertainty 




Muon Channel 


Electron Channel 




Mti < 450 GeV 


Mti > 450 GcV 


Mti < 450 GeV 


Mti > 450 GcV 


QCD multijct 


< 0.001 


< 0.001 


0.013 


0.013 


Jet energy scale 


0.003 


0.004 


0.008 


0.006 


b tag jet energy scale 


< 0.001 


< 0.001 


0.001 


0.001 


Pile-up jet energy scale 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


Jet reco efSciency 


< 0.001 


0.001 


0.003 


0.001 


Jet energy resolution 


0.006 


0.003 


0.009 


0.003 


Muon efficiencies 


< 0.001 


0.001 


(n.a.) 


(n.a.) 


Muon scales / resolution 


< 0.001 


0.001 


< 0.001 


< 0.001 


Electron efficiencies 


(n.a.) 


(n.a.) 


0.001 


< 0.001 


Electron scales / resolution 


< 0.001 


< 0.001 


0.001 


0.001 


b tag scale factors 


< 0.001 


0.001 


0.003 


0.002 


LAr hole uncertainty 


0.001 


0.001 


0.002 


0.001 


VF+jets normalisation 


0.001 


0.004 


0.004 


0.003 


VF+jets shape 


< 0.001 


< 0.001 


< 0.001 


0.001 


Z+iets normalisation 


< 0.001 


< 0.001 


0.001 


< 0.001 


Z+jets shape 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


Single top 


< 0.001 


< 0.001 


0.001 


< 0.001 


Diboson 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


Charge mis-identification 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


b tag charge 


0.001 


0.001 


0.001 


0.001 


Luminosity 


< 0.001 


0.001 


0.001 


< 0.001 


Combined 


0.007 


0.007 


0.019 


0.015 



Table B. 2. : List of all systematic uncertainties taken into account for the measurement of the 
top charge asymmetry before unfolding. 
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C. SVD Unfolding 



As an additional cross-check, the inclusive unfolding has been performed with the SVD unfolding 
procedure (for details, refer to Chapter [8|), using otherwise identical analysis parameters. The 
same binning and the same response matrix as determined for the Bayesian iterative unfolding 
(c.f. Chapter I10.3P has been used. 

A closure test has been performed to verify that the SVD unfolding approach can be used 
to recover an arbitrary asymmetry present in the true distribution. Similar to the Bayesian 
iterative unfolding, pseudoexperiments using statistically independent sets of pseudodata cor- 
responding to the statistics expected in the used data sample were performed. Ensemble tests 
were conducted to confirm the linearity of the unfolding response in the true value of the cho- 
sen charge asymmetry observable, and the dependency of the obtained results on the 
regularisation of the unfolding procedure was studied. 

The average unfolded value of obtained from the respective sets of pseudodata as a 
function of the injected true value of Aq^^ can be found in Figure ICTT] for different regularisation 
strengths corresponding to choices of the regularisation parameter of r = 2, 3, 4, 5 and 6. Since 
r corresponds directly to a fraction of the number of bins, it can not exceed the number of 
bins chosen for the unfolding by construction |193j . A straight line fit has been performed as 




True True A^, 

Figure C.I.: The obtained inclusive asymmetry after SVD unfolding as a function of 

injected true asymmetry for different regularisation parameters r for both 

the muon+jets channel (left) and electron+jets channel (right). 

described in Chapter 110.31 As can be seen in the respective fits, a slope close to one can be 
achieved in all cases for r = 6, indicating a proper correspondence of the unfolded asymmetry 
to the respective true value, independent of the strength of the injected asymmetry. 

Since the SVD unfolding does not rely on an iterative procedure, no convergence criterion was 
imposed to choose a proper regularisation. Instead, the parameter r was chosen such that the 
remaining unfolding bias was minimised, i.e. such that the slope parameter a of the linear fit 
in the calibration was as close to unity as possible. As can be seen in Figure ICTH the optimal 
choice was given by r = 6 for both the muon+jets and electron+jets channel. 
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C. SVD Unfolding 



Additional closure tests have been performed as described in Chapter 110.31 using ensembles 
of pseudodata to quantify any remaining bias from the unfolding at the chosen regularisation 
strengths. A corresponding systematic uncertainty was assigned to the unfolded result accord- 
ingly from the residuals of the respective pull distributions (c.f. Chapter [9]). 

The unfolded integrated asymmetries for the inclusive unfolding, using r = 6 in the SVD 
unfolding process, were determined to be 

A'^^ = -0.005 ± 0.034 (stat.) ± 0.024 (syst.) 

in the muon channel and 

A^°f = -0.056 ± 0.043 (stat.) ± 0.029 (syst.) 

in the electron channel, respectively. In addition, the obtained distributions after unfolding can 
be found in Figure IC.2[ 
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Figure C.2.: Unfolded distribution of \yt \ — \yt \ using an SVD unfolding procedure, normalised 
to unity for both the muon+jets channel (left) and electron+jets channel (right). 
The uncertainties include both statistical and systematic shape components. 



A summarised list of all systematics and their contribution to the overall systematic uncer- 
tainty can be found in Table IUT] All systematic uncertainties based on a replacement of the 
unfolding matrix which were of the same order of magnitude or lower than the respective MC 
statistics uncertainty could not be resolved to full extent due to the inherent fluctuations from 
limited statistics in the response matrix. In those cases, the largest of the systematic uncer- 
tainties in both channels was used for the final systematic uncertainty on the unfolded charge 
asymmetry. The resulting combined systematic uncertainties were 0.024 in the muon-l-jets chan- 
nel and 0.029 in the electron -|-jets channel for the inclusive unfolding. 

The Blue method was used to combine the measurements after performing the SVD unfolding 
in the muon-|-jets and electron-|-jets channel. The same assumptions about correlations as 
described in Chapter 110.41 were used. A combined value of 

= -0.024 ± 0.027 (stat.) ± 0.024 (syst.) 

was obtained, where the relative weight of the muon-l-jets channel result was 64.3%. 

The obtained results were compatible with the results from the Bayesian unfolding shown in 
Chapter fnOl 
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Absolute systematic uncertainty 




Muon Channel 


Electron Channel 


QCD multijet 


0.002 


0.012 


Jet energy scale 


0.007 


0.013 


b tag jet energy scale 


< 0.001 


< 0.001 


Pile-up jet energy scale 


0.001 


0.002 


Jet reco efficiency 


0.003 


0.001 


Jet energy resolution 


0.008 


0.004 


Muon efficiencies 


0.001 


(n.a.) 


Muon scales / resolution 


0.001 


< 0.001 


Electron efficiencies 


(n.a.) 


0.001 


Electron scales / resolution 


< 0.001 


0.002 


b tag scale factors 


0.003 


0.005 


PDF uncertainty 


< 0.001 


< 0.001 


LAr hole uncertainty 


0.004 


0.001 


ISR and FSR 


0.011 


0.008 


ti modelling 


(0.008) 0.011 


0.011 


Parton shower / fragmentation 


(0.001) 0.010 


0.010 


Top mass 


0.006 


(0.002) 0.006 


M^+jets normalisation 


0.005 


0.008 


M^+jets shape 


< 0.001 


0.004 


Z+jets normalisation 


< 0.001 


< 0.001 


Z+jets shape 


0.001 


0.005 


Single top 


< 0.001 


< 0.001 


Diboson 


< 0.001 


< 0.001 


Charge mis-identification 


< 0.001 


< 0.001 


b tag charge 


0.001 


0.001 


MC statistics 


0.004 


0.006 


Unfolding bias 


0.001 


0.004 


Luminosity 


0.001 


0.001 


Combined 


0.024 


0.029 



Table C.I.: List of all systematic uncertainties taken into account for the unfolding proce- 
dure in the measurement of the top charge asymmetry. The numbers in brackets 
denote the uncertainties before using the largest uncertainty of both channels as 
conservative estimate. 



Ill 



D. Additional Studies of Unfolding and 
Systematics 



As described in Chapter 110.31 the unfolding procedure employed used a strong regularisation in 
order to achieve a convergent state and to reduce the remaining bias after unfolding. However, 
this implies large statistical uncertainties on the obtained results. Furthermore, as described in 
Chapter O the QCD multijet background contribution systematic uncertainty was conservatively 
assumed to be 100 % and no explicit shape uncertainty was included despite the availability of 
more advanced estimates. 

As an additional cross-check, the unfolding procedure has been repeated, discarding the ex- 
plicit requirement of convergence of the Bayesian iterative procedure. Instead, smaller values of 
A^it were used, implying lower statistical uncertainties, and the expected remaining bias of the 
unfolded value A'^^ was extracted from pseudoexperiments and corrected for. This, however, 
implies larger assumptions about the underlying physics since the remaining bias is determined 
and corrected for based on the observable Ac only. Consequently, no fully model-independent 
calibration can be performed and a bias not covered by the uncertainties can remain. Further- 
more, the described approach does not allow for the extraction of the full unfolded distribution 
and the associated bin-by-bin uncertainties since the calibration is performed with respect to 
the integrated asymmetry observable Ac only. 

A closure test has been performed to calibrate the unfolded asymmetry A^^ with respect to 
the true asymmetry ^4^"^. Similar to the procedure described in Chapter 110.31 pseudoexperi- 
ments using statistically independent sets of pseudodata corresponding to the statistics expected 
in the used data sample (after background subtraction) were created based on a Poissonian fluc- 
tuation of the respective reconstructed distributions of \yt\ — \yt\. Ensemble tests were performed 
to confirm the linearity of the unfolding response in the true value of the chosen charge asym- 
metry observable, and the dependence of the obtained results on the regularisation of the 
unfolding procedure was studied. 

The average value of A^^ obtained from the respective sets of pseudodata as a function 
of the injected true asymmetry, A^"®, can be found for both the inclusive unfolding and the 
simultaneous unfolding in \yt \ — \yt\ and M^f in Figure iD. II for different regularisation strengths, 
using A'lt = 2, 4, 6, 7, 8 and 10 iterations in the Bayesian unfolding. A straight line fit has been 
performed as described in Chapter 110.31 As can be seen in the respective fits, the slopes differ 
from one, in particular for M^j < 450. 

In addition, the respective expected statistical uncertainty on A^^ is shown in Figure lD.2l as a 
function of Nn, taking into account the parameters for the slope and offset a and b, respectively, 
from the straight line fits of the calibration curves. For such small changes in Nn, the expected 
statistical uncertainty after correction increases only slightly. 

The number of iterations for the unfolding of the used data set was chosen such that the 
statistical uncertainty was minimised, while requiring that any injected asymmetry in the tested 
range was recoverable within the expected statistical uncertainties. This approach, however, led 
to a remaining bias of the unfolded asymmetry Aq^^ with respect to the true asymmetry A^"'^. 
Consequently, this bias was corrected for by taking into account the slope and offset parameters a 
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True A(, True A^, 



Figure D.I.: The obtained overall inclusive asymmetry after unfolding, A^^, as a function 
of injected true asymmetry Aq^'^ for different regularisation parameters Nu for 
both the muon+jets channel (left) and electron+jets channel (right). The top 
row shows the respective distribution for the inclusive measurement, while the 
lower rows show the corresponding distributions for M^f < 450 GeV and M^f > 
450 GeV, respectively. For the simultaneous unfolding in \yt \ — \yt\ and M^j, a 
cut on the event reconstruction likelihood log L was applied to improve the M^^ 
resolution of the selected events. 
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Figure D.2.: Expected statistical uncertainties on as a function of the regularisation pa- 
rameter Nit for different injected true asymmetries A^"® for both the muon+jets 
channel (left) and electron+jets channel (right). The top row shows the respec- 
tive distribution for the inclusive measurement, while the lower rows show the 
corresponding distributions for M^f < 450 GeV and M^^ > 450 GeV, respectively. 
For the simultaneous unfolding in \yt \ — \yt\ and -Mjf, a cut on the event re- 
construction likelihood logL was applied to improve the M^-i resolution of the 
selected events. The arrows indicate the chosen values for iVit. 
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and b of the straight hne fit, respectively, to obtain a corrected value of the unfolded asymmetry, 

A unf,corr • i 

' gi'^^n by 

I unf 



/lunl _ h 

^unfjCorr 

~ a ' 



(D.l) 



Consequently, the corresponding statistical and systematic uncertainties after correction, cr^^''^"^^ , 
can be determined by 

_unf 



unf,corr '^Ap 



a 



(D.2) 



Following the described requirements, regularisation parameters of A'^it = 6 for the muon+jets 
and Nit = 4 in the electron+jets channel were chosen for the inclusive unfolding, while for the 
simultaneous unfolding in \yt \ — \yt\ and M^f, A^it = 7 and Nn = 6 were chosen for the muon+jets 
and electron+jets channel, respectively. The corresponding slope and offset parameters extracted 
from the straight line fit in the calibration can be found in Table ID. II for completeness. 



Unfolding 


muon+jets channel 


electron+jets channel 




a (slope) 


b (offset) 


Nu 


a (slope) 


6 (offset) 


Inclusive 


6 


0.841 ± 0.002 


-0.0009 ± 0.0001 


4 


0.739 ± 0.002 


0.0004 ± 0.0001 


Mti < 450 GcV 


7 


0.711 ± 0.004 


-0.0021 ± 0.0002 


6 


0.664 ± 0.004 


0.0012 ± 0.0002 


Mtt > 450 GeV 


7 


0.974 ± 0.003 


0.0007 ± 0.0002 


6 


0.902 ± 0.003 


-0.0032 ± 0.0002 



Table D.l.: Slopes and offsets from the linear fit in the unfolding calibration. The parameters 
were obtained for a linear fit of the average unfolded value of Aq^^ as a function of 
the true value obtained from sets of pseudoexperiments. For the simulta- 

neous unfolding in \yt\ — \yt\ and Mjj, a cut on the event reconstruction likelihood 
log L was applied to improve the M^j resolution of the selected events. 



The unfolded integrated asymmetries in A'^^ are shown in Table ID.2I for both the inclusive 
unfolding and for the simultaneous unfolding in \ yt\ — \yt\ and M^f, using the chosen regularisation 
strength in the unfolding process. 



Unfolding 


Observed 


Predicted (mc@nlo) 


Muon+jets Channel 


Inclusive 

Mu < 450 GeV 

Mtt > 450 GeV 


-0.020 ± 0.027 (stat.) ± 0.022 (syst.) 
-0.025 ± 0.054 (stat.) ± 0.034 (syst.) 
-0.015 ± 0.035 (stat.) ± 0.021 (syst.) 


0.0056 ± 0.0003 (stat.) 
0.0024 ± 0.0004 (stat.) 
0.0086 ± 0.0004 (stat.) 


Electron+jets Channel 


Inclusive 

Mtt < 450 GeV 
Mtt > 450 GeV 


-0.063 ± 0.035 (stat.) ± 0.029 (syst.) 
-0.122 ± 0.069 (stat.) ± 0.046 (syst.) 
-0.046 ± 0.046 (stat.) ± 0.034 (syst.) 


0.0056 ± 0.0003 (stat.) 
0.0024 + 0.0004 (stat.) 
0.0086 ± 0.0004 (stat.) 



Table D.2.: Unfolded values of the charge asymmetry observable A^^^ for the muon+jets and 
electron+jets channel. The results for the inclusive measurement and the respec- 
tive results for the simultaneous unfolding in \yt\ — \yt \ and M^f for M^f < 450 GeV 
and Mtt > 450 GeV, taking into account the correction for unfolding bias, are 
shown. For the simultaneous unfolding, a cut on the event reconstruction like- 
lihood logL was applied to improve the Mtt resolution of the selected events. 
Furthermore, the respective mc@nlo predictions are shown. 
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D. Additional Studies of Unfolding and Systematics 



A summarised list of all systematics and their contribution to the overall systematic uncertain- 
ties can be found in Table IDTsI for the inclusive unfolding and in Table ID^ for the simultaneous 
unfolding in \yt\ — \yt\ and 71%. 





Absolute systematic uncertainty 




Muon Channel 


Electron Channel 


QCD multijct 


0.005 


0.011 


Jet energy scale 


0.008 


0.013 


b tag jet energy scale 


< 0.001 


< 0.001 


Pile-up jet energy scale 


0.001 


0.002 


Jet reco efficiency 


0.002 


0.001 


Jet energy resolution 


0.009 


0.003 


Muon efficiencies 


0.002 


(n.a.) 


Muon scales / resolution 


0.001 


0.001 


Electron efficiencies 


(n.a.) 


0.001 


Electron scales / resolution 


0.001 


0.002 


b tag scale factors 


0.004 


0.006 


PDF uncertainty 


< 0.001 


< 0.001 


LAr hole uncertainty 


0.003 


< 0.001 


ISR and FSR 


0.005 


0.013 


tt modelling 


(0.002) 0.007 


0.007 


Parton shower / fragmentation 


(0.004) 0.011 


0.011 


Top mass 


0.004 


0.002 


M^+jets normalisation 


0.007 


0.010 


ly+jets shape 


< 0.001 


0.005 


^-|-jets normalisation 


< 0.001 


< 0.001 


Z+jets shape 


0.002 


0.003 


Single top 


< 0.001 


< 0.001 


Diboson 


< 0.001 


< 0.001 


Charge mis-identification 


< 0.001 


< 0.001 


b tag charge 


0.001 


0.001 


MC statistics 


0.004 


0.005 


Luminosity 


0.002 


0.002 


Combined 


0.022 


0.029 



Table D.3.: List of all systematic uncertainties taken into account for the unfolding proce- 
dure in the measurement of the top charge asymmetry. The numbers in brackets 
denote the uncertainties before using the larger uncertainty of both channels as 
conservative estimate. 

All systematic uncertainties based on a replacement of the unfolding matrix which were of 
the same order of magnitude or lower than the respective MC statistics uncertainty could not 
be resolved to full extent due to the inherent fluctuations from limited statistics in the response 
matrix. Hence, for those cases, the largest of the systematic uncertainties in both channels 
was used for the final systematic uncertainty on the unfolded charge asymmetry. The resulting 
combined systematic uncertainties were 0.022 in the muon-l-jets channel and 0.029 in the elec- 
tron-l-jets channel for the inclusive unfolding. For the simultaneous unfolding in \yt\ — \yt\ and M^f 
the resulting total systematics were 0.034 (M^f < 450 GeV) and 0.021 (M^j > 450 GeV) in the 
muon-Fjets channel, and 0.046 (M^f < 450 GeV) and 0.034 {Mu > 450 GeV) in the electron-l-jets 
channel, respectively. 

The Blue method was used to combine the measurements after performing the unfolding 



116 





Absolute systematic uncertainty 




Muon Channel 


Electron Channel 




Mft < 450 GeV 


Mft > 450 GeV 


Mtt < 450 GeV 


Mil > 450 GeV 


QCD multijet 


0.010 


0.002 


0.021 


0.020 


Jet energy scale 


0.007 


0.009 


0.021 


0.011 


b tag jet energy scale 


0.001 


0.001 


0.001 


0.002 


rile-up jet energy scale 


0.001 


0.001 


0.002 


0.001 


Jet reco efBciency 


< 0.001 


0.001 


0.007 


0.001 


Jet energy resolution 


0.009 


0.005 


0.011 


0.007 


Muon efficiencies 


0.001 


0.001 


(n.a.) 


(n.a.) 


Muon scales / resolution 


0.002 


0.001 


< 0.001 


0.001 


Electron efficiencies 


(n.a.) 


(n.a.) 


0.001 


0.001 


Electron scales / resolution 


< 0.001 


0.001 


0.005 


0.003 


b tag scale factors 


0.001 


0.003 


0.007 


0.004 


PDF uncertainty 


0.001 


0.001 


0.001 


0.002 


LAr hole uncertainty 


0.001 


0.001 


0.002 


0.002 


ISR and FSR 


0.021 


0.007 


0.024 


0.018 


tt modelling 


0.016 


(0.004) 0.016 


0.005 


(0.003) 0.005 


Parton shower / fragm. 


(0.003) 0.006 


0.006 


(0.001) 0.012 


0.012 


Top mass 


0.010 


(0.003) 0.010 


0.003 


0.005 


VF+jets normalisation 


0.003 


0.006 


0.009 


0.006 


VF+jets shape 


^ u.uu± 




n nni 


n nni 


Z+jets normalisation 


< 0.001 


< 0.001 


0.001 


< 0.001 




0.007 


< 0.001 


0.013 


0.003 




< 0.001 


< 0.001 


0.001 


< 0.001 


Diboson 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


Charge mis-identification 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


b tag charge 


0.001 


0.001 


0.001 


0.001 


MC statistics 


0.008 


0.005 


0.010 


0.007 


Luminosity 


0.001 


0.001 


0.002 


0.001 


Combined 


0.034 


0.021 


0.046 


0.034 



Table D.4-: List of all systematic uncertainties taken into account for the unfolding proce- 
dure in the measurement of the top charge asymmetry. The numbers in brackets 
denote the uncertainties before using the larger uncertainty of both channels as 
conservative estimate. 

for the described alternative approach in the muon-l-jets and electron -|-jets channel, taking into 
account the respective systematic uncertainties and the associated correlations. The same as- 
sumptions about correlations as described in Chapter 110.41 were used and a combined value 
of 

A'^^ = -0.035 ± 0.021 (stat.) ± 0.021 (syst.) 

for the inclusive unfolding was obtained, where the relative weight of the muon-l-jets channel 
result was 65.7%. The combination of the results for the simultaneous unfolding in \yt\ — |yf| 
and M^f yielded 

A'^^'^Mu < 450 GeV) = -0.058 ± 0.043 (stat.) ± 0.035 (syst.) , 
A^^Mu > 450 GeV) = -0.025 ± 0.028 (stat.) ± 0.022 (syst.) , 

where the relative weight of the muon-l-jets channel was 65.6 % and 69.4 %, respectively. 
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D. Additional Studies of Unfolding and Systematics 



As expected, the resulting statistical uncertainties are significantly lower than the ones ob- 
tained in the procedure described in Chapter [TOl This reduction is due to the smaller number of 
iterations used. The obtained central values from the two different approaches were compatible 
within the respective statistical uncertainties. 

The systematic uncertainties in the individual channels were of the same order for both ap- 
proaches for the inclusive unfolding. Larger deviations were observed for the simultaneous 
unfolding in \yt\ — \yt\ and A%, where the statistical component in the evaluation of the system- 
atic uncertainties was larger. These statistical fluctuations are expected to be reduced by using 
a smaller amount of iterations in the unfolding procedure. This effect can be seen in particular 
for the (dominant) contributions involving the replacement of the response matrix with samples 
with a smaller number of events. For most of these cases, the obtained systematic uncertainties 
are significantly lower than the ones obtained in the procedure described in Chapter [TOl 

Consequently, the unfolding procedure used in the main part of this analysis, despite being 
more conservative, is expected to be significantly more stable, in particular due to the reduced 
model dependency and the requirement of a convergent unfolding process. 

In addition, the conservative systematic uncertainties assumed for the QCD multijet normal- 
isation in the muon+jets channel have been replaced by a combined normalisation and shape 
systematic uncertainty as described in Chapter 16.3. 1[ All other parameters of the analysis have 
not been changed with respect to the nominal procedure described in Chapter 110.31 This ap- 
proach yielded uncertainties on the measurement which were significantly lower than for the 
assumption of a 100 % normalisation uncertainty, as expected. Uncertainties of 0.0007 for the 
inclusive unfolding, and 0.002 (M^j < 450 GeV) and 0.001 {M^^ > 450 GeV) for the simultaneous 
unfolding in \yt\ — \yt\ and Mjf were obtained (compared to 0.0011, 0.018 and 0.004, respectively, 
as shown in Chapter llOp . 
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E. b Tag Weighting Control Plots 



It was verified that the direct apphcation of b tag efficiencies to distributions obtained without 
the requirement of at least one b tagged jet yields effective distributions comparable to the 
ones obtained with the full event selection. A comparison was performed for the distributions 
of \yt\ — lutl both for the background contributions (c.f. Figure [ETT]) only and for the data 
distribution after background subtraction (c.f. Figure lE?2]) . The distributions are in excellent 
agreement within the shown statistical uncertainties. 
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Comparison of the distribution for \yt\ — \yt\ for the Monte Carlo background 
contribution for the full event selection (blue) and without b tag requirement, but 
with b tag efhciencies applied directly in a reweighting approach (red, dashed). 
The bin-by-bin statistical uncertainty is reduced by up to 25%. 
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Comparison of the distribution for \yt \ — \yt\ for data after background subtrac- 
tion for the full event selection (blue) and without b tag requirement, but with 
b tag efficiencies applied directly in a reweighting approach to the subtracted 
background contribution (red, dashed). 
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